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SUMMARY 


The  present  project  deals  with  the  fuel  spray  modelling  in 
multifuel  engines  using  the  droplet  behaviour  approach.  The 
analysis  given  in  The  First  Annual  Report  showed  that  the  spray 
is  affected  by  air  flow  field,  heat  transfer,  fuel  atomization, 
droplet  behaviour,  fuel  evaporation  and  chemical  kinetics  . 

The  present  work  deals  with  the  flow  field, spray  modelling  and 
fuel  atomization  inside  multifuel  engine  swirl  chambers  . 

Prediction  of  the  flow  field  is  based  on  the  solution  of 
the  finite  difference  from  the  governing  differential  equations 
for  the  transport  of  mass  ,  momentum  and  energy  .  The  turbu¬ 

lence  kinetic  energy  K  and  its  dissipation  rate  €  model  of  turbu¬ 
lence  has  been  used  .  A  special  technique  has  been  used  to 

satisfy  compressibility,  streamline  curvature  ,  non-isotropy  re¬ 
sulting  from  the  operating  conditions  and  swirl  chamber  geometry, 
local  continuity  and  reduction  of  the  mass  source  at  each  grid 
node  to  zero  . 

Experiments  are  performed  in  a  set  up  that  allowed  variation 
of  swirl  chamber  geometry  and  operating  conditions  to  measure  the 
air  velocity  components  and  turbulence  fluctuating  velocities 
using  a  temperature  compensated  hot  wire  anemometer  .  A  care¬ 

ful  comparison  of  the  air  velocity  and  the  predicted  values  over 
a  wide  range  of  conditions  has  been  included  to  establish  genera¬ 
lized  equations  for  the  coefficient  of  turbulence  .  Furthermore, 

results  for  the  effect  of  engine  speed  and  swirl  chamber  geometry 
on  the  flow  velocity  and  turbulence  intensity  are  obtained  and 
discussed  . 

A  new  technique  for  the  measurement  of  air  flow  velocity 
inside  reciprocating  engines  based  on  the  microphone  condenser 
has  been  established  .  Experiments  are  performed  on  a  calib¬ 

rating  set  up  to  clearify  the  effect  of  medium  density  on  the 
velocity  measurements  .  Comparison  with  the  results  obtained 

by  hot  wire  anemometer  confirms  the  accuracy  of  the  microphone 
condenser  anemometer  . 

Prediction  of  the  air  flow  field  with  the  suggested  model 
which  includes  the  derived  coefficient  of  turbulence  shows  a 
fairly  good  agreement  with  the  measurements  .  Reasonable 

values  of  the  turbulence  mixing  length  at  wide  range  of  conditions 
have  been  obtained  .  The  solution  of  the  fuel  trajectories  , 

size  and  temperature  histories  is  carried  to  evaluate  the  source 
terms  in  the  gas-phase  equation  .  Experiments  have  been  planed 

to  photograph  the  spray  for  various  conditions  to  validate  the 

model  . 

Analytical  and  experimental  investigation  of  the  atomization 
quality  of  multifuel  intermittent  sprays  has  been  carried  out  . 

The  effect  of  fuel  properties  and  operating  conditions  on  size 
distribution  and  mean  diameter  at  any  instant  in  the  vicinity 
of  the  nozzle  of  a  pintle  type  is  studied  using  an  automatic  slide 
sampling  technique  which  is  built  in  a  simulating  model  for  multi¬ 
fuel  engine  combustion  chamber  .  The  physical  properties  of 
multifuel  are  firstly  determined  and  expressions  for  viscosity  , 
surface  tension  and  density  are  developed  . 
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Droplet  size  distribution  in  number  fractions  has  been 
determined  at  various  conditions  ,  It  is  found  that  the 

distribution  is  affected  by  fuel  pressure  and  properties  and 
slightly  affected  by  medium  density  .  Unique  comulative 

distribution  curve  is  obtained  .  Dimensionless  mathema¬ 

tical  expressions  based  on  seven  of  the  known  statistical 
distribution  functions  relating  the  droplet  size  distribution 
and  mean  droplet  diameter  with  physical  properties  ,  nozzle 
geometry  and  working  conditions  are  developed  and  compared. 
Relations  are  also  determined  for  different  types  of  mean 
droplet  diameter  as  well  as  the  maximum  diameter  . 
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Variation  of  fuel  viscosity  with  temperature  rise  for 

clear  and  multifuel. 

17c 

Fig. 21 

Variation  of  fuel  density  with  temperature  rise  for 

clear  and  multifuel. 
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Fig. 22 

Droplet  size  distribution  of  clear  fuel  sprays  injected  173 
by  a  pintle  type  injector  with  injection  pressure  118  bar. 

Fig- 2  3 

Droplet  size  distribution  of  diesel  gasolene  mixtures. 
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Fig. 24 

Droplet  size  distribution  at  various  medium  densities. 
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Fig . 25 

Droplet  size  distribution  at  various  injection  pressures . l?o 

Fig . 26 

Droplet  size  distribution  for  diesel  fuel  sprays  inside 
a  cold  engine  model  at  various  injection  pressures. 
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Fig. 27 

Droplet  size  distribution  for  diesel  fuel  sprays  inside 
a  cold  engine  model  for  various  nozzle  diameters. 
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Fig. 28 

Variation  of  sauter  mean  diameter  with  the  quantity  of 

diesel  fuel  in  multifuel  at  various  injection  pressures. 
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Fig. 29 

Effect  of  injection  pressure  on  sauter  mean  diameter  for  182 
clear  and  multifuel. 

Fig . 30 

Effect  of  medium  density  on  Sauter  mean  diameter  for  clear  fuel 
at  injection  pressure 

Fig. 31 

Comparison  between  Sauter  mean  diameter  D^2  and  other 
mean  diameters. 
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Fig. 32 

Comulative  frequency  distribution  of  droplet  size  for 
multifuel  sprays  with  dimensionless  droplet  diameter. 
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Fig. 33 

Chi-square  distribution  curves  for  different 
freedom  against  experimental  values. 

degrees  of 
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Fig. 34 

Comparison  of  various  distribution  equations 

periraental  values. 

with  ex- 
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Fig . 35 

Graphical  representation  of  equation  44. 
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Fig. 36 

Comparison  of  calculated  and  experimental  values  of 

Sauter  mean  diameter. 
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Discharge  coefficient 
Constant  pressure  specific  heat 
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Coeffecients  of  the  turbulence  model 
Velocity  coefficient 
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Pressure  correction , bar 
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Gas  constant , kj/kg . k 

Radial  co-ordinate 
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Coefficients  of  linearized  source  term 

Relative  swirl  chamber  volume  ratio 

Temperature , K 

Ambient  temperature 
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Density  coefficient 
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,m/s 
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Ejected  air  velocity, m/s 
Friction  velocity  =/7^7c,m/s 
Velocity  vector  at  any  point 
Dimensionless  velocity  =u/ut 
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e  Dissipation  rate  ofturbulent  kinetic  energy ;  ,  rn/s3 

e  Compression  ratio 

0  Angular  distance  co-ordinate 

Turbulence  mixing  length* iu 
k  Von  Karman  constant 

A  Turbulence  model  constant 

u  Viscosity, kg/m. s 

v  Kinematic  yiscosity  ^i/p^ipVa 

p  Density,  kg/m3 

E  Summation 

a  Prandtl/Schnidt  nuroher 

ak  c  h  Prandtl/ Schmidt  number  for  k,  c  and  h  respect iyely 

t  Shear  stress  ,N/ra2 

p  Arbitrary  dependent  variable 

ij;  Crank  angle  degree 

u>  Vorticity 
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CHAPTER  1 


INTRODUCTION 


The  development  of  multifuel  engines  is  required  to  solve 
problems  arising  from  the  disproportion  between  consumption  and 
production  of  diesel  fuel  and  from  the  use  of  petroleum  products 
in  petrochemical  industries.  This  development  can  be  achieved 
by  improving  the  combustion  characteristics  of  various  types  of 
fuel'  either  clear  or  mixture.  The  improvement  of  combustion 
characteristics  can  be  achieved  by  shortening  the  ignition  delay 
period,  increasing  air  motion  and  turbulence.  Therefore,  the  air 
movement  in  the  combustion  chamber  of  diesel  engines  has  a  great 
influence  on  the  improvement  of  multifuel  engine. 

A  diesel  engine  with  swirl  combustion  chamber  was  found  to 
be  more  suitable  for  combustion  of  multifuel  due  to  the  induced 
highly  turbulent  flow  during  the  compression  stroke  and  the  in¬ 
jection  period.  The  induced  swirl  motion  of  air  carries  the  ato¬ 
mized  fuel  away  from  the  Jet  and  forms  a  nearly  homogeneous  mix¬ 
ture  with  the  start  of  combustion.  Thereafter  reversal  of  flow 
takes  place  and  the  gases  stream  out  of  the  swirl  chamber  into  the 
main  engine  cylinder  where  they  continue,  their  burning.  The  flow 
velocity  and  turbulence  in  the  swirl  chamber  are  affected  by  en¬ 
gine  speed,  swirl  chamber  volume  ratio  and  port  area  ratio. 

Prediction  and  measurement  of  flow  and  heat  transfer  in  a 
motored  diesel  engine  cylinder  have  been  made |5, 6 , 33|using  mathe¬ 
matical  modelling  of  turbulent  flow  field.  Few  studies |20-22|have 
been  made  for  the  prediction  and  measurement  of  flow  field  in  diesel 
engine  swirl  chamber.  However,  the  determination  of 
the  effect  of  different  factors  affecting  it  have  not  received  de¬ 
tailed  study.  The  flow  inside  these  swirl  chambers  is  governed  by 
the  differential  conservation  equations  of  energy,  mass  and  momen¬ 
tum  for  turbulent,  flow.  These  governing  equations  are  obtained  by 
decomposing  the  velocity,  pressure,  density  and  enthalpy  into  a 
mean  and  random  value.  The  turbulence  diffusion  fluxes  which  appear 
in  the  governing  equations  are  modelled  using  Boussinesq  |g|  app¬ 
roach  in  terms  of  the  mean  flow  gradients  and  eddy  dif fusivities . 
These  eddy  dif fusivities  are  determined  by  solving  additional  dif¬ 
ferential  conservation  equations  of  the  kinetic  energy  of  turbulence 
K  and  its  dissipation  rate  e  .  This  approach  has  been  given  in 
detail  in  the  first  annual  report.  In  this  approach,  the  values  of 
the  turbulence  constants  have  been  considered  unchanged  which  may 
lead  to  missleading  value  of  the  mixing  length  during  the  compres¬ 
sion  stroke  and  inside  the  swirl  chamber.  Actually  streamline  wall 
curvature  may  increase  or  decrease  the  intensity  of  the  turbulence 
mixing  depending  upon  the  degree  of  the  wall  curvature,  and  it  can 
strongly  affect  the  skin  friction  and  the  heat  transfer  rates. 
Braishow  I9 \  has  showed  that  even  very  mild  streamline  curvature 
can  cause  significant  effects  on  the  turbulent  boundary  layer. 
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Suggestions  have  been  made  in  the  new  approach  to  modify  the  length 
scale  equation  for  rotating  flow  inside  swirl  chambers  by  making 
the  constant  of  the  t  -  equation  a  function  of  the  flux  Richardson 
number  based  on  the  extra  production  of  lateral  fluctuations  V2  due 
to  the  body  forces.  In  the  same  time  the  constant  Cp  in  the  tur¬ 
bulence  kinetic  energy  equation  has  been  determined  for  incompres¬ 
sible  isotropic  flow.  Experimental  results  show  that  the  use  of 
the  recommended  value  can  give  unreliable  results  especially  for 
turbulence  length  scale  as  it  depends  upon  the  local  turbulence 
conditions.  The  constant  Cp  is  multiplied  in  this  work  by  a  correc¬ 
tion  function  which  is  a  function  of  the  combustion  chamber  geometry 
and  engine  speed  to  modify  the  turbulence  equation  to  suit  nonisot¬ 
ropic  compressible  flow. 

The  aim  of  the  first  part  of  the  present  work  is  to  develop 
better  modelling  of  the  turbulence  flow  field  in  swirl  chambers. 
Moreover,  experimental  measurements  of  the  random  flow  velocity 
have  been  planned  to  satisfy  careful  comparison  over  a  wide  range 
of  regimes  with  the  predicted  values.  This  facilitates  the  deter¬ 
mination  of  reliable  expression  relating  the  constants  of  tur¬ 
bulence  model  with  the  engine  operating  conditions  and  swirl  cham¬ 
ber  geometrical  parameters.  An  experimental  set  up  has  been  built 
to  facilitate  measurements  of  the  air  velocity  components  inside 
the  swirl  chamber.  A  swirl  chamber  has  been  built  up  in  one  cy¬ 
linder  of  a  two  cylinder  4 -stroke  water  cooled  diesel  engine.  The 
other  cylinder  has  been  used  for  motoring  the  engine  to  the  requir¬ 
ed  regime.  The  air  velocity  components  have  been  measured  by  two 
measuring  techniques:  the  microphone-condenser  anemometer  and  tem¬ 
perature  compensated  hot-wire  anemometer.  A  calibrations  set  up 
has  been  built  up  to  calibrate  both  probes  since  their  output  sig¬ 
nals  depend  upon  air  pressure  and  temperature.  The  output  signals 
have  been  analysed  to  determine  the  mean  velocity,  random  velocity 
and  the  derived  statistical  values.  The  experimental  results  fa¬ 
cilitate  the  determination  of  the  turbulence  model  constants. 

The  solution  of  the  governing  equations  for  the  transport  of 
mass,  momentum,  energy,  kinetic  energy  of  turbulence  and  its  dis¬ 
sipation  rate  has  been  done  by  solving  their  finite-difference  form 
using  digital  computer.  Comparison  have  been  carried  out  with  the 
experimental  results. 

During  the  past  years  comprehensive  works  were  carried  out 
at  Cairo  University  about  spray  modelling.  Firstly  the  spray  cal¬ 
culation  model  was  based  on  the  solution  of  the  equation  of  motion 
of  droplets  and  air.  The  interaction  between  the  droplet  and  air 
was  considered  by  averaging  the  over  flow  processes  on  a  scale  smal¬ 
ler  than  the  droplet  size  and  employing  typical  correlations  for 
the  droplet  drag,  heat  and  mass  transfer.  The  results  of  this  work 
which  have  been  reported  in  several  publications  1 21 1  were  suffi¬ 
ciently  encouraging  to  extend  the  model  to  include  the  prediction 
of  air  motion  and  turbulence  interaction  between  fuel  and  air. 

Spray  modelling  based  on  the  combination  of  the  discrete  droplet 
model  with  a  multidimensional  gas  flow  prediction  method  is  the  aim 
of  the  second  part  of  the  present  work.  However,  this  model  is 
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still  awaiting  sufficiently  detailed  information  about  the  droplet 
size  distribution  and  droplet  yelocity  distribution. 

The  need  of  the  atomization  characteristics,  which  are  re¬ 
quired  for  spray  modelling,  has  been  a  motivating  factor  for  a' 
vast  amount  of  research.  It  is  the  purpose  of  the  third  part  of 
the  present  work  to  illustrate  and  discuss  the  problem  of  atomiza¬ 
tion  and  results  obtained  for  clear  and  multifuel. 

The  solution  of  the  flow  field  using  the  new  approach  and 
the  measurement  of  velocity  pattern  to  define  the  constants  of  the 
governing  equations  have  been  conducted  by  Prof.Elkotb,  M.M. ,  Dr. 
Abou-Elleil,  M.M.  and  Eng,  Salem,  S. 

The  solution  of  the  droplet  behaviour  has  been  conducted  by 
Prof.  Elkotb,  M.M.  and  Dr.  Abou-Elleil. 

The  atomization  characteristics  of  multi  fuel  and  the  effect 
of  various  parameters  have  been  determined  experimentally.  Gene¬ 
ralized  expressions  have  been  determined  for  the  calculation  of 
the  droplet  size  distribution  and  the  mean  droplet  diameter. 

Prof .Elkotb, M.M.  and  Eng.Abdalla,  M.A.  are  responsible  for  this 
part  of  the  report.  The  manuscript  of  this  annual  report  has  been 
written  by  Prof .Elkotb, M.M.  and  revised  by  Dr.Aly,  S.L.  Experimen¬ 
tal  set  ups  for  the  investigation  of  droplet  behaviour,  evaporation 
rate  and  chemical  rate  haye  been  built  up.  Their  details  together 
with  the  results  will  be  discussed  in  the  third  Annual  Report. 


CHAPTER  2 


MATHEMATICAL  MODELLING  OF  FLOW  IN  SWIRL 
CHAMBERS  OF  MULTIFUEL  ENGINES 


2.1.  Inroduction 

The  development  of  multifuel  engine  can  be  achieved  by  im¬ 
proving  the  burning  characteristics  of  fuel  by  shortening  its  igni¬ 
tion  delay.  This  could  be  improved  by  increasing  turbulence  inside 
the  combustion  chamber.  Diesel  engines  with  swirl  combustion  cham¬ 
bers  are  more  suitable  for  the  burning  of  different  types  of  fuel 
due  to  its  ability  to  evaporate  and  mix  heavy  fuels  as  well  as  light 
fuels  and  to  work  with  high  compression  ratios  |  22  j  .  Swirl 

chambers  are  also  characterized  by  its  ability  to  control  combus¬ 
tion  of  fuel  during  period  of  rapid  pressure  rise  and  to  obtain 
better  burning  characteristics .  The  study  of  pure  phenomena  in 
subclasses  will  help  the  solution  of  the  more  complex  combustion 
phenomenan  resulting  from  multifuel  combustion. 

Researches  during  the  last  few  years  show  that  the  air  move¬ 
ment  in  the  combustion  chamber  of  diesel  engine  in  general  and  in 
swirl  chambers  in  particular  has  great  influence  on  the  problem  of 
controlling  the  air-fuel  mixing  and  combustion  process.  Moreover, 
the  mixing  process  controles  combustion  which  inturn  controles  the 
rate  of  pressure  rise,  noise  and  pollutant  emissions.  |  5  j  .  It 

is  therefore  very  important  to  measure  the  flow  field  in  the  mo¬ 
tored  diesel  engine  and  to  study  the  effect  of  engine  operating  con¬ 
ditions  and  swirl  chamber  geometry  upon  the  flow  field  inside  it. 

Many  researches  on  the  prediction  and  few  measurements  of 
flow  have  been  carried  out  in  a  motored  diesel  engine  cylinder. 
Prediction  was  carried  out  using  mathematical  modelling  of  the  tur¬ 
bulent  flow  in  diesel  engine  cylinder  with  assumptions  which  keep 
the  results  far  from  the  actual  events  when  comparing  these  re¬ 
sults  with  the  measured  ones.  Prediction  of  flow  field  was  carried 
out  theoretically  by  numerical  solution  of  the  Navier  -  stokes 
equation  applied  to  a  turbulent  flow  motion  by  introducing  the  time- 
averaged  properties  together  with  statistical  correlations  involv¬ 
ing  fluctuating  velocities  and  temperatures  in  the  conservation 
equations.  Since  there  is  no  direct  way  of  knowing  the  magnitude 
of  these  terms,  approximation  or  modelling  of  their  effect  in  terms 
of  quantities  which  can  be  determined  was  used.  These  modelling 
terms  depend  upon  the  local  turbulence  conditions  especially  in  swirl 
chambers  which  have  not  received  detailed  study  till  now. 

2.2.  Review  of  Previous  Work 


Several  theoretical  studies  and  few  experimental  ones  have  been 
made  to  determine  the  flow  field  inside  engine  combustion  chambers 
and  engine  main  cylinder.  These  theoretical  studies  can  be  divided 
into  three  main  categories.  In  the  first  category  flow  systems  are 
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treated  as  being  in  a  quasi-steady  state  and  account  being  taken  of 
time  mean  average  properties  of  velocity,  pressure  temperature  and 
chemical  species  concentration.  ElkotbtlB|  used  this  theory  in 
the  calculation  of  the  flow  field  inside  the  cylinder  of  two-stroke 
engine  in  which  he  suggested  that  the  central  part  of  the  flow  under¬ 
went  a  solid  vortex  motion  and  the  part  near  cylinder  wall  underwent, 
a  nearly  free  vortex  motion.  He  also  used  this  theory  in  calculating 
the  flow  field  in  swirl  chambers  |20l  .  The  model  in  this  work  was 
based  on  the  equality  of  the  moment  of  momentum  resulting  from  the 
ejection  velocity  to  the  swirl  chamber  and  the  moment  of  momentum  of 
the  rotating  air  inside  the  swirl  chamber. This  model  can  define 
the  average  tangential  yelocity.  Ogasowara  lb)|  and  OZ  |  |  have  use 

this  theory  in  calculating  the  trajectory  of  'fuel  spray  inside  a 
diesel  engine  combustion  chamber.  However  in  many  cases  of  gas  motio 
tangential  separation  surfaces  arise  and  the  flow  on  either  side  of 
this  surface  is  considered  a  jet.  The  jet  may  be  moving  in  the  same 
direction  or  in  opposite  directions.  The  instability  of  these  tan¬ 
gential  separation  surfaces  causes  eddies  on  it,  which  move  in  dis¬ 
orderly  fashion  both  along  and  across  the  stream.  This  brings  about 
a  transverse  transfer  of  momentum,  heat  and  constituents  between 
neighbouring  jets  which  form  a  boundary  layer  of  finite  thickness 
with  a  continueous  distribution  of  yelocity,  temperature  and  species 
concentration.  At  very  low  Reynolds  numbers,  this  boundary  layer 
may  be  laminar.  In  these  works,  no  account  being  taken  of  the  tur¬ 
bulence  exchange  flow  which  affects  the  accuracy  of  the  flow  field 
calculations . 


At  the  present  time  several  theories  of  turbulence  are  known. 
Prandtl's  old  theory  of  free  turbulence  in  which  he  arrived  at  the 
well  known  formula  rXy  =  p  3u/3y  and  for  pure  turbulent  nonviscous 
motion  tXy  -=  -pu'v1  which  indicates  that  the  turbulence  shear 
stress  is  defined  by  the  lateral  transfer  of  the  longitudinal  com¬ 
ponent  of  mo  ntura.  He  assumed  also  that  the  transverse  fluctuating 
velocities  v*  are  proportional  to  the  longitudinal  fluctuating  ve¬ 
locity  u‘.  The  longitudinal  fluctuating  velocity  is  determined  in 
turn  from  the  fact  that  during  a  transverse  displacement  of  a  mol 
of  fluid  it  preserves  its  initial  velocity  u  until  its  transverse 
displacement  ceases,  and  the  change  over  the  distance  of  a  trans¬ 
verse  mole  length  equals  u'=  l  ! du/dy j,  where  i  is  called  Prandtl 
mixing  length.  Finally,  the  prandtl  formula  of  turbulent  shearing 
stress  was  written  in  the  following  form  : 
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xy 
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The  mixing  length  was  assumed  constant  in  the  transverse  direction 
due  to  the  absence  of  walls  and  proportion  to  a  characteristic  length 
in  the  longitudinal  direction  {l  —  cx} .  Taylor,  after  noticing  a 
contradiction  in  Prandtl's  theory  of  free  turbulence,  assumed  that 
the  turbulent  shearing  stress  in  a  stream  must  be  determined  by  the 
lateral  transfer  of  vorticity  and  not,  as  Prandtl  assumed,  by  momen¬ 
tum  transfer.  For  two  dimensional  flow,  the  mean  velocity  is  given 
by 

—  _  1  ,3  u  3  v 
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neglecting  J  v/J  x  for  a.\ially  symmetric  flow  the  mean  vorticity  is 
given  by 

.1  =  3  u/3  y 

If  the  mean  free  path  of  a  fluid  particle  in  a  turbulent  flow  equals 
during  a  transverse  transfer  of  a  particle  from  one  layer  to 
another,  then  vorticity  will  be  equal  to 

=  9  =  I  £ 

“U.  t  3  y  2  *T 

from  which  the  apparent  shearing  stress  is  obtained  in  the  follow¬ 
ing  form 

t  =  +  j  p  22  (3  u/3  y)  2 


where 


=  f/2 


Thus,  Taylor  removed  an  inconsistency  of  Prandtl* s  theory  resulting 
from  the  neglection  of  local  instantaneous  pressure  gradients  that 
substantially  influence  momentum  exchange  but  do  not  affect  vorticity 
transfer.  This  hypothesis  as  Prandtl  hypothesis  needs  one  experimen¬ 
tal  constant  for  the  mixing  length. 

Both  theories  gave  no  more  accurate  results  for  the  velocity 
profile  especially  in  cases  that  involved  temperature  profiles  and 
heat  transmission.  This  situation  forced  Prandtl  to  establish  a  new 
theory  of  free  turbulence  in  which  he  compared  the  general  expression 
for  turbulent  shear  xXy  =  -  pu'v'  with  the  formula  proposed  by 
Boussinesq  in  analogy  with  the  Newtonian  law  of  friction  in  laminar 
flow, 

x  =  pu.  3u/3y  ,  to  obtain  an  expression  for  vt  in  the 


3  u/3  y 


following  form  ; 


v  t  =  Z  3  u/3  y 

where  the  constancy  of  the  mixing  length  £  in  a  transverse  cross 
section  of  the  jet  causes  the  kinematic  viscosity  along  the  cross 
section  to  vary,  becoming  zero  at  the  boundaries  of  the  mixing  re¬ 
gion  when  3 u/3 y  =  0. 

In  the  second  category  of  mathematical  calculation  of  flow  field, 
the  new  Prandtl  theory  has  been  used.  Boni  |5|  used  this  theory  for 
the  calculation  of  the  velocity  profile  in  a  divided -chamber ,  stra¬ 
tified  charge  engine  using  a  correlation  developed  by  Ferri  (  )  for 

the  turbulent  viscosity.  Diwkar  |15|  has  also  calculated  the  velo¬ 
city  profile  in  the  cylinder  of  an  engine  during  the  four  strokes 
by  numerical  solution  of  the  two  dimensional  inviscid  flow  field 
equations.  Griffen  1 3 5 1  reported  similar  calculations  for  a  viscous 
flow  by  solving  the  complete  Navier-stokes  equations.  These  studies 
have  been  carried  out  at  low  Reynolds  number . 


Prediction  of  time  mean  average  velocity,  pressure,  temperature 
and  concentration  in  turbulent  combust  ion  system  using  stokes-Navier 
equation  can  be  made  provided  the  turbulent  flux  tensor  t,  the  tur¬ 
bulent  enthalpy  flux  vector  and  the  turbulent  chemical  species 

flux  vector  Jj  are  specified.  Turbulent  exchange  coefficients  p, fu,  T , 
are  generally  used.  The  main  difficulty  in  any  prediction  method  ^ 

is  lack  of  knowledge  of  these  turbulent  fluxes  and  exchange  coefficients. 
Without  assuming  isotropic  turbulence  the  p  , r.  associated  with 

the  turbulence  exchange  coefficient  require  individual  specification. 
Assuming  isotropic  flow  i.e  assuming  rh  and  F,  to  bear  constant  ratios 
to  the  isotropic  turbulent  viscosity,  which  gives  the  Prandtl  and 
Schmidt  numbers 

°h  =  p/Ih  '  °j  =  ^/rj 

will  give  one  unknown  p  in  the  governing  equation  system.  A  simple 
turbulence  model  for  p  would  then  allow  predictions  to  be  made.  In 
Prandtl1 s  mixing  length  as  mentioned  before  he  assumed  that  the  isot¬ 
ropic  viscosity  p  could  be  specified  in  terras  of  local  gradients  of 
mean  quantities  via  a  mixing  length  £  and  constant  mixing  length 
parameter  A  i.e.  £  =A  times  the  width  of  the  mixing  region  where 
A  =  0.007  to  0.13. 


In  the  third  category  of  flow  field  calculations,  more  recent 
work  on  calculating  turbulent  flows  has  postulated  that  the  turbulence 
may  be  adequately  described  by  the  turbulent  viscosity  p  =  pk^  L. 
Therefore,  two  differential  equations  are  proposed  for  the  kinetic 
energy  K  and  any  other  variable  Z  =  KL  where  L  is  the  turbulence 
length  scale.  These  two  equations  for  "k"  and  "z"  are  developed  by 
a  combination  of  physical  reasoning  and  intuitive  guess  work  and  is 
called  Energy-Length  model.  Another  two  differential  equations  has 
been  reported  by  Launder  and  spalding  l5ol  •  This  two  equation  model 
of  turbulence  is  based  on  the  solution  of  additional  partial  differen¬ 
tial  equations  for  turbulent  kinetic  energy  "k"  and  its  dissipation 
rate  "e"  as  follows  ; 


DK 
P  Dt 


3  ,Pt  3k  ,  ,3  u.  2  n 

J7  (57  37  1  +  vt{Ty]  "  CD  pe 


Dt 


3 _  A  3_e  . 

3y  sy1 


+  e  I  C  — -  -  C 

K  l3y  ^2 


p2K, 


which  contain  three  constants  CD  ,  and  C2-  They  suggested  that  the 
value  of  CD  near  the  wall  is  determined  by  the  relation  : 

c£  =  x/pk 


and  the  other  constants  from  these  relations 
C2 

C1  “ 

D 

=  2  C, 


A2/ (a 


C*  ) 
e  1 


and  on  the  base  of  a  wide  range  of  experimental  data  for  an 


incompressible  flow  they  reported  the  following  recommen¬ 
ded  values  ; 

CD  =  0.08  ,  Oj  =  1.44  ,  c2  =  0.18,  '  =  0.4 

Many  researchers  132 ,  i  1 ,  3  1  ,  28  ,  l  ^  rely  on  a  two  equation  model  of  tur¬ 
bulence  as  reported  by  Jones  and  Launder  |45  I  which  use  new  values 
for  the  turbulence  model  constants  as  ; 

CD  =  1.0  ,  Cj  =  1.44  and  C2  =  1.92 

They  studied  the  laminar  and  turbulent  flow  fields  in  internal  com¬ 
bustion  engines  using  a  computer  code  called  RPM  (Reciprocating 
Piston  Motion) .  They  compared  their  results  with  experimental  results 
obtained  by  others  _.:d  showed  qualitative  agreement  during  some  re¬ 
gimes.  No  information  has  been  reported  about  the  effect  of  operating 
conditions  and  combustion  chamber  geometry  upon  the  results. 

Ramos  et  al  |67  |  used  the  two  equation  k-traodel  of  turbulence 
with  the  same  constants  to  calculate  the  flow  field  in  a  motored 
axisymmetric  reciprocating  internal  combustion  engine  with  a  centrally 
located  valve  which  opens  and  closes  instantaneously.  They  included 
in  their  study  the  compressibility  effects  upon  turbulence.  No  com¬ 
parison  with  experimental  data  had  been  reported  to  validate  the  model. 
Elkotb  et  al  1 2 2 \  have  used  also  the  k-e  model  of  turbulence  taking 
into  consideration  the  compressibility  effect  to  calculate  the  flow 
field  rotating  .inside  a  swirl  chamber  of  a  diesel  engine.  Compari¬ 
son  has  been  given  with  experimental  -esults  which  showed  that  rotating 
flow  has  a  strong  effect  on  turbulence. 

Various  extensions  of  this  theory  to  flows  with  swirl  have  been 
proposed.  This  has  been  done  by  linking  the  r-0  shear  with  the  r-x 
viscosity  of  the  axial  equation  and  allowing  for  the  nonisotropy  of 
the  viscosity  by  taking  it  proportional  to  the  second  invariant  of 
the  mean  flow  of  deformation  tensor  as  follows  ; 
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The  factor  (1  +  XS  sx)  accounts  for  the  change  in  the  length  scale  due 
to  swirl  and  is  analogous  to  the  Monin-Oboukhov  formula 

a  =  HQ  (1  -  6^) 


where , 


(2w7r2) 


>  (fw) 

Richardson  number  R.  =  - - — ~  — - 

0  u/J  r )  2  -t  i  r  K-^  [vi/ r) 


$  is  an  adjustable  parameter. 

The  above  formula  is  suggested  to  consider  the  effect  of  streamline 
curvature  and  centriputal  acceleration  on  miking  length.  It  is  inten¬ 
ded  in  this  work  to  modify  the  model  of  turbulence  to  include  the 
nonisotropy  and  stream  line  curvature  effects  for  the  calculation  of 
the  flow  field  in  swirl  chambers. 

2.3.  Plan  of  the  Present  Work 


It  is  clear  from  the  previous  discussion  that  the  accuracy  of 
predictions  is  affected  by  the  isotropy  assumption  for  swirling  flow 
and  compressibility  of  medium  as  the  case  of  swirl  chambers  of  diesel 
engines . 

Recent  studies  of  the  flow  prediction  inside  the  cylinder  of 
internal  combustion  engines  were  carried  out  at  low  Reynolds  number 
and  weak  swirling  flow.  Mathematical  models  have  been  used  to  pre¬ 
dict  the  properties  of  turbulence  flow  motion  using  fine  grids  to 
render  truncation  error  insignificant.  No  consideration  of  the  phy¬ 
sical  properties  in  combustion  chamber  is  taken,  which  includes  the 
continuous  variation  of  the  medium  pressure,  temperature  and  heat 
transfer  during  flow.  Moreover,  the  combustion  chamber  configuration 
and  engine  operating  conditions  affect  the  streamline  curvature  and 
centripetal  acceleration  which  affect  directly  the  turbulence  mixing 
length.  Till  now  the  constants  in  the  current  k-e  model  of  turbulence 
have  not  been  determined  for  swirling  compressible  flows  or  extended 
to  portray  nonisotropy.  Thus  turbulence  model  await  more  extensive 
and  complete  experimental  data  before  Justifying  their  use  in  the  pre¬ 
diction  of  the  boundary  layer  type  of  flow  in  swirl  chamber  of  diesel 
engines . 

Since  recent  works  dispute  isotropy  assumption  for  swirling 
flows,  it  is  intended  in  the  present  work  to  predict  and  measure  the 
flow  field  in  motored  diesel  engine  swirl  chamber  and  to  compare  the 
predictions  and  measurements  carefully  over  a  wide  range  of  data  to 
give  reliable  information  about  the  constants  of  the  turbulence  model 
and  their  relation  with  the  combustion  chamber  geometry  and  engine 
operating  conditions. 

Predictions  are  based  on  the  solution  of  the  finite  difference 
form  of  the  governing  differential  equations  for  the  transport  of 
mass,  momentum  and  energy  as  well  as  the  solution  of  the  differential 
equations  of  both  kinetic  energy  of  turbulence  and  its  dissipation 
rate.  Suggestions  have  been  made  to  modify  the  k-e  model  of  turbulence 
to  include  compressibility  of  medium,  nonisotropy  for  rotating  flow, 
combustion  chamber  configuration  and  geometry,  and  engine  operating 
conditions. 

Experimental  measurements  of  the  random  air  velocity  are  carried 
out  by  two  different  probes  :  a  hot  wire  and  a  condenser-microphone 
anemometers.  An  experimental  set  up  has  been  built  to  facilitate  the 
measurements  of  air  velocity  components  inside  the  swirl  chamber  with 
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highest  possible  accuracy.  Special  calibration  set  up  has  also  been 
built  up  to  facilitate  the  calibration  of  both  hotwire  and  micro¬ 
phone  condenser  anemometers  and  to  investigate  the  effect  of  air 
flow  temperature  and  pressure  upon  their  response. 


2,4.  The  Governing  Conservation  Equations 

The  flow  is  governed  by  the  differential  conservation  equations 
of  mass,  momentum  and  energy.  The  flow  is  assumed  to  be  two-dimen- 
tional  for  cylindrical  swirl  chamber  with  small  width.  All  equations 
are  time  dependent  since  there  is  a  transient  turbulent  flow  change 
with  the  crank  angle  positions.  The  effect  of  turbulence  appears  in 
the  form  of  correlations  of  fluctuating  quantities,  which  are  often 
referred  to  as  turbulent  fluxes. 


The  equations  governing  the  flow  in  a  cylindrical  swirl  chamber 
are  obtained  by  decomposing  the  velocity,  pressure,  density  and  en¬ 
thalpy  into  a  mean  and  a  fluctuating  value  respectively.  Finally, 
it  is  assumed  that  the  turbulent  fluxes  of  momentum  (i.e.  Reynolds 
stresses)  and  thermal  energy  are  calculated  through  scalar  eddy  dif- 
fusivities  which  are  in  turn  determined  by  solving  additional  diffe¬ 
rential  conservation  equations  for  the  kinetic  energy  of  turbulence  K 
and  its  dissipation  rate  c. 

The  turbulent  momentum  and  energy  equations  differ  from  the 
laminar  versions  only  by  the  appearance  of  the  Reynolds  stresses  and 
the  temperature-velocity  correlations.  Both  have  been  modeled  using 
the  Boussinesq  ) 7 )  approach  in  terms  of  the  mean-flow  gradients. 

The  model  has  been  modified  to  introduce  compressibility ^effect  by 
analogy  with  laminar  stress  tensor.  The  equations  governing  the  tur¬ 
bulent  flow  field  in  polar  coordinates  can  be  written  as  follows  : 


i-  Continuity  Equation  : 

It  +  ?  97  (rpv)  +  r  3TT  (pul 
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ii-  Radial  Momentum  Equation  : 
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Tangential  Momentum  Equation  : 
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iv-  Energy  Equation  : 
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Since  no  combustion  has  been  considered,  the  enthalpy,  h,  is  the 
sensible  enthalpy  which  has  been  taken  equal  to  : 

h  =  Cp  T 

V"-  Turbulence  Kinetic  Energy  Equation  : 

The  turbulence  kinetic  energy  equation  is  a  modified  form  of 
the  turbulence  kinetic  energy  equation  used  by,  Launder  and  Spalding 
(12)  and  (13)  for  incompressible  flows.  The  introduction  of  compress¬ 
ibility  effects  has  been  done  by  considering  terms  containing  the 
divergence  term,  V.if.  The  equation  for  the  turbulence  Kinetic  energy 
is  obtained  by  multiplying  the  instantaneous  momentum  equation  by  the 
instantaneous  velocity  and  subtracting  the  product  of  the  mean  ve¬ 
locity  times  the  mean-momentum  equations  and  then  taking  the  average 
of  this  result.  The  final  result  is  obtained  in  the  following  form  : 
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where; 

cr  is  the  turbulent  prandtl /schraidt  Number  (-1.0) 

K 

The  term  before  the  last  one  on  the  right-hand  side  accounts 
for  the  effect  of  compressibility  in  the  K-equation,  namely  : 
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The  last  terra,  CD  p«.  ,  contains  the  factor  f  which  account  for 

the  effect  of  compressibility  mixing  length  scale  resulting  from  the 
change  of  the  engine  operating  conditions  and  combustion  chamber  geomet¬ 
ry.  The  constant  CD  was  suggested  by  Launder,  and  spalding  (12) for 
turbulence  in  the  vicinity  of  a  wall,  in  which  convection  and  dif¬ 
fusion  of  energy  are  nearly  always  negligible;  so  there  is  a  balance 
between  the  production  and  dissipation  energy  and  obtained  the  re¬ 
lation  : 
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They  recommended  a  value  equal  to  0.08  for  the  constant  CD  as  a  re¬ 
sult  of  wide  range  comparison  between  experimental  results  obtained 
for  incompressible  flow  and  prediction  calculations.  Experimental 
results  carried  out  on  reciprocating  engines  show  that  the  use  of 
this  value  can  give  unreliable  results  especially  for  turbulence 
length  scale  as  it  depends  upon  the  local  turbulence  condition.  A 
correction  of  this  kind  is  necessary  if  the  mean  velocity  is  to  be 
predicted  correctly  with  the  k—  e  model. 


vi-  Dissipation  Rate  of  Kinetic  Energy  Equation  < 

The  equation  for  the  dissipation  rate  of  the  turbulent  kinetic 
energy  in  the  modified  form  can  be  written  as  follows  ; 
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a£  is  the  turbulent  prandtl/Schraidt  Number  (-  1.3).  The  first 

terra  on  the  right  hand  side  in  the  e  -equation  contains  the  com¬ 
pressibility  effect  and  the  part  (I-C3  R^)  accounting  for  the  effect 
of  streamline  curvature  and  centripetal  acceleration  on  mixing  length 
and  allowing  for  the  nonisotropy  of  the  flow.  This  term  was  also 
suggested  by  Lilly  '154!  and  Launder  1 5 2 1  in  their  work.  Moreover, 
the  second  term  on  the  right  handside  in  the  e  —equation  contains 
the  factor  of  accounting  for  the  change  of  engine  operating  condi¬ 
tions  and  combustion  chamber  geometry.  The  constant  C3  can  be  de¬ 
termined  by  optimization  considering  data  for  a  wide  range  of 
flow  in  engine. 


2.5  ESTIMATE  OF  EXCHANGE  COEFFICIENTS 

i-  Estimate  of  C2  from  decay  of  turbulence  behined  grid  s 

Behind  the  grid  the  decay  of  turbulence  is  function  only  of 
the  streamwise  coordinate  x  ; 
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Hence  ; 


C2  =  CD(ra  -  n,/ 2) 


where  ra  and  n  are  the  exponents  of  K  and  mixing  length  *  in 
the  definition  of  e ,  so  for  : - 


e  =  K3/2/Jt 


n  =  -1) 


•  '  C2  2  CD 

where  CD  is  determined  experimentally  from  the  relation 
CD  =  (typk)2 

ii-  Estimate  of  from  near'^wall  turbulence  s 

As  suggested  before  for  turbulent  flows  near  walls  t 

x/  pk  =  CD^“  (7) 

Launder  and  spalding)51\  raeasurments  near  a  wall  indicated 
that  the  mean  velocity  gradient  yaries  with  distance  from  the  wall 
as  follows  : 


=  4VP  /  Ky 


where  k  -  constant  determined  from  their  experimental  data  and 
equal  to  0.4  . 

Substituting  equations  (71  and  (8)  in  equation  (6)  of  the 
dissipation  rate  of  turbulence  kinetic  energy  and  simplifying  - 
the  following  relation  can  be  obtained  c 
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2.6  GENERAL  GOVERNING  EQUATION 


Within  the  above  frame  work,  the  complete  equation  set  may  be 
compactly  represented  in  terms  of  a  single  general  equation  for  an 
arbitrary  dependent  variable  $  in  the  following  form 


3-c^>+  ?  |r(rpv^+  ?  |^(pU^“  ?  fr (rr<j,tr1'  ^2  h(rof$=s* 


where :  <j>  :  stands  for  radial  velocity  vr  tangential  velocity  uf 

total  enthalpy  h,  K  and  c. 

t  :  time 

p  ;  density 

F,  :  the  effective  diffusivity  coefficient. 

<  source/sink  terms  for  variable  <j> . 


Equation  (9)  represents  also  the  continuity  equation  by  replacing 
ij>  by  1 . 


The  definition  of  r  .  and  S*  are  given  in  table  (1) .  The  new 
quantities  appearing  in  this  table  include  the  effective  viscosity 
U  ff  which  is  defined  as  : 


1Jeff  =  u  +  pt 

and  the  turbulent  viscosity  is  connected  with  k  and  the  turbulence 
energy  dissipation  rate  e  by  the  following  algebriac  equation  ; 


Pt  =  CA*  p  k2 
where  Cp  ~  0. Q  9  f 

This  allows  for  the  nonisotropy  of  the  viscosity  in  swirling  flows 
for  the  various  combustion  chamber  geometry  and  engine  operating 
conditions . 


Table  1  Defination  of  r,  and  S 
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2.7 


Initial  Conditions 


At  the  inlet  plane  to  the  swirl  chamber  the  inlet  velocity, 
density  and  temperature  vary  with  crank  angle  degree.  They  are  specifi¬ 
ed  by  solving  an  ordinary  differential  equation  for  the  flow  from  the 
main  chamber  to  the  swirl  chamber,  or  vise  versa,  through  the  tangen¬ 
tial  port,  as  reported  by  Elkotb  |20|  and  can  be  summarized  as  follows: 
By  equating  the  rate  of  mass  transfer  obtained  from  the  continuity 
equation  which  has  the  following  form  : 
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with  that  obtained  from  the  energy  equation,  assuming  polytropic  com¬ 
pressible  flow,  which  has  the  following  form  : 
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an  ordinary  differential  equation  of  the  following  form  can  be  ob¬ 
tained  : 
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and  c  =  compression  ratio. 


as  =  V=A 
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instantaneous  volume. 


(16) 


The  solution  of  this  differential  equation  using  Rung-Kutt  4 
method  gives  the  instantaneous  value  of  the  parameters  X,  which  in  turn 
and  with  the  help  of  equations  (13)  and  (15)  the  dimensionless  main 
chamber  pressure  and  the  dimensionless  swirl  chamber  pressure  at  any 
crank  angle  position  can  be  determined. 

The  velocity  of  air  ejected  into  the  swirl  chamber  is  deter¬ 
mined  from  the  energy  equation  when  assuming  polytropic  compressible 
flow  which  has  the  following  form  ; 
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and  when  substituting  eqns.(13),  (14)  and  (15)  it  becomes 
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2 . 8  Numerical  Solution  of  Equations 


The  governing  set  of  parabolic  equations  (eqns  1-6)  which  are 
written  in  a  single  general  equation  (g)  are  solved  using  a  finite 
difference  procedure.  The  solution  allows  calculations  of  unsteady 
compressible,  laminar  and  turbulent  flows  with  heat  transfer.  In 
addition,  the  calculations  of  temperature  and  pressure  fields  have 
been  made.  So  the  swirl  chamber  is  overlaid  with  a  grid  and  circles 
Fig.  (1) ,  at  which  the  scalar  variables  h,k,e  and  p  are  stored;  while 
the  velocities  are  located  raid-way  between  the  pressures  which  deter¬ 
mine  them. 


Imaging  each  variable  stands  in  its  own  control  volume  then 
equation  (9)  can  be  integrated  over  this  control  volume  for  a  time 
increment  3t.  The  finite  difference  equations  for  the  governing 
equations  have  been  derived  and  can  be  compactly  represented  by  the 
following  single  equation  : 
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where , 


V  :  is  the  cell  volume 


The  fs  are  weighting  factors  whose  specification  is  determined 
from  the  choice  of  spatial  differencing  approximation.  This  in  turn 
influnces  both  accuracy  and  computational  stability.  It  follows  the 
following  specification  j- 
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pew  !  *s  ttie  local  factor  number,  defined  by  Runchal 
as  follows, 
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and  similar  formulae  exists  for  the  remaining  f's. 

Finally  equation  (20)  can  be  written  in  the  follows  form  : 

%  +  Mp  -  y  +p  -  i  Ac  +  $  %°  *  su  1211 

where, 

■  't  "  i  Ac 

and  S  p  is  given  in  table  (  2)  ,  and  S  are  coefficients  of  linearized 
source  terms  obtained  by  integrating  sA  Ac  are  defined  as  the  in¬ 
fluence  coefficients  which  give  the  combined  effect  of  diffusion  and 
convection. 

Table  (2)  Definition  of  the  Quantity  S  for  the  individual 
Variables.  p 
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The  first  step  in  solving  the  finite  difference  equations  is 
to  solve  the  energy  equation  using  the  old  fields  of  n,u,v,p  and 
Op  1 3 1 )  ,  to  yield  the  new  enthalpy  field  (hn)  and  hence  the  new  tem¬ 
perature  field  (Tn) .  At  each  time  step,  the  solution  algorithm  proceeds 
as  follows  : 

★  v*» 

With  guessed  values  of ^pressure, P  ,  the  momentum  equations  are  solyed 
to  determine  V  and  u  n  .  In  general  these  velocities  will  not  satis¬ 
fy  the  continuity  equation  but  will  produce  a  net  mass  source  at 

each  grid  node.  The  pressure,  density,  and  velocities  are  then  correct¬ 
ed  so  as  to  reduce  the  mass  source  at  each  grid  node  to  zero  in  the 
following  manner  : 
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*  :  indicates  tjuossed  or  preliminary  value 

p'  :  pressure  correction 
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Dr  and  Du  are  evaluated  from  the  rclavant  momentum  equation  : 
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Now  equations  (23)  to  (251  are  substituted  incontinuity  eqn  then  ; 
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which  can  be  written  in  the  following  form 
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where;  subscripts  E,W,N,S,  and  e,w,n,s  are  neighbouring  nodes  and 
corresponding  nodes  mid-way  between  both  nodes. 

Equation  (27)  is  solved  to  yield  the  p'field  which  is  required 
to  correct  p,v  and  u  fields.  The  new  K  and  r  are  then  compared  from 
their  differential  equation  (21)- 


Equation  (27)  is  then  a  modified  form  of  the  pressure  correction 
equation,  developed  here,  qiving  a  test  convergence  and  reduces  the 
total  computational  time. 

The  conservation  and  local  pressure  correction  equations  are 
solved  by  using  a  Gauss-elimination  line  by-line  double  sweep  technique. 
While  equation  (21)  converges  very  fast,  equation  (27)  needs  more 
sweeps  and  sometimes  more  than  one  iteration  to  reddce  the  residual 
mass  errors  to  an  acceptable  level.  The  time  increment  <5t  which  is 
used  in  the  present  programme  is  2.2.2xlO~4  which  is  equivalent  to  2° 
crank  angle  degree. 

In  addition  to  the  local  pressure  correction  p’,  a  global 
pressure  correction  p'  is  computed  based  on  overall  continuity  balance 
and  is  added  to  the  existing  pressure  field  before  solving  the  pressure 
correction  equation  (27 1. 

The  flow  diagram  for  the  solution  is  given  below 
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CHAPTER  3 


EXPERIMENTAL  SET  UP  AND 
MEASURING  TECHNIQUE 


3.1.  INTRODUCTION 

Experimental  information  about  flow  inside  swirl  chambers  is 
limited  due  to  the  difficulties  involved  in  measuring  turbulent 
shear  flow.  Shear  layers  are  formed  when  two  streams  are  brought  to¬ 
gether  at  different  velocities.  The  two  streams  of  fluid  are  generat¬ 
ed  from  the  introduced  air  through  the  port  and  that  of  the  flowing 
air  inside  the  swirl  chamber.  The  resulting  velocity  gradients  cause 
shear  forces  to  arise.  At  high  Reynolds  number  the  forces  generated 
are  severe  enough  to  cause  turbulent  flow. 

Because  of  the  strong  current  interest  in  air  flow  inside  swirl 
combustion  chambers  this  work  is  directed  to  determine  the  effect  of 
working  conditions  and  combustion  chamber  parameter  on  it.  This  is 
required  for  the  validation  of  the  suggested  mathematical  model  which 
leads  to  better  calculation  of  the  fuel  spray  behaviour. 

Therefore,  an  experimental  set  up  has  been  built  up  to  facili¬ 
tate  measurements  of  the  air  velocity  components  inside  a  diesel  en¬ 
gine  swirl  chamber  with  highest  possible  accuracy.  Two  different 
methods  are  used  to  measure  the  turbulent  flowfield  inside  the  swirl 
chamber.  The  first  one  is  a  constant  temperature  hot-wire  anomometer 
(DISA  55M  system)  fitted  with  a  platinum  wire  probe.  The  second 
method  is  the  microphone  condenser  type  which  is  based  on  the  measure¬ 
ments  of  the  dynamic  pressure  of  the  flow.  As  the  flow  inside  the 
swirl  chamber  is  unsteady  and  measuring  instruments  are  affected  by 
the  medium  pressure,  signals  of  the  velocity  components  detected  by 
the  hot  wire  anemometer  and  microphone  condenser  as  well  as  the  gas 
pressure  variation  have  to  be  recorded  with  time.  It  was  necessary 
to  record  data  over  a  quite  sufficient  period  to  obtain  better  samp¬ 
ling  of  the  flow. 

A  special  calibration  set-up  has  been  built  to  determine  the 
effect  of  pressure  and  temperature  on  the  response  of  the  two  anemo¬ 
meters. 

3.2.  EXPERIMENTAL  SET-UP 

The  experimental  investigation  has  been  carried  out  on  a  2- 
cylinder,  4-stroke  water  cooled  medium  speed  diesel  engine>manufac- 
tured  at  the  Egyptian  Helwan  Factory.  The  data  of  this  engine  is 
given  in  Appendix  -2 . 

The  experimental  swirl  chamber  has  been  built  up  in  one  of  the 
engine  cylinder  while  the  other  cylinder  has  been  used  for  motoring 
the  engine  to  the  required  regime.  The  engine  is  directly  coupled 
with  hydraulic  dynamometer  for  speed  variation  and  equipped  with 
various  instruments  required  for  measuring  the  main  engine  parameters 
as  well  as  the  air  velocity  components  and  pressure.  The  test  cylin¬ 
der  head  has  been  modified  and  a  swirl  chamber  has  been  constructed 
and  attached  instead  of  the  original  chamber.  The  tested  swirl  chamber 
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tic 1 1  i  t  a  t  os  the  variation  of  the  swirl  chamber  volume  ratio  and  the 
port  area  ratio.  The  intake  and  exhaust  manifolds  of  the  test  cylin¬ 
der  are  separated  from  those  of  the  second  motor incj  cylinder. 

The  following  requirments  are  considered  during  the  desiyn  of 
the  set  up  : 

1-  Possibility  of  changing  the  engine  speed. 

2-  Possibility  of  changing  the  port  area  of  channel  connecting  the 
main  cylinder  with  the  swirl  chamber  and  its  inclination. 

3-  Possibility  of  changing  swirl  chamber  volume  ratio. 

4-  Control  of  variables  to  ensure  repeatability. 

5-  Data  recording  and  processing. 

A  schamatic  drawing  of  the  experimental  set  up  is  shown  in  Fig. 
(*’.)•  The  experimental  set  up  is  equipped  with  the  following  instru¬ 
ments  : 

i~  X-probe  hot  wire  anemometer, 

ii-  microphone  condenser, 

iii-  Piezoelectric  transducer, 

iv-  Two  magnetic  transducers, 

v-  ' 'sc.  i  I  loseope 

vi-  Magnetic  tape  recorder  for  recording  various  signals. 

3.2.1.  Test  Swirl  Chamber 

The  air  velocity  pattern  inside  the  swirl  chamber  is  affect¬ 
ed  by  the  entering  inlet  moment  of  momentum  of  the  air  1 2 4  |  ,  which 
is  affected  by  the  swirl  chamber  relative  volume,  the  connected  chan¬ 
nel  relative  area  and  its  inclination.  Figure  (  *  )  shows  a  scheme 
of  the  designed  swirl  chamber  which  satisfy  easy  variation  of  the 
relative  swirl  volume  ratio,  port  area  ratio  and  inclination  angle. 

A  cylindrical  swirl  chamber  with  swirl  chamber  volume  ratio 
of  0.772  and  diameter  of  48.5  cm,  has  been  constructed  and  built  up 
with  compression  ratio  of  17.  The  swirl  chamber  has  an  optically 
flat  fused  quartz  window  on  one  of  its  side.  It  is  also  fitted  with 
a  cooling  jacket  supplied  with  circulating  water  to  control  the  cham¬ 
ber  wall  temperature.  Two  oblique  washers  are  fitted  on  either  sides 
of  the  quartz  plate  to  avoid  stress  concentration.  The  quartz  plate 
was  fastened  by  a  self  locking  screwed  culler  to  facilitate  easy 
and  quick  assembly.  Connecting  ports  of  different  sizes  have  been 
constructed  in  a  conically  shapped  plugs  to  supply  air  to  the  swirl 
chamber  from  the  main  chamber1.  By  these  plugs,  it  is  possible  to 
adopt  connecting  ports  with  various  configurations.  Sealing  between 
the  plugs  and  their  seat  is  satisfied  by  conical  fitting.  The 
diameters  of  the  connecting  ports  have  been  changed  between  8  and  15mm 
to  satisfy  variation  of  moment  of  momentum. 

The  swirl  chamber  is  equipped  with  the  hot-wire  anemometer 
probe  and  the  microphone  -  condenser  probe  to  measure  the  air  ve¬ 
locity  components.  The  radial  movement  of  any  of  these  probes  is 
possible.  A  piezoelectric  pressure  transducer  is  also  fitted  to 
the  swirl  chairbcr  for  recording  the  pressure  variation  inside  the 
chamber. A  ferrous  -konstantan  thermo  couple  is  fixed  to  measure 
the  surface  temperature. 
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3.2.2.  Measuring  Instruments 

The  accurate  determination  of  experimental  results  to  va¬ 
lidate  the  predicted  ones  requires  sensitive  and  accurate  measur¬ 
ing  instruments. 

i-  Air  Velocity  Measurements 
a-  Hot  Wire  Anemometer 

A  temperature  compensated  hot-wire  anemometer  is  employed  to 
measure  both  the  mean  and  the  fluctuating  velocity  components.  The 
use  of  hot-wire  anemometer  is  recomended  due  to  its  small  dimensions, 
low  inertia  and  small  disturbance  created  by  the  instrument  itself. 

The  theory  of  the  hot-wire  anemometer  as  given  in  references  j  38,40,41,19| 
is  based  upon  the  heat  transfer  between  the  flow  and  a  constant  tem¬ 
perature  wire.  The  constant  -  temperature  hot-wire  anemometer  type 
DISA  55M  is  employed  using  a  platinum  wire  of  20u  in  diameter  to 
withstand  mechanical  stresses.  The  platinum  wire  is  able  to  with¬ 
stand  the  effect  of  dust  and  oil  particles  as  well  as  the  required 
higher  temperature  coefficient. 


The  variation  of  the  electric  current  passing  through  the  hot  - 
wire  due  to  air  velocity  variation  is  amplified  and  displayed  on  the 
screen  of  a  two  channel  cathod-ray  oscilloscope.  A  time  display  unit 
is  used  to  satisfy  the  x-motion  of  the  velocity  signal.  The  time 
display  unit,  which  is  fixed  on  the  crankshaft,  consists  of  two 
discs  and  two  magnetic  transducers.  The  first  disc,  which  is  made 
of  steel,  is  prepared  with  a  slot  to  define  the  TDC .  The  magnetic 
transducer  which  is  fixed  opposite  to  this  disc,  generates  electri¬ 
cal  current  corresponding  to  the  crank  motion.  The  slot  of  the  disc 
is  fixed  at  a  position  corresponding  to  the  TDC.  The  TDC  ,  is  ad¬ 
justed  within  +0.5  degree.  The  second  plastic  disc  whicn  is  pre¬ 
pared  with  steel  liners  on  the  rim  correspond  to  5, 10,.... 360  degree 
of  the  crank  angle,  together  with  the  magnetic  transducer  display¬ 
ing  on  the  second  channel  of  the  oscilloscope  the  time  base.  The 
signals  of  the  flow  velocity  and  time  are  photographed  by  a  movie 
camera. 

b-  Microphone-Condenser  Anemometer 

A  microphone  -  condenser  probe  is  developed  to  detect  both  the 
mean  and  fluctuating  values  of  velocity.  The  probe  receives  pres¬ 
sure  fluctuations  and  transfers  them  to  electrical  quantities  from 
which,  with  suitable  assumptions  and  simplifications,  mean  velocities 
and  velocity  fluctuations  u,  v,  u'  and  v*  can  be  derived. 

The  basic  electric  circuit  of  the  used  Microphone  condenser 
is  shown  in  Fig.  (  4  ) ,  in  which  the  air  pressure  affects  the  conden¬ 
ser  diaphragm  causing  change  of  its  capacity.  A  variable  electric 
current  output  proportional  to  the  acting  pressure  will  be  obtained. 
The  electric  circuit  of  the  Microphone  condenser  is  fed  with  voltage 
through  a  high  frequency  generator.  The  frequency  of  the  passinq 
current  is  chosen  not  less  than  10  times  the  frequency  of  the  cyclic- 
flow,  resulting  from  the  piston  motion. 

A  water  cooled  microphone  -  condenser  probe  was  constructed  as 
shown  in  Fig. (  5  )  with  two  different  probes;  one  for  the  tangential 
component  and  the  other  for  the  radial  one. 
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c-  Pressure  Measurement 


The  variation  of  pressure  inside  the  swirl  chamber  is  measur¬ 
ed  using  piezoelectric  transducer  and  a  charge  amplifier.  The  output 
electrical  signals  from  the  transducer  and  the  time  display  are  fed 
to  the  oscilloscope  at  the  same  time  with  the  recording  of  the  crank 
angle  degree.  The  transducer  has  been  dynamically  calibrated  directly 
on  the  engine  by  using  a  maximum  pressure  indicator  (disindicator 
instrument).  The  linear  response  of  the  transducer  facilitates  the 
calibration . 

3.3.  CALIBRATION  SET-UB  AND  RESULTS 

3.3.1.  Calibration  Set-Up 

The  flow  in  the  diesel  engine  is  unsteady  resulting  in  a 
change  in  velocity,  pressure  and  temperature.  This  may  affect  to  a 
certain  extent  the  response  of  hot  wire  and  microphone  condenser 
anemometers.  Therefor  a  special  calibration  set-up  has  been  cons¬ 
tructed  to  calibrate  both  the  hot-wire  and  microphone  -  condenser 
anemometers  with  the  change  of  pressure,  velocity,  and  temperature. 

The  general  layout  of  the  calibration  set-up  is  shown  in  Fig. 

(6  ) .  The  following  requirements  have  been  satisfied  during  the 

design  of  the  set-up  : 

1-  Possibility  of  changing  air  flow  rate  (air  velocity) . 

2-  Possibility  of  changing  air  pressure. 

3-  Possibility  of  changinn  air  temperature. 

4-  Possibility  of  changing  turbulence  level. 

Compressed  air  from  a  reciprocating  compressor  is  supplied  to 
the  test  section  through  a  pipe  of  5  cm  diameter.  Its  velocity  can 
be  measured  by  using  a  calibrated  orifice  meter,  which  is  fixed 
after  the  heater  at  a  distance  of  100  cm  from  its  output  section. 

An  electronic  anemometer  has  been  used  to  measure  the  pressure  across 
the  orifice  meter.  A  thermocouple  and  a  pressure  gauge  are  used  to 
measure  the  air  temperature  and  pressure.  The  anemometers  are  fixed 
at  a  distance  10  cm  from  the  orifice  meter.  A  wire  grid  is  mounted 
directly  before  the  orifice  meter  in  order  to  change  the  air  turbu¬ 
lence  level  in  the  test  section. 

3.3.2.  Calibration  of  Hot-Wire  Anemometer 

Before  starting  the  calibration  procedure  of  the  hot-wire 
anemometer,  it  has  been  rotated  to  define  the  normal  position  of  the 
flow  which  has  been  detected  by  the  maximum  voltage  reading.  This 
reading  corresponds  to  the  air  velocity  and  pressure.  Several  read¬ 
ings  have  been  taken  at  various  pressures  keeping  the  air  velocity 
constant  and  this  has  been  repeated  at  several  air  velocities.  Dur¬ 
ing  calibration  the  following  parameters  have  been  recorded. 

-  initial  volt  of  the  hot  wire  anemometer  VQ 

-  Mean  voltage  of  hot  wire  anemometer  V 

-  Pressure  drop  across  orifice  meter  AP 

from  which  the  mean  velocity  is  obtained 

by  the  following  relation  U  =  m/s 

-  Ambient  air  temperature  °C  tQ 

-  Air  pressure  bar  P 

-  RMS  of  random  voltage  of  hot  wire  v 

anemometer . 


9_0scilloscope  10.Angle  scale  1J.Pressure  gauge  12_Wire  grid  13 -Orifice  meter 


The  experimental  results  have  beer  illustrated  in  Fig. (7  ).  The  re¬ 
lation  between  pressure  times  the  velocity  as  a  inaction  of  the  hot 

wire  anemometer  response  can  be  written  in  the  following  form  ; 

PU  =  A (V  -  VQ)n 

where  A  and  n  -  constants  to  be  determined  from  experimental  results. 

Therefore,  the  experimental  results  are  illustrated  on  a  log-log  scale 
as  shown  in  Fig. (  8)  and  the  constants  have  been  determined  equal  to: 

A  =  32  and  n  =  1.8 

This  equation  is  valid  for  pressures  greater  than  1  bar.  Regarding 
the  turbulence  level  of  the  flowing  air  it  is  found  from  the  experi¬ 
mental  results  that  the  variation  of  the  air  pressure  and  temperature 
has  no  effect. 

Studying  the  effect  of  air  temperature  on  the  response  of  the 
hot  wire  anemometer  it  is  found  that  the  change  of  temperature  up 
to  the  maximum  temperature  of  air  at  the  end  of  compression  stroke 
cannot  affect  its  compensation  when  using  platinum  wire  and  compen¬ 
sating  ratio. 

3.3.3.  Calibration  of  Microphone  Condenser  Anemometer 

It  has  been  found  from  the  experimental  results  that  the 
response  of  the  microphone  -  condenser  anemometer  is  affected  by 
both  pressure  and  temperature. 

i-  Effect  of  pressure  and  temperature  on  mean  velocity  : 

At  atmospheric  pressure  the  microphone  -  condenser  is  tested 
at  different  air  velocities  for  a  temperature  range  of  28°C  to  106°C. 
The  change  of  mean  output  signal  volt  for  different  air  velocities 
is  illustrated  in  Fig. (  g  )  from  which  it  is  found  that  the  temperat¬ 
ure  variation  has  a  negligable  effect  on  the  mean  velocity  response 
of  the  microphone  -  condenser.  Keeping  the  air  temperature  constant 
at  28  °C,  the  microphone  -  condenser  has  been  calibrated  at  different 
air  velocities  and  pressures  ranging  from  1  to  35  bar.  A  sample  of 
the  results  for  a  pressure  1  baris  shown  in  Fig. (  9 )  which  can  be 
represented  by  the  following  relation  : 

n ' 

U  =  A'  (u) 

where ; 

U  :  air  flow  velocity  m/sec. 
u  :  mean  output  from  raic  (mv) 

A',n':  calibration  constants  for  Mic. 

The  constants  A' ,n'  have  been  determined  for  different  air 
pressures  and  the  results  are  illustrated  in  Fig.(  lo)»which  can  be 
represented  by  the  following  relations  : 

A'  =  1. 999-5. 2xlO-2(p)+  4.052xl0-4  (p2) 

n'  =  0.6  +  2 . 95xlO~3 (p)  -  4 .023xl0~5 (p2 ) 

where , 


re  calibration  curve 


Log  PU 


Log (V- Vo  ) 

Fig. (  8  ) Determination  of  hot-wire  calibration 
constants 
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for  microphone  condenser  pr  be  at 
ssure  and  different  temperatures 


Calibration  factor  A 


=  1.999-5.2  x  Iff  p  ♦  4  .052  x  10  P 


r 


p  :  the  absolute  pressure  in  bar. 

b-  Effect  of  Pressure  and  Temperature  on  fluctuating  Velocity 

Mean  values  of  pressure  and  velocity  could  be  measured  as  ex¬ 
plained  before.  For  fluctuating  velocity  it  was  found  that  both 
pressure  and  temperature  affect  the  response  of  the  Microphone 
condenser . 


The  principle  of  the  microphone-condenser  is  based  upon  the 
measuring  of  the  velocity  fluctuations  from  the  dynamic  pressure 
fluctuations  as  stated  by  lenze  |5 3 1  .  The  basic  equation  for  the 

calculation  of  the  velocity  fluctuations  from  their  dynamic  pressure 
fluctuations  can  be  written  at  the  inlet  to  the  probe  as  follows  : 


total 


=  P+P'  =  i  j  (p+p)  |(u+u * ) 2+y * 2+w * 2 


il 


+  P 


static 


Cl) 


where 

P  :  mean  pressure  in  bar 
P1:  'luctuating  pressure  in  bar 

y',v'and  w' fluctuating  velocity  components  in  polar  coordinates 
_  r,  6,z. 

u  :  mean  velocity  component  in  radial  direction, 
p  :  mean  density  in  Kg/ra^  ^ 

p's  fluctuating  density  in  Kg/m  . 

Equation  (1)  can  be  written  as  : 

p+p<  =  £  p |u2+2uu'+u,2+v'2+v'2  |+.j  p | u2+2uu  *+u ' 2+y ' 2+w' 2  j  +-PgtatC21 


Considering  the  fluctuating  velocity  u'  and  neglecting  the  other  two 
fluctuating  velocities  and  the  terra  1/2  p'u’^,  the  dynamic  pressure 
can  be  subdivided  into  the  mean  pressure  and  the  fluctuating  one  in 
the  following  forms  respectively  : 


P 


1 

J 


p  u2  +  _P 


stat . 


(3) 


and  in  the  following  form  : 

1  /— 2  1  -  2  1  — • 

P*  =  pu  u*  +2  pu  +  J  pul  +  puu' 


(5) 


Rearanging  the  fluctuating  pressure  as  follows 


p  i  _  p~£2  (H_l  +  “ 1 ~  + 

u 


2  p 


u  ’ ' 
2  u2 


+ 

u  p 


(6) 


Taking  the  root  mean  square  for  both  sides  of  equation  (6)  and  de- 
viding  the  results  by  equation  (4)  the  following  relation  can  be 
deduced 


ru  - 


2  +  a,2  .  TP 

+  2a+T 


2  +a 

?a+2 


(7) 


where , 

a  =  Tp/  Tu  (8) 

Tu  =  /  u 


(9) 


and 


\/  p  '  2  /  p  -  /  T 


(10) 


T 

P 

Equation  (7)  can  be  used  for  the  calculation  of  velocity  fluc¬ 
tuations  u'  ar  v’  from  the  corresponding  pressure  fluctuations  if  the 
value  of  a  is  determined. 

Since  a  =  T  /  T 

P  u 

\£2  V77  0  C  S  (11) 

T  u  y^T2  '  f  T 

then  the  constant  "a”  could  he  obtained  using  the  mean  velocity  and 
mean  temperature  which  cauld  be  easily  measured.  Calculating  the 
value  of  "a"  along  the  compression  and  expansion  strokes,  it  is 
found  that  its  value  ranges  between  0.01  and  0.2. 

From  the  measurments  of  pressure,  temperature  and  velocity, 
these  two  extreme  values  of  "a"  are  constructed  in  Fig. (  ll)  ,  and 
it  is  noticed  that  there  is  a  small  difference  between  these  two 
values,  which  indicates  that  the  change  in  temperature  has  a  slight¬ 
ly  small  effect  on  the  value  of  velocity  fluctuations. 

3.4.  EXPERIMENTAL  TECHNIQUE 

The  aim  of  the  present  work  is  to  determine  the  effect  of  work¬ 
ing  conditions  and  geometrical  parameters  of  the  swirl  chamber  in  a 
diesel  engine  upon  air  flow  pattern.  This  is  required  to  modify  the 
mathematical  model  to  correct  the  turbulence  mixing  length  due  to 
stream  line  curveture  and  flow  nonisotropy  resulting  from  swirl  cham¬ 
ber  geometry  and  operating  conditions. 

The  following  experimental  data  have  been  determined  at  various 
conditions  defined  by  swirl  chamber  volume  ratio,  port  area  ratio, 
and  engine  speed  : 

1-  Average  and  random  air  velocity  components  defined  by  the  hot¬ 
wire  and  microphone  condenser  anemometers. 

2-  Combustion  chamber  pressure  variation, 

3-  Maximum  swirl  chamber  pressure. 

4-  Crank  angle  degree  as  a  time  base, 

5-  Top  dead  centre. 

6-  Engine  speed. 

The  two  position  technique  defined  by  the  x-probe  hot-wire  an 
emometer  has  been  used  to  calculate  the  radial  and  tangential  com¬ 
ponents  of  the  instantanous  velocity  at  each  location.  One  of  the 
wires  is  located  in  a  plane  perpendicular  to  both  radial  and  tan¬ 
gential  directions,  as  shown  in  Fig.(  12. a).  In  this  position  the 
wire  response  has  a  superimposed  effect  of  both  radial  and  tangential 
velocity  components.  The  second  wire  is  located  perpendicular  to 
the  first  wire  such  that  it  will  be  only  normal  to  the  radial  v  - 
locity  component.  The  hot  wire  response  in  this  case  results  only 
from  the  radial  velocity  component.  The  later  wire  can  determine  the 
radial  mean  velocity  Vr  and  its  fluctuating  velocity  V^.  .  The 
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tangential  velocity  can  be  determined  by  subtracting  from  the  response 
of  the  first  wire  the  instantaneous  radial  velocity.  Thereafter  the 
mean  tangential  velocity  and  its  fluctuation  can  be  determined. 

The  same  technique,  which  is  the  two  position  shown  in  Fig. (  i3  ), 
has  been  used  by  the  microphone  condenser  to  determine  the  mean  tangen¬ 
tial  and  radial  velocities  and  their  fluctuations. 

The  experimental  work  has  been  carried  out  at  various  swirl 
chamber  diameters  and  tangential  part  diameters  satisfing  relative 
swirl  chamber  volume  ratios  and  relative  port  area  ratios  given  in 
Table  ( 3  )  • 


Table  (3)  Variation  of  swirl  chamber  ratio  and  port  area  ratio. 


Swirl  chamber  diameter  mm 

48.5 

45.0 

41.5 

38.0 

34.5 

Swirl  chamber  valume  ratio 

.772 

.64 

.  545 

.457 

.376 

Connected  channel  diameter  mm 

8 

10 

12 

14 

16 

Port  area  ratio 

.0041 

.0064 

.0092 

.0125 

.0164 

The  ergine  speed  has  been  changed  also  between  800  and  1500  RPM. 

One  of  the  previous  parameters  has  been  changed  while  the  other  para¬ 
meters  are  kept  constant  within  the  following  values. 


Compression  ratio 
Swirl  chamber  volume  ratio 
Tangential  port  area  ratio 
Engine  speed 


17,0 

0.772 

0.009 

1000  rpm 


3.5.  ERROR  ANALYSIS 


It  is  known  that  the  actual  value  of  any  measurements  can  not  ex¬ 
actly  be  obtained.  The  determination  of  the  possible  maximum  errors 
which  is  expected  during  the  readings  of  the  measured  quantity  is  re¬ 
ported  in  this  section. 

3.5.1.  Systematic  Errors 

These  kind  of  errors  are  created  from  the  imperfection  of 
measuring  methods,  lack  of  accuracy  and  personal  errors. 

Systematic  errors  can  be  calculated  as  follows  :- 
Let(G)  be  any  function  of  x^,X2»....  ,  xn 

G  =  G(x1,x2, . ,  xn)  (  I  2  ) 

where,  x^x2,...,xn  are  independent  variables 

G  may  be  any  one  of  the  measured  quantities. 

Then  the  error  in  any  one  of  the  measured  quantities  can  cal¬ 
culated  from  the  following  relation  ; 

; 

i=l  3x. 


dG 


(  »3  ) 
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applying  this  equation  to  the  main  measured  quantities,  their  error  can 
be  evaluated  : 


1-  Flow  velocity  :- 

The  pressure  drop  in  the  orifice  meter  of  calibration  set  up  has 
been  measured  by  water  manometer  of  accuracy  1% 


U  = 


Hence , 


6U 

U 


2gH 

0  • 
Kair 

1  6H 


(  14  ) 


=  +  0.005 


2  -7  300 

and  the  maximum  error  in  velocity  U  is  +  0.5%. 

2-  Error  in  hot-wire  reading 


PU  =  PU  (V,V  ) 
and  the  error  becomes 


6PU 

6V 

0 

> 

<o 

PU 

V 

Vo 

0.093 

,  0.093 

m  i_  _  _ ■ _ 

10.02 

7.68 

=  +  0.0214 

The  maximum  error  in  velocity  measurements  by  hot-wire  is  +  2.14%,  which 
eleminates  other  errors  resulting  from  cooling  effect,  cprrosion  tem¬ 
perature  compensation  and  unhomogeneous  distribution  on  wire. 


3-  Error  in  Microphone^condenser  reading 
Vs  -  U<Vvolt> 


6U  <5V  0.0032 

U  V  0.112 


=  +  0.0285 


The  maximum  error  in  velocity  measurements  by  microphone  condenser  is 
+  2.85%. 


3.5.2.  Accidental  Errors 


It  includes  the  errors  originated  from  accidental  unknown 
causes.  Such  errors  are  too  complex  in  nature  to  be  tracted.  It  can 
be  estimated  by  the  help  of  theory  of  probability. 


Let  x 


1 '  x2 


be  n  observations  of  a  variable  x  at  the  same 


regime  then  the  most  pro-Mble  value  of  x  is  given  by 


x  - 


x2  * 


r 

T. 

i=l 


(  ) 


Hence  the  mean  value  of  the  deviations  from  x  is 


e 


n 

+  >: 

i=  1 


Ai2/n 


1^  ) 


1 


where , 


-  x 


and  the  most  probable  error  which  oecurs  with  probability  50%  for  a 
single  observation  is  given  by  : 

n  _ 

e_,  =  +  0.67  Z  Ai  / {n-1 )  (  H  > 

P  i=l 

The  most  probable  error  of  the  measured  parameters  in  all  tests 
was  calculated  from  equation  (3.30)  for  10  observations  of  the  same  re 
gime.  It  was  found  that  the  probable  error  for  the  swirl  chamber  maxi 
mum  pressure  and  engine  speed  are  +  1.81%  and  +  0.6%  respectively. 


CHAPTER  4 


PRESENTATION  AND  DISCUSSION  OF  RESULTS 

This  work  locks  in  the  modelling  of  flowfield  inside  a  motored 
diesel  engine  swirl  chamber,  by  developing  a  mathematical  model  ca¬ 
pable  to  predict  the  air  velocity  inside  swirl  combustion  chamber. 

For  this  purpose  it  is  required  to  investigate  experimentally  the 
effect  of  various  geometrical  parameters  of  the  swirl  chamber  and  en¬ 
gine  working  conditions  on  the  flowfield  to  validate  the  computation¬ 
al  procedure  with  a  wide  range  of  reliable  information. 

This  chapter  includes  a  careful  comparison  of  the  air  flow  ve¬ 
locity  and  the  predicted  yalues  over  a  wide  range  of  experimental  data. 
Generalized  dimensionless  equations  for  the  coefficient  of  the  turbu¬ 
lence  model  which  depends  on  the  operating  conditions  of  engine  and 
swirl  chamber  geometry  are  derived.  Furthermore,  results  for  the 
effect  of  engine  speed  and  swirl  chamber  geometry  on  flow  velocity 
and  turbulence  intensity  are  predicted  and  discussed. 

4 . 1  EFFECT  OF  SWIRL  CHAMBER  GEOMETRY  ON  THE  TURBULENCE 
COEFFICIENT  CD 

The  value  of  the  turbulence  coefficient  Cp  at  various  radial 
distance  has  been  calculated  from  the  experimental  results  using  the 
following  equation  s 

CD  =  (t/!'K) 2  (  1) 

where 

T  =-■  pu  1 V 1 
K  -  \  Cu,2+v,2+w'2} 

The  measurements  have  been  carried  out  for  several  operating  conditions 
keeping  the  reference  values  constant  as  mentioned  in  Table  ( 4  )  and 

varying  only  one  parameter. 

Table  (  4  )  Reference  Operating  and  Geometrical 
Conditions 


Engine  speed 

Swirl  volume 

Ejection  port 

Measurement 

(r.p.m) 

ratio 

area  ratio 

point 

1000.0 

0.772 

0.0092 

At  T.D.C.  plane 

4  radius  12.2mm 

A  sample  of  the  calculation  results  of  the  turbulence  coefficient  CD, 
along  radius,  at  engine  speed  1000  r.p.m.  swirl  volume  ratio  0.772, 
ejection  port  area  ratio  0.0092,  plane  4,  and  top  dead  center,  is  il¬ 
lustrated  in  Figure  (  14  ) .  It  is  clear  from  the  figure  that  the  ex¬ 
perimental  results  fall  on  a  parabolic  curve,  which  can  be  represent¬ 
ed  with  adequate  accuracy  by  the  following  relation  : 

CD  =  A1  (r/Rs)nl  (  2) 

Plotting  the  experimental  results  on  a  log-log  graph  the  values  Oj.  A^ 
and  nj  have  been  calculated  and  equal  1.16  and  -0.08  respectively. 

The  equation  of  the  dimensionless  coefficient  CD  becomes  ; 


-  oj  - 


Dimensionless  radial  distance  r/Rs 


o5 


-0.08 

CD  =  1.16  <r/Rs)  (  3) 

It  is  clear  from  Fig.  (14  )  that  tne  coefficient  CD  has  different 
values  along  the  radial  distance  as  it  depends  upon  the  local  tur¬ 
bulence  condition  obtained  from  the  variation  of  the  ratio  between 
the  turbulence  shear  stress  rand  the  turbulence  generation  pK.  The 
coefficient  has  a  maximum  value  of  about  1.42  near  the  swirl  chamber 
centre  .at  which  the  shear  stress  r  is  higher  than  the  turbulence 
generation  pK  and  it  decreases  behind  the  swirl  chamber  wall  reach¬ 
ing  a  value  of  1.165  due  to  the  effect  of  chamber  wall  on  the  turbu¬ 
lence  generation  pK. 

Studying  the  effect  of  engine  speed  on  the  coefficient  Cq  for  the  same 
mentioned  parameters  at  a  distance  12.2  mm  from  the  swirl  chamber 
centre,  it  is  found  as  shown  in  Fig.  (15  )  that  the  results  fall  on  a 
parabolic  curve,  which  can  be  represented  by  the  following  relation: 

-0.88 

CD  =  1.2255  (N/Nq)  (  4) 

It  is  clear  from  the  figure  that  the  turbulence  shear  stress  r  is 
higher  than  the  turbulence  generation  at  low  speeds  than  at  high 
speeds.  This  indicates  that  the  turbulence  of  flow  inside  the  swirl 
chamber  increases  with  the  increase  of  engine  speed.  It  is  indicated 
also  that  the  engine  speed  afects  the  coefficient  CD  signif icantly 
while  Cq  along  the  radial  distance  has  a  smaller  variation. 

Studying  the  effect  of  ejection  port  area  ratio  (Aj/A)  upon  the  tur¬ 
bulence  coefficient  CD,  it  was  required  to  calculate  it  from  the 
experimental  results  obtained  by  hot-wire  and  microphone-  conden¬ 
ser  anemometers.  A  sample  of  results  are  shown  in  Fig. (16)  where 

is  plotted  versus  the  dimensionless  ejection  port  area  ratio. 
It  is  clear  from  the  figure  that  the  ejection  port  area  has  a  strong 
effect  upon  the  coefficient  Cq  and  that  it  increases  with  the  increase 
of  ejection  port  area.  This  shows  that  the  turbulence  shear  stress  t 
is  higher  than  the  turbulence  generation  pK  and  therefore  the  turbu¬ 
lence  of  the  flowfield  inside  the  swirl  chamber  increases  with  the 
increase  of  ejection  port  area.  Similar  to  the  effect  of  both  radial 
distance  ratio  and  engine  speed  ratio,  it  is  found  that  the  experimen¬ 
tal  results  tall  in  a  parabolic  curve  of  the  following  form  : 

0.704 

CD  =  33.359  (A.  /A)  (  5) 

Regarding  the  effect  of  swirl  chamber  volume  ratio  (V  /V_)  the  re¬ 
sults  of  calculation  of  the  turbulence  coefficient  are  plotted  in 
Fig.  (17  ) ,  from  which  it  is  clear  that  the  experimental  results  fall 
on  a  parabolic  curve  of  the  following  form  : 

0.229 

CD  =  1.303  (Vs/Vc)  (  6) 

It  is  noticed  from  the  figure  that  the  value  of  the  coefficient  Cp 
increases  with  the  increase  of  the  swirl  volume  ratio,  which  means 
that  using  low  swirl  chamber  volume  ratio  creates  turbulence  gene¬ 
ration  much  higher  than  when  using  high  swirl  chamber  volume  ratio. 

The  generalized  dimensionless  equation  of  the  effect  of  the  radi  i 
distance,  engine  speed,  ejection  port  area  ratio  and  swirl  chambi. . 
volume  ratio  on  the  turbulence  coefficient  can  be  obtained  by  plott¬ 
ing  (r/Rs) “°*°8 (N/N0) 88 (Aj/A) ° • 7°4 (Vs/Vc) °* 229  versus  CD  as  shown 


^'9  (  15  Variation  of  coefficier.  Cq  with  engine  speed- 
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uimensiomess  port  area  ratio  (Aj/A) 
i6  ) Variation  of  coefficient  Cg  with  ejection  port  area  ratio* 


imensionless  swirl  volume  ratio 


in  Fig. (18  ).  By  means  of  the  least  square  method,  it  is  found  that 
the  following  relation  fits  adequatly  the  combined  effect  of  the 
above  dimensionless  groups  on  CD  : 

-0.08  -0.88  -0.704  0.229 

CD  *  33.47  (r/Rs)  (N/No)  (A  -,/A)  (As/Ac) 

(7) 

Actually  the  coefficient  Cq  could  be  expressed  as  a  product  of  the  con¬ 
stant  value  of  Cqq  (Cdq=1.0)  and  a  factor  (F)  which  depends  on  the 
swirl  chamber  geometry  and  engine  speed  as  follows  : 


=  Cr 


8) 


where,  Cn  =  1.0 

F:  correction  factor  determined  from  experimental  results. 

Then  the  correction  factor  takes  the  following  form  : 

F=  3  3.47  (r/Rs)-0'08  (N/No)  "°*88  (A./A)0-704  (V g/V^  0,229 


9) 


The  coefficients  Cq  and  Cp  were  derived  for  near  wall  condition  while 
the  coefficient  C2  was  derived  from  the  decay  of  turbulence  behind  a 
grid.  The  coefficient  C2  can  be  written  in  the  following  form  : 


Cj  =  2  Cn  F 

o 

while  the  coefficient  can  be  written  as  follows  : 

l 

2 


=  <VCD  F  *  * 


/°e  (Cf 


F) 


(  10) 


(  ID 


The  coefficient  depends  on  the  local  turbulence  and  distance  and 
at  the  same  time  involves  the  other  two  coef f icients .  Therefore, 
the  coefficient  Cp  is  the  adjustable  coefficient  for  the  mathematical 
model  to  validate  the  experimental  data.  The  fine  adjustment  of  the 
model  could  be  obtained  by  the  coefficient  whcih  includes  the 
constant  A  as  it  depends  on  the  experimental  data  and  operating  con¬ 
ditions  |  50 1  . 

In  order  to  test  the  turbulence  model  incorporating  these  coefficients 
calculated  from  equations  {  7)  through  (  11)  a  wide  range  comparison 

between  predicted  and  measured  results  has  been  carried  out  in  the 
next  paragraph.  The  comparison  has  been  carried  out  also  with  the 
previously  predicted  results  { 2 2  j  using  the  following  turbulence  con¬ 
stants  : 


=  1.0 


C2  =  1.92 


and 


C1  =  1.44 


If.  2  COMPARISON  BETWEEN  EXPERIMENTAL  AND  PREDICTED  RESULTS 

The  K-  and  e  -  equations  incorporating  the  experimental  coeffi¬ 
cients  together  with  momentum  equation  in  radial  and  tangential  ’irec- 
tion  and  continuty  equations  are  solved  using  the  finite  difference 
procedure  mentioned  before.  The  computation  is  skirted  at  the  beginn¬ 
ing  of  the  compression  stroke.  The  instantaneous  inlet  velocities  to 
the  swirl  chamber  are  calculated  by  the  previous  computational  method 


equation  ( 


discussed  in  chapter  2.  Cc mput  ations  are  carried  cut  mainly  i  u  a 
diesel  enyine  haviny  a  swirl  chamber  of  48.5  mm  diameter  correspond¬ 
ing  to  a  swirl  volume  ratio  ot  0.772,  tanyenti.il  port  of  12  mm  dia¬ 
meter  corresponding  to  an  ejection  port  area  r  it io  of  0.0092,  enyine 
speed  of  1000  r.p.m.,  and  compression  ratio  of  :7.0.  The  tiro-  step 
of  computations  is  chosen  to  be  2°  crank  angle  interval. 

A  wide  range  comparison  of  the  present  predicted  results  during  the 
compression  and  expansion  strokes,  the  predicted  results  of  Re f .  2  2 
and  experimental  results  is  shown  in  Figures  ( 1 9  )  through  (  24  j 
for  tangential  velocity,  turbulence  intensity,  and  radial  velocity. 
The  constant  CD  has  been  -alculated  using  equation  (4.8)  and  the 
experimental  results  have  been  determined  by  the  hotwire  and  mic¬ 
rophone-condenser  anemometers . 

i-  Tangential  Velocity 

The  predicted  results  for  the  tangential  velocity  together  with 
the  experimental  results  are  determined  for  the  conditions  mentioned 
in  the  previous  paragraph  at  the  top  dead  centre  and  plane  4 .  The 
results  are  illustrated  in  Fig. (19  ),  which  shows  the  variation  of 
tangential  velocity  with  radius  at  the  end  of  compression  stroke. 
Figure  (20  )  shows  also  the  variation  of  the  tangential  velocity 
for  the  same  conditions  with  the  crank  angle  during  the  compression 
and  expansion  strokes.  It  is  clear  from  the  results  that  the  tangen¬ 
tial  velocity  along  the  radius  consists  of  two  zones,  the  solid  vor¬ 
tex  zone  at  the  centre  of  the  swirl  chamber,  and  the  free  vortex 
zone  near  the  swirl  chamber  wall.  A  quite  sufficient  agreement  is 
noticed  between  the  experimental  and  predicted  results  by  the  suggest 
ed  method  especially  near  the  swirl  chamber  wall.  Discrepancy  is 
noticed  near  the  centre  with  the  hot  wire  results  which  may  be  attri¬ 
buted  to  the  error  resulting  from  -the  determination  of  the  probe  in¬ 
clination  and  the  wire  length.  A  fairly  sufficient  agreement  is 
noted  in  the  free  vortex  zone  near  the  swirl  chamber  wall. 

Figure  (20  )  shows  the  variation  of  tangential  velocity  with  the  cran 
angle  during  compression  and  expansion  strokes.  A  quite  sufficient 
agreement  is  noticed  during  the  compression  stroke  between  the  ex¬ 
perimental  results  and  the  predicted  ones  obtained  by  the  suggested 
model  in  this  work. 


ii-  Turbulence  Intensity 

The  turbulence  intensity  has  been  calculated  from  the  experimen 
tal  and  predicted  results  by  the  following  relation  : 

1  =  /I 

and  the  results  at  various  radial  distances  and  ci  ink  angle  degress 
during  compression  and  expansion  strokes  are  illustrated  in  Figs. 

(21  )  and  (22),  respectively.  A  more  satisfactory  agreement  along 
the  radial  distance  is  noted  between  the  experimental  and  computed  re 
suits  by  the  suggested  model  than  that  obtained  by  fixed  turbulence 
constants.  It  is  noticed  from  Fig. (  21)  that  turbulence  intensity 
has  a  maximum  value  of  about  1.51  near  the  swirl  chamber  centre  bile 
a  minimum  turbulence  intensity  of  about  0.9  occurs  at  a  radial  dis¬ 
tance  ratio  of  r/Rs  =  0.6.  Figure  (  22'  shows  the  turbulence  in 
tensity  variation  with  the  crank  angle  at  r/R  -  0.503,  from  which 


Tangential  velocity 


Dimensionless  radial  distance  r/Rs 
Fig.  (  19  )Comparison  of  predicted  tangential  velocity  by  the 
suggested  and  previous  method  l  I  along  the  radius 

with  the  experimental  results  at  TDC  keeping  Vs/Vc  =0.772, 
Aj /A  =0.0092  and  N=1500r.  p.m.  constant. 
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Fig  (  20  )Comparison  of  predicted  tangential  velocity  by  the  suggest*’ 
and  previous  methodl  1  during  compression  and  expansion  strokes 
with  the  experimental  results  at  radial  distance  ratio  0.503  keeping 
Vs /Vc  =0.772  ,  A  j /A  =0  0092  and  N  =  1500  r.  p.m.  constant. 
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Dimensionless  radial  distance  r /Rs 
Fig.  (  21  Comparison  of  predicted  turbulence  intensity  ty  the 
suggested  and  previous  method  l  I  along  the  radius  with 

the  experimental  results  at  TDC  keeping  Vs/Vc  =0.772 ,  Aj /A  = 
0.0092  and  N  =  1500  r.  p.m.  constant- 
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Fig.  (  22  )  Predicted  and  experimental  variation  of  the  turbulence 
intensity  during  compression  and  expansion  strokes  at  radial 
distance  ratio  0.503  keeping  Vs/Vc  =0.772 ,  A j /A  =  0.0092  and  N  = 
1500  r.p.m.  constant. 
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it  is  noted  that  a  quite  good  agreement  is  obtained  during  the  com¬ 
pression  stroke  with  the  suggested  model.  hit  tie  deviation  is  not¬ 
iced  during  the  expansion  stroke  which  may  be  attributed  to  the  trun¬ 
cation  error  resulting  f r>  >.i  the  mathematical  iteration  solution. 

The  turbulence  intensity  has  a  very  low  value  at  the  beginning  of 
compression  stroke  and  increases  rapidly  during  the  compression 
stroke  reaching  a  maximum  value  of  1.5  at  about  loo°  AT DC .  There¬ 
after  it  decreases  again  during  the  expansion  stroke. 

i i i-  Radial  Veloci  by 

Figures  (23  )  and  (24  )  show  a  comparison  between  the  pre¬ 

dicted  and  measured  results  of  the  radial  velocity  at  various  ra¬ 
dial  distances  and  crank  angles  during  compression  and  expansion 
strokes  respectively.  Figure  (  23  )  shows  the  variation  of  radial 

velocity  at  various  radial  distances  at  planes  4  and  8  at  TDC,  good 
agreement  is  noticed  along  the  radial  distance  between  the  predictec 
and  measured  values.  There  is  no  discripancy  between  the  predicted 
results  computed  with  the  suggested  model  and  the  fixed  turbulence 
constants  model,  which  may  be  attributed  to  the  lower  value  of  the 
radial  velocity  component.  The  radial  velocity  lias  a  maximum  value 
of  about  8  m/s  at  about  2  mm  from  the  swirl  chamber  centre,  and  dec¬ 
reases  slowly  with  the  increase  of  the  radial  distance,  reaching 
zero  value  at  the  swirl  chamber  wall.  The  same  behaviour  of  the  radial 
velocity  is  noticed  at  plane  8  where  a  maximum  value  of  2.2  m/s  is 
obtained  near  the  centre.  This  indicates  that  the  presence  of  a 
source  at  the  centre  causes  decay  of  the  radial  velocity  as  it  travels 
towards  the  chamber  wall. 

Figure  (  24  )  show  the  variation  of  the  radial  velocity  during  com¬ 

pression  and  expansion  strokes.  A  quite  sufficient  agreement  is  not¬ 
iced  between  experimental  and  predicted  results.  The  change  of  the 
turbulence  coefficient  with  the  operating  condition  has  a  negligble 
effect  due  to  the  small  value  of  the  radial  velocity.  The  radial  ve¬ 
locity  has  a  low  value  at  the  beginning  of  compression  stroke,  and 
increases  rapidly  with  the  motion  of  the  piston  toward  the  TuC  reach¬ 
ing  a  m  iximum  value  of  hl.4m/s  at  38°  BTDC .  Thereafter  it  decreases 
with  a  higher  rate  than  the  rate  of  increase  during  compression  stroke 
reaching  a  zero  value  nearly  at  40°  ATDC. 

From  the  previous  disscussion  and  analysis  of  the  experimental  and 
predicted  results,  it  can  be  easily  detected  that  the  use  of  the  K-t. 
turbulence  model  with  coefficients  depending  on  the  operating  con¬ 
ditions  gives  better  results  and  satisfactory  agreements  with  the 
experimental  data.  Actually  the  use  of  the  K-t  turbulence  model  with 
the  empirical  values  of  these  constants  gives  not  only  deviations  from 
the  experimental  results  but  also  unlogic  values  for  the  turbulence 
mixing  length.  However  little  discripancy  is  noticed  with  the 

experimental  results  which  can  be  attributed  to  the  errors  result¬ 
ing  from  the  calibration  of  anemometers  and  nonisotropy  of  the  flow 
inside  the  swirl  chamber.  Besides  uncertainty  can  exist  as  a  result 
of  the  unaccurate  determination  of  the  velocity  vector  inclination 
angle.  The  measured  values  are  clearly  senstive  to  this  inclination 
angle . 
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Fig. (  24  )  Comparison  of  predicted  radial  velocity  by  the  suggested^and 

previous  method  f  I  during  compression  and  expansion  strokes  with 

the  experimental  results  at  radial  distance  ratio  0.503  keeping  Vs/Vc  = 
0.772  ,  A  j /A  =  0.0092  and  N=  1500  r.  p.m.  constant- 


4.3. 


EFFECT  OF  SWIRL  CHAMBER  PARAMETERS  AND  ENGINE  SPEED 
ON  THE  FLOW  FIELD 


The  flow  field  inside  the  swirl  chamber  is  affected  by  the  in¬ 
let  moment  of  momentum  of  the  ejecting  air  from  the  main  cylinder. 
Consequently , the  flow  field  is  affected  by  the  swirl  chamber  volume 
ratio,  the  relative  port  area,  and  the  engine  speed.  The  study  of 
the  effect  of  these  parameters  is  very  important  for  the  spray  motion 
investigation  and  for  the  validation  of  the  computational  model . 

4.3.1.  Flow  Field  During  Compression  and  Expansion  Strokes 

The  air  flow  field  inside  the  swirl  chamber  is  affected  by 
the  inlet  moment  of  momentum  of  the  aii  entering  it,  which  changes 
during  the  compression  and  expansion  strokes.  This  moment  of  momentum 
depends  on  the  pressure  ratio  between  the  main  cylinder  and  swirl 
chamber  (Pc/Ps)  and  consequently  on  ejected  air  velocity.  Figure 
(  25  )  shows  the  variation  of  the  pressure  ratio  and  the  ejected  ve¬ 

locity  versus  crank  angle  during  compression  and  expansion  strokes, 
from  which  it  is  clear  that  the  maximum  ejected  air  velor  ' ty  from 
the  main  cylinder  occurs  at  about  22°  BTDC .  Thereafter  ne  ejected 
velocity  decrea-os  so  long  as  the  crank  travels  towards  i_ne  TDC.  The 
flow  then  reversesits  direction  to  he  from  the  swirl  chamber  to  the 
main  cylinder  during  the  expansion  stroke.  Another  maximum  velocity 
is  noticed  at  about  20°  AT DC  during  expansion  stroke. 

The  computation  is  started  at  zero  time  at  which  the  compression 
stroke  starts.  The  computations  are  carried  out  at  2°  crank  angle 
intervals,  the  resulting  velocity  field  computed  at  various  crank 
angles  is  shown  in  Figs.  (  26  )  and  (27  ) ,  by  way  of  vectors  represen¬ 

ting  the  magnitude  and  direction  of  the  local  fluid  motion. 

In  time  sequence,  one  can  observe  the  features  of  the  flow  evolution. 

At  crank  angle  20°  ABDC ,  it  can  be  seen  that  the  flow  velocities  are 
low,  and  i'  almost  pure  radial  flow  directed  from  air  jet  to  the 
opposite  wall  surface.  After  some  time,  for  example  at  100°  ABDC, 
the  ejected  velocity  increases  and  the  built  up  pressure  increases 
causing  complete  rotational  flow  in  the  bulk  of  the  swirl  chamber. 
Before  the  initiation  of  the  complete  rotational  flow,  for  example 
at  60°,  100° ,  140°  and  160°  ABDC,  a  vortex  is  noticed  at  the  right  hand 
side  of  the  tangential  port  jet.  Before  the  end  of  the  compression 
stroke  the  radial  pressure  gradients  strengthen  to  the  point  of  induc¬ 
ing  strong  flow  recirculation  which  is  observed  at  20°  BTDC.  Dur¬ 
ing  expansion  stroke  the  jet  velocity  at  the  tangential  port  begins 
to  reverse  its  direction  from  swirl  chamber  to  the  main  cylinder.  At 
the  beginning  of  the  expansion  stroke  the  recirculating  flow  is  un¬ 
able  to  reverse  its  direction  to  flow  through  the  tangential  port 
inspite  of  the  absence  of  a  jet  flow  into  the  swirl  chamber,  as  shown 
in  Figure  (27  ) •  Before  the  complete  reverse  of  the  flow  a  vortex 

is  noticed  at  the  opposite  side  of  the  port,  and  low  flow  velocities 
can  be  noticed  during  the  expansion  stroke.  The  recirculating  flow 
begins  to  reverse  its  direction  at  60  degree  ATDC  and  at  140  degree 
AT DC  it  corapletly  changes  its  direction  to  rush  out  from  the  tangential 
port . 

The  predicted  and  measured  tangential  velocity  versus  crank  an>  le 
degree  during  compression  and  expansion  strokes  at  radial  distances 
of  2.2,  7.0,  12.2,  and  17.7mm  from  the  swirl  chamber  centre  are  shown 
in  Fig. (28  )•  The  solid  lines  represent  the  predicted  results  and 

the  back  signs  represent  the  experimental  results  obtained  by  hot-wire 


Crank  angle  degree  ij> 

Fig. (25  ) Predicted  ejection  velocity  and  pressure  ratio  during  compression 
and  expansion  strokes  for  A j /A  =0.0092 , Vs /Vc  =0.772  and  N  =  1000r.p.m. 


Fig  (26)  Velocity  vector  plot  during 
compression  stroke  for Vs/Vcs0772, 
NslOOOr  p.m  and  Aj/A  a  0.0092 


Crank 


Fig(  27  )  Veloci ty  vector  plot  during 
expansion  st r oke  fcr  Vs/Vc  =  0 772, 
NslOOOr.p.mand  Ai/As  0.0092 


Tangentiai  velocit  u  m^s 


Fig(  28  )Predicted  and  measured  tangential  velocity  during  compression 
and  expansion  strokes  in  plane  A  at  various  radii  kee  pi  n  g Vs/Vc=0.772; 
Aj  /  A  =  0.0092  and  N  =  1000  r.p.m  constant. 


anemometer . 


At  the  beginning  of  the  compression  stroke,  it  is  noticed  that  the 
tangential  velocity  component  has  a  negative  sign  due  to  the  anticlock¬ 
wise  flow  direction  which  continues  for  Severn  1  crank  angles.  As  the  crank 
angle  increases  during  the  compression  stroke,  the  tangential  velocity 
increases  in  the  clockwise  direction,  reaching  its  maximum  value  at 
about  2°  BTDC ,  and  then  a  slight  decrease  of  the  tangential  velocity 
occurs  at  about  10°  ATDC .  This  may  be  attributed  to  the  continuous 
decrease  of  the  ejected  velocity  till  zero  value  at  TDC  resulting  from 
the  decrease  of  the  pressure  gradient  across  the  port  as  shown  in  Fig. 

(  2rt  )  .  A  maximum  tangential  velocity  during  expansion  stroke  can  be 
noticed  in  Fig.(  28  )  at  about  22°  ATDC.  This  differs  from  the  angle 
at  which  the  maximum  ejected  velocity  occurs  during  the  expansion 
stroke.  The  same  behaviour  is  obtained  during  the  compression  stroke 
where  the  maximum  tangential  velocity  occurs  at  about  2°  BTDC  while 
the  maximum  ejected  air  flow  velocity  from  the  main  cylinder  to  the 
swirl  chamber  occurs  at  about  22°  BTDC.  The  maximum  value  of  the 
tangential  velocity  changes  from  plane  to  plane  through  the  swirl  chamber. 

The  predicted  and  measured  tangential  velocity  versus  the  radial  dis¬ 
tance  at  different  crank  angle  degrees  during  compression  and  expan¬ 
sion  strokes  are  shown  in  Fig.  {  29  ) .  The  tangential  velocity  can 
be  noticed  to  be  rather  small  along  the  radial  distance  at  20°  ABDC , 
and  *  ion  begins  to  increase  at  higher  crank  angles  during  compression 
stroke.  For  example  at  80°  BTDC  (  100°)  the  tangential  velocity 

has  a  negative  value  (anticlockwise  direction)  near  the  swirl  chamber 
centre  ana  a  closkwise  direction  near  the  swirl  chamber  wall  as  the 
complete  circulating  flow  is  not  reached  at  all  locations.  Before 
the  end  of  compression  stroke  the  tangential  velocity  distribution 
along  the  radial  distance  decomposes  to  a  solid  vortex  motion  and  a 
nearly  free  vortex  motion  near  the  swirl  chamber  wall.  The  free  vor¬ 
tex  motion  region  increases  as  the  crank  travels  towards  the  TDC  as 
shown  in  Fig. (  29  ) /  and  then  begins  to  decrease  as  the  crank  travels  to¬ 
wards  the  BDC  during  the  expansion  stroke  as  shown  at  80°  ATDC  and 
160°  ATDC. 

The  distribution  of  the  radial  velocity  versus  the  crank  angle  degree 
during  the  compression  and  expansion  strokes  at  different  radial  dis¬ 
tances  is  shown  in  Fig.  (30  ).  The  radial  velocity  distribution  ver¬ 

sus  radial  distance  for  different  crank  angles  is  shown  also  in  Fig. 

(  31  ) .  The  radial  velocity  has  a  low  value  at  the  beginning  of  com¬ 
pression  stroke,  and  increases  slightly  with  the  travel  of  the  crank 
towards  the  TDC  and  a  maximum  value  can  be  noticed  at  20°  BTDC  at  a 
radial  distance  of  2.2mm.  Thereafter,  a  rapid  decrease  occurs  dur¬ 
ing  the  expansion  stroke.  The  radial  velocity  diminishes  also  near 
the  swirl  chamber  wall  where  the  velocity  is  purely  tangential.  It 
is  noted  from  Fig.  (  31  )  that  a  source  is  created  at  the  centre  re¬ 

sulting  from  the  vorticies  existing  at  the  swirl  chamber  centre  with 
high  shear  stress. 

The  contours  of  the  dimensionless  turbulence  intensity  distribution, 
defined  by  the  ratio  ^2K/3/Cpm,  during  compression  stroke  are  shown 
in  Fig. (  32  )  It  is  clear  from  the  figure  that  the  turbulence  levels 

at  the  beginning  of  the  compression  stroke  are  very  low,  and  have  a 
maximum  value  of  0.082  in  the  centre  at  20°  ABDC.  During  the  ex¬ 
pression  stroke  the  inflow  to  the  swirl  chamber  becomes  much  stron¬ 
ger  and  begins  to  separate  ~ near  the  fcarugent i a  1  port  to  form  two  eddies 
on  either  sides  of  the  port.  For  example,  at  1000- ABDC~tTfe " "turbulence 


Fig(  3o)Predicted  and  measured  radial  velocity  during 
compression  and  expansion  strokes  in  planed  at  various 
radii,  keeping  Vs/Vc  =  0.772, Aj/A=00092  an d  N  =  1000  r-p.m. 


intensity  has  a  value  of  0.475  at  the  right  side  of  the  port  and 
0.54  at  the  left  side.  The  magnitude  of  these  turbulence  intensities 
generated  during  the  compression  stroke  due  to  the  hiqh  shear  stress 
and  pressure  gradient  on  both  sides  of  the  port  increases  so  long  as 
the  crank  travels  towards  the  TDC .  This  is  due  to  the  combined  ac¬ 
tion  of  fluid  inertia,  pressure  gradients  and  wali  shear.  The  normal 
stresses  produced  at  the  swirl  chamber  centre  causes  significant 
turbulence  intensity  increase  until  the  expansion  stroke  begins,  Fig. 

(  33  ) .  At  TDC  it  is  noticed  that  the  turbulence  generation  is  part¬ 

icularly  strong  in  the  shear  layers  along  the  port  and  has  a  value  of 
1.56  adjacent  to  the  entry  lip.  In  addition  there  are  significant 
turbulence  intensities  in-  the  swirl  chamber  which  have  been  built  up 
during  the  compression  stroke.  Near  the  chamber  walls  the  turbulence 
intensity  is  low  due  to  the  influe. ice  of  the  damping  factor.  The 
turbulence  generated  during  compression  stroke  reaches  a  maximum  value 
of  1.4  near  the  port  entry  at  TDC.  This  value  of  turbulence  intensity 
begins  to  decrease  during  the  expansion  stroke  and  has  a  value  of  1.23 
at  20°  AT DC .  The  vorticies  which  have  been  generated  during  the  com¬ 
pression  stroke  begin  to  disappear  in  the  expansion  stroke,  as  shown 
in  Figure  (33  ) ,  although  a  high  turbulence  intensity  exists  at  the 
centre  having  a  maximum  value  of  l.S  at  60°  ATDC .  Thereafter  it  begins 
to  diminishes  during  the  crank  travel  towards  BDC  where  small  vor¬ 
ticies  appear,  and  their  values  decay  due  to  tne  decrease  of  the  shear 
stress . 

Figure  (  34  )  shows  the  variation  of  the  turbulence  intensity  with  the 
crank  angle  degree  during  compression  and  expansion  strokes  at  plane 
4.  It  is  clear  from  the  figure  that  the  turbulence  intensity  is  very 
low  at  the  beginning  of  the  compression  stroke,  reaching  a  maximum 
value  at  about  20°  ATDC  where  the  outflow  from  the  swirl  chamber  to 
the  main  cylinder  is  maximum.  Thereafter  the  turbulence  intensity 
decreases  again  during  the  expansion  stroke.  High  value  of  turbu¬ 
lence  intensity  obtained  near  the  chamber  centre  (r=2.2  mm)  with  a 
value  of  about  1.84. 

Variation  of  turbulence  intensity  versus  radial  distance  in  plane  4  at 
different  crank  angles  is  shown  in  Fig. (  35  ) .  It  is  noticed  from 
the  figure  that  the  turbulence  intensity  has  a  low  values  along  the 
radial  distance  at  20°  ABDC  and  begin  to  increase  at  100°  ABDC .  A 
maximum  value  can  be  noticed  near  the  chamber  wall  and  centre  due  to 
high  shear  stresses  generated  in  these  regions.  At  TDC  (<p=  18 0°) 
a  maximum  turbulence  intensity  of  about  1.35  occurs  near  chamber  wall 
and  centre  and  a  minimum  value  of  about  0.74  occurs  at  a  radial  dis¬ 
tance  of  about  12  mm.  At  80°  ATDC  a  maximum  value  of  turbulence  in¬ 
tensity  of  1.55  is  noticed  near  the  chamber  wall  (r/Rs  =  0.74)  and 
a  maniraura  value  of  about  1.15  is  noticed  at  a  radial  distance  ratio 
of  about  0.25.  It  can  be  concluded  that  a  significant  change  of  the 
turbulence  intensities  along  the  radial  distance  specially  at  TDC, 
at  which  a  highly  difference  exists  between  the  maximum  values  occur- 
ing  at  chamber  wall  and  centre.  The  minimum  value  occurs  at,  n  arly 
the  midway  between  chamber  wall  and  centre. 

A >3. 2.  Effect  of  Port  Area  Rati o  on  the  Flow  Field 

The  air  flow  field  inside  the  swirl  chamber  is  affected  by 
the  port  area  ratio  as  defined  by  the  ratio  of  the  port  cross-s  tion- 
al  area  (A  j)  to  the  piston  area  (A) .  Therefore  the  investigation  of 
the  influence  of  the  port  area  ratio  on  the  flow  field  and  consequently 
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Fig(  34  JPredicted  and  measured  turbulence  intensity  during 
compression  and  expansion  strokes  in  p!ane4at  various  radii  keeping 
Vs/Vc  =  0  772  Aj/As  0.0092  and  N  s  1000  r- p.  m  constant . 
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Fig<  35  ) predicted  and  measured  turbulence  intensity  in 
plane4at  various  radii  ard  different  crank  angles  keeping 


Vs/ Vc  =  0.772,  Aj  /  A  =  0  0092  and  NslOOOr  p.  m.  constant. 


on  the  turbulence  model  coefficients  is  required  to  validate  the  pre¬ 
diction  model. 

The  calculation  shows  that  the  decrease  of  the  relative  port  area 
ratio  will  give  a  remarkable  increase  of  the  ejected  air  velocity 
through  the  port,  consequently  an  increase  of  the  flow  velocities  in¬ 
side  the  swirl  chamber  is  obtained. 

The  effect  of  the  relative  port  area  ratio  (A^/A)  on  the  liow  field 
velocities  inside  the  swir  i  chamber  is  shown  in  Fig.  (  36  )  .  A  re¬ 
markable  decrease  of  the  flow  velocities  can  be  not  iced  with  the  in¬ 
crease  of  the  port  area  ratio.  Tn  addition  the  Lw. .  eddies  which  had 
been  formed  during  the  compression  stroke  on  both  sides  of  the  port  dec- 
crease  with  the  increase  of  the  relative  port  area  ratio.  A  strong 
flow  recirculation  at  low  relative  port  area  ratio  is  shown  in  the 
figure.  This  may  be  attributed  to  the  existance  of  significant  ra¬ 
dial  pressure  gradients  when  using  low  port  area  ratio.  The  strong 
recirculating  flow  generates  a  large  tangential  velocity  as  shown  in 
Figs.  (  37  )  and  (  38  ) .  In  both  figures  it  is  clear  that  the  port 
area  has  no  effect  on  the  crank  angle  or  the  radius  at  which  the  maxi¬ 
mum  tanqential  velocity  occurs. 

Figure  (  37  )  shows  the  experimental  and  predicted  results  of  the 
effect  of  the  oort  area  ratio  (A; /A)  on  the  tangential  velocity  during 
compression  and  expansion  strokes  at  engine  speed  of  1000  r.p.m.  It 
is  clear  that  wit.  the  decrease  of  the  port  area  ratio  a  remarkable 
increase  in  tne  tangential  velocity  component  is  noticed,  reaching  a 
value  of  about  102  m/s  with  port  area  ratio  of  0.0041  at  2°  ATDC .  It 
then  decreases  again  during  expansion  stroke.  Another  maximum  velo¬ 
city  is  obtained  during  expansion  stroke  which  can  be  noticed  at 
about  26°  ATDC  and  has  a  maximum  value  of  about  97  m/s  at  port  area 
area  ratio  of  0.0041.  The  slight  decrease  of  the  tangential  velocity 
at  this  instances  is  obtained  due  to  the  reverse  of  the  flow  direction. 

Figure  (  38  )  shows  the  effect  of  the  port  area  ratio  on  the  distrib¬ 
ution  of  the  tangential  velocity  component  along  the  radial  distance 
of  plane  4.  It  is  noticed  that  both  the  solid  and  free  vortices  which 
are  built  up  at  the  end  of  compression  stroke  (TDC) ,  have  larger  va¬ 
lues  at  lower  port  area  ratio.  A  strong  free  vortex  near  the  chamber 
wall  is  noticed  at  port  area  ratio  of  0.0041. 

The  effect  of  the  relative  port  area  ratio  on  the  radial  velocity  dis¬ 
tribution  during  compression  and  expansion  strokes  is  shown  in  Fig. 

(  39  ) ,  which  indicates  lower  values  of  the  radial  velocity  as  com¬ 
pared  to  the  values  of  the  tangential  one. 

Figure  (  4o  ) shows  the  predicted  turbulence  intensity  contours  at  dif¬ 
ferent  port  area  ratio  of  0.0041,  0.0064,  0.0092,  0.0125  and  0.0164, 
at  20°  BTDC  and  at  engine  speed  of  1000  r.p.m.  The  increase  of  the 
ejected  air  velocity  from  the  main  cylinder  with  the  decrease  of  the 
port  area  ratio  has  a  great  effect  on  the  creation  of  a  strong  turbu¬ 
lence  intensity  field  in  t.,o  swirl  chamber  which  is  clearly  evident 
in  Fig. ( <J0  ) .  Regions  of  intensive  turbulence  are  noticed  near  the 

lip  of  the  connected  port  due  to  the  high  flow  shear  stresses  produced 
with  the  inlet  ejected  flow.  The  turbulence  level  has  a  maximui  value 
of  about  2.4  near  the  port  at  port  area  ratio  of  0.0041.  An  increase 
of  the  turbulence  intensity,  which  is  formed  in  the  central  zone  of 
the  swirl  chamber,  takes  place  when  the  port  area  ratio  decreases. 


Fig.  { 36  /Velocity  vector  plot  for  different 
port  area  ratios  at  2Cf  BTDC  keeping 
Vs/ Vc= 0.772  and  N=1000  r.  p.  m  constant . 
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Fig.(  38  )Predicted  and  measured  tangential  velocity  in  plane4at 
TDC  and  various  radii  for  different  port  area  ratios  keeping 
Vs/Vc  s0.772  and  Ns.1000 r.p.m.  constant. 


F/G  21/5 


AD-Allb  818  CAIRO  CENTER  FOR  COMBUSTION  (EGYPT) 

SPRAY  MODELLING  FOR  MULTIFUEL  ENGINES. (U) 
JUL  62  M  M  ELKOT0 


UNCLASSIFIED 


DA-ER0-79-G-0017 

NL 


-  9&  - 


0.0  0.25  0.5  0.75  1.0 


Dimesionless  radial  distance  r/Rs 
Fig.(  35  )Predicted  and  measured  tangential  velocity  in  plane4at 
TDC  and  various  radii  for  different  port  area  ratios  keeping 
Vs/Vc  =0.772  and  N  =  1000r.p.m.  constant. 
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Fig(  39  )  Predicted  and  measured  radial  velocity  in  plane  4  at  r/  Rs=0.503 
and  various  cran  k  angle  for  different  port  area  ratios,  kee  p  ing  Vs/Vc= 
0  772  and  N  =  1000  r.  p.  m-  constant. 


However  the  turbulence  intensities  elsewhere  in  the  swirl  chamber  have 
larger  values  at  low  port  area  ratio.  It  has  a  value  of  0.7  near  the 
swirl  chamber  wall  at  port  area  ratio  of  0.0125  and  increases  to  a 
value  of  1.78  at  port  area  ratio  of  0-0041. 

The  effect  of  the  port  area  ratio  on  the  turbulence  intensity  dis¬ 
tribution  during  compression  and  expansion  strokes  is  shown  in  Fig. 

(  41  >  •'  It  is  clear  from  the  Figure  that  the  turbulence  intensity 

h>-  a  maximum  value  reaches  about  3.0  with  port  area  ratio  of  0.0041 
at  about  60°  ATDC .  The  increase  of  the  turbulence  intensity  during 
compression  and  expansion  strokes  with  the  decrease  of  port  area  ratio 
may  be  attributed  to  the  produced  air  flow  jet  from  or  to  the  swirl 
chamber.  The  experimental  results  in  Fig.  (41  )  show  that  the  crank 

angle  at  which  the  maximum  turblenoe  intensity  occurs  (\j  =  60°  ATDC) 
is  independent  of  the  port  area  ratio.  The  predicted  results  show  a 
slight  deviation  from  the  experimental  ones  at  high  port  area  ratio. 

Figure  (  42  )  shows  the  effect  of  the  port  area  ratio  on  the  turbu¬ 

lence  intensity  distribution  along  the  radial  distance.  It  is  clear 
from  the  figure  that  for  all  ratios  of  port  area  the  maximum  turbulence 
intensity  exists  near  the  chamber  centre  due  to  the  double  effect  of 
both  shear  stresses  and  flow  inertia  generated  during  compression 
stroke  which  forms  a  strong  recirculation  zone.  Another  maximum  tur¬ 
bulence  intensity  appears  near  the  chamber  wall  for  a  1  i.  port  area  ra¬ 
tios  due  to  nigh  shear  stresses  resulting  from  the  e  pansion  of  the 
air  jet.  The  minimum  value  of  the  turbulence  intensity  is  0.35  occur- 
ir.g  at  a  port  area  ratio  of  0.0164  and  a  radial  distance  of  about 
12.6  r.mi. 

i.u.3.  Effect  of  Swirl  Volume  Ratio  on  the  Flow  Field 

The  effect  of  swirl  chamber  volume  ratio  as  defined  by  the 
ratio  of  swirl  chamber  volume  (Vs)  to  the  total  clearance  volume  (Vc) 
is  studied  for  five  swirl  chamber  volume  ratios  of  0.772,  0.54,  0.545, 
0.457  and  0.376,  and  at  engine  speed  of  1000  r.p.m.,  constant  com¬ 

pression  ratio  of  17.0,  and  port  area  ratio  of  0.0092,  to  investigate 
its  effect  on  both  flow  velocities  and  turbulence  intensity.  The  study 
indicates  that  by  increasing  the  swirl  volume  ratio  (Vs/Vc)  the  ejec¬ 
ted  flow  velocity  and  consequently  the  flow  field  velocities  and  tur¬ 
bulence  intensity  increase. 

Figure  (  43  )  shows  the  predicted  results  of  the  effect  of  the  swirl 
volume  ratio  (Vs/Vc)  on  the  flow  field  velicities.  It  is  clear  that  the 
velocities  decrease  with  the  decrease  of  the  swirl  volume  ratio,  and 
consequently  a  corresponding  decrease  of  the  two  eddies  formed  on  both 
sides  of  the  connected  port  is  obtained.  It  is  clear  also  that  the 
swirl  volume  ratio  has  a  higher  effect  on  the  flow  field  velocities 
than  the  port  area  ratio.  This  effect  can  be  noticed  by  comparing  the 
results  illustrated  in  Figs.  (  44  )  and  (37  ) .  This  comparison  shows 

that  a  maximum  tangential  velocity  of  102  m/s  is  reached  when  using  a 
port  area  ratio  of  0.0041,  while  a  maximum  tangential  velocity  of  60 
m/s  is  reached  when  .using  a  swirl  volume  ratio  of  0.772.  The  same 
results  is  obtained  when  comparing  the  results  illustrated  in  Figs. 

(  45  )  and  (  38)  for  the  tangential  velocity  distribution  along  the 

radial  distance  . 

Yhe  experimental  and  predicted  results  indicate  also  that  the  radial 
distance  ratio  at  which  the  max ’’muni  tangential  velocity  occurs  (about 
r/Rsa  0.75  1  I  «  ;  u  m  •  :  ■  ■  r<  .’a,  i  ■  swirl  chamber  volume  ratio. 
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Fig  (41  )  Predicted  and  measured  turbulence  intensity  in  plane  4 
at  r/Rs  =0-503  during  compression  and  expansion  strokes  f  or  different 
port  area  ratio  keeping  Vs/Vcr0.772  and  N  =  1000r.p.m.  constant. 
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Fig(42  )Predicted  and  measured  turbulence  intesity  in  plane4 
at  TDC  and  vaious  radii  for  different  port  area  ratios  keeping 
Vs/ Vc  =  0  772  and  N=  lOOOr.p.m  constant. 
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Fig.(  44  )Predicted  and  measured  tangential  veloci  ty  during  compression 
and  expansion  strokes  in  plan e  4  at  r/Rs  =0  503  for  different  swirl  volume 
ratios  keeping  Aj/A  =  00092  and  N=1000  r.  p.m.  constant. 
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Fig. (45  JPredicted  and  measured  tangential 
velocity  in  plane  A  and  TDC  for  various  radii 
and  different  swirl  volume  ratioskeeping  Aj/A 
0.0092  and  N=  1000  r.  p.m.  constan f  • 


The  of f I’i'L  of  the  swirl  chamber  volume  ratio  (Vs/Vr)  on  the  radial  ve¬ 
locity  distribution  during  compression  and  expansion  strokes  is  shown 
in  Fig.  (  4 (>  ),  from  which  it  can  be  seen  that  i  no  effect  of  swii  1 

chamber  volume  ratio  on  the  radial  v  '  >city  component  is  smaller  than 
that  of  the  port  area  ratio. 

Figure  (  )  shows  the  predicted  turbulence  intensity  for  various  relative 

swirl  chamber  volume  ratios.  It  is  clear  from  the  figure  that  the 
turbulence  intensity  increases  in  all  swirl  chamber  space  with  the  in¬ 
crease  of  the  swirl  chamber  volume  ratio,  taking  into  consideration  that 
the  compression  ratio  ’as  a  const.".!  va  1  m'  of  ]~.o  for  all  re  1  at  i  ve  swirl 
chamber  volume  ratios.  A  slight  increase  of  the  turbulence  intensity 
value  near  the  swirl  chamber  wall  is  noticed  with  the  increase  of  the 
swirl  chamber  volume.  Values  ol  about  0.7  and  0.91  are  obtained  at 
0.545  and  0.772  swirl  volume  ratio  respectively . 

The  two  eddies  generated  at  both  sides  of  the  part  begin  to  decrease 
and  spread  around  the  inlet  port  and  the  right-hand  side  eddy  has  a 
value  of  about  0.75  at  a  swirl  chamber  volume  ratio  of  0.457.  It  is 
clear  from  Fig. (  47  .)  also  that  the  maximum  turbulence  value  is  obtain¬ 
ed  near  the  centre  of  the  chamber  with  high  shear  stress.  Figure 
(  48 )  illustrates  the  effect  of  the  swirl  chamber  volume  ratio  on  the 

turbulence  intensity  distribution  along  the  radial  distance  of  plane 
4  at  BTDC .  It  is  clear  from  this  figure  that  the  maximum  turbulence 
intensity  of  1.35  occurs  near  the  centre  when  using  a  swirl  chamber 
volume  ratio  of  0.772.  A  minimum  turbulence  intensity  of  0.485  occurs 
at  about  0.65  radial  distance  ratio  when  using  0.357  swirl  volume  ratio. 

Figure  (  49  )  shows  the  effect  of  the  swirl  chamber  volume  ratio  on 
the  turbulence  intensity  distribution  during  compression  and  expansion 
strokes  in  plane  4  at  12.2  mm  radial  distance  and  1000  r.p.m.  engine 
speed.  It  is  clear  from  this  figure  that  the  turbulence  intensity 
generated  during  the  compression  stroke  increases  with  the  increase 
of  the  swirl  chamber  volume  ratio,  and  takes  on  a  maximum  value  at 
about  60°  ATDC . 


4.3.4.  Effect  of  Engine  Speed  on  the  Flow  Field 


The  effect  of  engine  speed  on  the  flow  field  in  the  swirl 
chamber  has  been  studied  by  varying  the  engine  speed  from  800  to  1500 
r.p.m  while  keeping  all  other  factors  constant.  Figure  (  50  )  shows 
the  effect  of  engine  speed  on  the  flow  I ield  inside  the  swirl  chamber. 
It  is  clear  from  the  figure  that  the  velocities  at  each  point  in¬ 
creases  with  the  increase  cf  the  engine  speed,  and  the  eddies  generated 
in  both  sides  of  the  port  are  strongly  increased  with  the  increase  of 
the  engine  speed. 


The  effect  of  the  engine  speed  on  the  tangential  velocity  distribution 
at  various  crank  angles  during  compression  and  expansion  strokes,  and 
at  various  radial  distances  is  shown  in  Figs. (  51  )  and  (  52  )  res¬ 
pectively.  It  is  clear  from  Fig.(  51  ’  that  the  tangential  velocity 
increases  with  the  increase  of  engine  speed  during  both  compression 
and  expansion  strokes.  The  decrease  of  the  tangential  velocity  occur- 
ing  during  the  reverse  of  flow  of  the  ejecting  air  to  the  main  cylin¬ 
der  increases  with  the  increase  of  engine  speed.  The  second  maximum 
tangential  velocity  occuring  at  about  22°  ATDC  is  independent  on  engine 
speed.  The  maximum  tangential  velocity  component  has  a  value  of  about 
92  m/s  at  1500  r.p.m. 
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Fig  (  46  )  Predicted  and  measured  radial  velocity  in  plane  4at  r/Rs-0.503 
and  various  crank  angles  for  di  fferen  t  swirl  volume  rati  os  keeping 
Aj  /  A  =  0  0092  and  N  =  1 000  r  p.  m  const  an  t. 
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Dimensionless  radial  distance  r/Rs 
Fig. ( 48  )Predicted  and  measured  turbulence  intensity  in  plan4ut 
TDC  and  various  radii  for  different  swirl  volume  ratios  keeping 
Aj/A  =  00092  and  N  =  1000  r.p.m  constant. 
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Fig(  5o)Velocity  vector  plot  for  different 

engine  speed  at  20°  B  T  DC  keeping 
Vs/Vc=0772  And  A) /A  =0  0092  constant  . 
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Fig .(  JPredicted  and  measured  tangential  velocity  during  compression 
and  expansion  strokes  in  plane  4  for  r/Rs=0  503  at  various  engine 


speeds  keeping  Vs/Vc  =  0772  and  Aj/A  =  00092  constant- 
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Fig. (52  )  Predicted  and  measured  tangential  velocity  in  plane  4 

at  TDC  for  different  engine  speeds  keeping  Vs/Vc  =  0.772  and 
A  j  /  A  =0.0092  constant 


Fiyure  (  '-'2  )  shows  also  that.  the  tangential  velocity  increases  with 
the  increase  of  the  engine  speed  along  the  radial  distance,  and  that 
the  position  of  the  maximum  value  is  independent  on  the  engine  speed. 

It  has  a  value  of  about  123  m/'s  for  1300  r.p.m.  engine  speed  at  <bout 
0.75  radial  distance  ratio  (r/Rs) .  Tt  is  noticed  that  at  this  crank 
angle  (TDC)  all  tangential  velocities  for  all  engine  speeds  are  de¬ 
composed  into  two  regions,  solid  vortex  region  near  the  chamber  centre, 
and  free  vortex  one  near  the  chamber  wall.  These  two  regions  are 
strenthened  with  the  increase  of  engine  speed. 

The  effect  of  the  engine  speed  variation  on  the  radial  velocity  dur¬ 
ing  compression  and  expansion  strikes  is  shown  in  Fig. (  53  ).  It  is 
clear  from  the  figure  that  the  radial  velocity  increases 
witli  the  the  increase  of  engine  speed  reaching  a  maximum  value  of 
about  11.3  m/s  at  1500  r.p.m.  *  about  40°  BTDC .  and  then  decreases 
again  to  zero  at  the  end  of  the  expansion  stroke  (BDC) . 

Figure  (  54  )  shows  the  effect  of  engine  speed  on  the  turbulence  in¬ 
tensity  field  inside  the  swirl  chamber.  It  is  clear  that  the  engine 
spedd  has^slightly  lower  effect  on  the  turbulence  intensity  than  the 
effect  of  relative  port  area  ratio,  and  relative  swirl  volume  ratio. 

The  two  eddies  which  are  formed  on  both  sides  of  the  port  change  from 
a  value  of  about  1.2  at  engine  speed  of  800  r.p.m.  to  about  1.1  at 
engine  speed  of  1000  r.p.m.  This  small  effect  of  the  engine  speed 
on  the  turbulence  intensity  is  attributed  to  the  effect  of  engine 
speed  on  the  kinetic  energy  of  turbulence  K  and  the  mean  piston  speed 
Cpn1  determining  the  turbulence  intensity.  The  effect  of  engine  speed 
on  the  turbulence  intensity  distribution  with  crank  angle  degree  and 
radial  distance  is  shewn  in  Figs  (  55  )  and  (  56  )  respectively. 

Figure  (  55  )  shows  the  effect  of  engine  speed  on  the  turbulence  in¬ 
tensity  distribution  during  compression  and  expansion  strokes  for  dif¬ 
ferent  engine  speeds  at  TDC  in  plane  4,  from  which  it  is  clear  that 
the  engine  speed  has  a  very  low  effect  on  the  turbulence  intensity. 

Figure  (  56  )  shows  the  effect  of  engine  speed  on  the  turbulence  in¬ 
tensity  distribution  along  radial  distance  of  plane  4,  from  which 
the  slight  increase  of  turbulence  intensity  with  the  increase  of  engine 
speed  is  clear.  It  is  clear  also  from  the  near  wall  results  that  the 
turbulence  intensity  values  are  independent  on  the  engine  speed, 
while  an  appreciable  increase  of  the  turbulence  intensity  appears  near 
swirl  chamber  centre.  A  value  of  about  1.51  at  engine  speed  of  1500 
r.p.m.  is  reached.  The  minimum  value  of  the  turbulence  intensity  is 
obtained  at  about  r/Rg  =  0.51  and  lias  a  value  of  about  0.71  at  engine 
speed  of  800  r.p.m.,  its  location  is  independent  on  the  engine  speed. 


Radial  velocity  v  mis 


r 


Crank  angle  degree  tji 


Fig  (  53)Predicted  and  measured  radial  velocity  at  plane  4  and r/Rs  = 

0503  during  comp ression  and  expansion  strokes  at  various  engine 
speeds  keeping  Vs/Vc  =  0772  and  Aj  /  A  =  000  9  2  const  an  t . 
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Fig(  5s) Predicted  and  measured  turbulence  intensity  at  r/Rs  =  0503  m 
plane4for  various  crank  angles  and  ditferent  engine  speeds  keeping 
Vs/Vc  =  0.772  and  Aj/A  =  0  0092  constant 


Dimensionless  radial  distance  r/Rs 
Fig(  56  ^Predicted  and  measured  turbulence  intensity  at  T  DC  in 
plane4for  different  radii  and  different  engine  speeds  k  ee  p  i  n  g 
Vs/Vc  =  0772  and  Aj/A=  0.0092  constant 
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constants  in  Tanasaw-Tes  i ma  Eqn.  1 4  ) 

2 

area,  m 

mass  transfer  coefficient 

discharge  coefficient. 

constants  of  the  turbulence  model 

absolute  velocity  of  fuel  droplet,  m/s. 

nozzle  diameter,  m 

ligament  diameter,  m. 

droplet  diameter,  m 

diffusion  coefficient. 

diameter  of  fuel  ligament,  m 

relevant  mean  diameter  depending  on  values  of  p  and  q 
such  as  D^0  -linear  mean  diameter,  D^-  Sauter  mean 
diameter,  jam. 

energy  lost  in  forming  new  ligament  J. 
activation  energy 

degree  of  freedom  for  chi-square  distributions 

percentage  by  volume  of  diesel  fuel  in  the  mixture. 

2 

heat  transfer  coefficient,  W/m  hr  C 
total  enthalpy,  W. 

constants  of  rate  Eqn .  (1,2)  and  Eqn .  (32). 

kinetic  energy  of  turbulence 

thermal  conductivity,  W/m  C. 

spray  penetration,  m. 

length  of  fuel  nozzle,  m. 

latent  heat  of  fuel,  J/Kg  fuel. 

theoretical  quantity  of  air  needed  for  complete 
combustion  of  fuel,  Kg/Kg  fuel, 
mass,  Kg. 

index  in  rate  Eqn.  (1) 

Ohnesorse  Number,  dimensionless  value  p  ?  / (p  d  O'  ) 
total  mass  of  fuel  injected,  Kg. 
total  mass  of  unprepared  fuel, Kg. 
number  of  fuel  droplets. 
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engine  speed,  RPM 

Nusselt  number,  hD/K,  dimensionless 
pressure,  bar 

pressure  drop  across  fuel  nozzle,  bar, 

normalized  chamber  density  (o  /o  )  p  ,  bar. 

1  a  a  c  o 

preparation  rate, 
partial  pressure  of  oxygen,  bar. 
prandtle  number,  p  c  /K,  dimensionless 
normalized  fuel  quantity  per  cycle,  dimensionless, 
radius ,m . 

Universal  qas  constant. 

instantaneous  heat  release  rate  or  reaction  rate. 
Rosin-Rammler . 

Reynolds  number,  CDp/p  ,  dimensionless 
gas  constant  of  fuel  vapour,  J/Kg  °C . 
source  vector  of  Eqn.  (7). 

volumetric  source  rate  of  1  arising  from  gas 
volumetric  source  rate  of  i>  arising  from  fuel 
Schmidt  number,  droplet,  dimensionless, 

Sherwood  number,  dimensionless 
time,  S. 

absolute  temperature , °K 
Tanasawa-Tesima 

tangential  velocity  of  air,  m/s. 

cumulative  volume  of  droplets  less  than  a  given  size 
D,m^ . 

radial  velocity  of  air,  m/s. 

incremental  value  of  drops  within  the  size  range 
3 

D , D+dD ,  m 

relative  velocity  of  droplet,  m/s. 

Weber  number 

index  in  rate  Eqn.  (1) . 

droplet  velocity  in  x  and  y  direction  . 
accelaration  of  droplet  in  x  and  y  direction, 
excess  air  factor,  m  /(L  m  ),  dimensionless. 

cl  C 

constants  of  T-T  Eqn. (41) . 

optimum  length/diameter  of  ligament 

K-i  -x 

Gamma  tunction  !'(K)  =  j  x  e  dx . 
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eddy  duffusivity  for  <t> . 

inclination  of  droplet  relative  velocity  to  x-axis. 
dissipation  rate  of  kinetic  energy  of  turbulence, 
angular  position 

crankangle,  degree,  arbitrary  dependent  variable, 
injection  angle,  degree 

2 

dynamic  viscosity,  kg  s/m  . 

2 

Kinematic  viscosity,  m  /s. 
density,  kg/m'3 
surface  tension,  N/m. 

constants  of  turbulence  model  for  enthalpy , kinetic 
energy  dissipation  rate  and  concentration. 


Subscripts 

a  air 

c  chamber 

f  surface 

i  arbitrary  value 

in  inlet  to  control  volume 

l  liquid  fuel 

o  initial 

out  outlet,  to  control  volume, 

p  piston  cavity 

r  radius 

s  saturation 

v  vapour 

x  x-direction 

y  y-direction 

free  streem. 
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CHAPTER  5 


FUEL  SPRAY  MODELS  INC 


5.1  INTRODUCTION 

The  fuel-air  mixture  creation  prior  to  js  we))  is  after 
auto  i  gn  i  t:  i  on  is  important  in  diost  1  engines.  It  fontrols  the  Late 
of  heat  release  and  consequently  rate  of  pressure  rise,  noise,  and 
pollutants  formation.  The  mixing  process  depends  upon  the  kinetic 
energy  of  both  fuel  and  air.  In  open  combustion  chambers,  since  the 
air  motion  is  small  the  fuel  is  injected  under  high  pressure 
through  a  nozzle  with  one  or  several  holes  to  ensure  adequate 
mixing  of  fuel  and  air.  Increasing  the  air  motion  during  the  suc¬ 
tion  stroke  by  a  mask  valve,  tangential  port  or  swirl  inlet  mani¬ 
fold  assists  in  decreasing  the  fuel  pressure.  The  air  motion  in 
quiescent  combustion  chambers  can  be  created  also  during  the  com¬ 
pression  stroke  by  using  bowl  piston.  The  energy  of  mixture  for¬ 
mation  in  swirl  chambers  is  mainly  generated  by  the  kinetic 
energy  of  the  air  swirl  motion  during  the  compression  and  expan¬ 
sion  strokes.  The  energy  of  the  fuel  spray  has  rather  smaller 
effects  on  mixture  formation  in  these  chambers. 

The  mixing  process  now  includes  transport  of  fuel  from  the 
jet  by  swirling  air  flow,  deflection  of  the  jet,  impingement  of 
the  jet  on  the  combustion  chamber  wall.  Besides,  the  impinged 
jet  is  spread  on  the  walls  of  the  combustion  chamber  forming  a 
film  of  very  small  thickness.  This  film  takes  its  required  heat 
of  evaporation  from  the  walls  and  surrounding  gases  by  convection. 
The  injection  and  mixing  of  fuel  with  air  in  the  combustion  cham¬ 
ber  needs  to  be  achieved  very  rapidly  if  the  engine  is  to  be  ef¬ 
ficient.  Many  investigators  have  studied  the  problem  of  mixture 
formation  and  combustion  in  diesel  engines,  but  most  of  them 
treated  the  problem  either  in  general,  or  by  formulating  empiri¬ 
cal  relations  for  the  individual  processes. 
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The  use  o£  modi’rn  computers  makes  i  h.>s  i ! ,  I  (<  the  calcula¬ 
tions  of  more  realistic  and  complex  engine  cycles  ,  and  these 
in  turn  require  mathematical  formulation  of  the  combustion  pro¬ 
cess.  A  major  purpose  of  early  studies  was  to  calculate  the  cy¬ 
linder  pressure  diagram  using  a  predetermined  heat  release  pat¬ 
tern  obtained  by  analysing  experimental  cylinder  pressure  diag¬ 
ram  of  similar  engines.  The  cycle  calculations  involve  adding 
energy  to  the  cylinder  contents  using  small  steps  and  consider¬ 
ing  uniform  cylinder  conditions  of  composition,  pressure  and 
temperature.  This  method  of  combustion  modelling  is  called 
single-zone  model.  It  suffers  from  severe  limitations  if  the 
effect  of  various  factors  on  combustion  process  are  to  be  inves¬ 
tigated.  Modifications  were  proposed  to  include  the  fuel-air 

1  2 

mixture  prep* ration  by  Austen  and  Lyn  '  in  which  the  fuel  injec¬ 
tion  process  is  divided  into  a  number  of  steps.  The  fuel  inject¬ 
ed  during  one  step  is  calculated  according  to  a  simple  mathemati¬ 
cal  formula.  Summing  all  such  steps  would  then  give  the  prepara¬ 
tion  rate  diagram.  However,  the  used  formulae  neglect  oxygen 
availability.  Whitehouse  et  al3'4  suggested  semi-empirical  for¬ 
mulae  for  the  preparation  and  combustion  rate  of  fuel.  The  rate 
of  preparation  was  assumed  depenendent  upon  the  total  surface 
•  iff  »  i<f  rue  dr  s  former  j  t  he  fuel  spray.  Considering  homo- 

\izo  i  i.st.r  ibution  of  size  equal  to  the  mean  drop¬ 
let  diameter  the  preparation  rate  is  represented  by 

P  =  k1m1(1‘x,mu  ro2  (l) 

and  the  reaction  rate  is  calculated  by  the  following  equation 

k2  P 

=  -N  ^__-2  Exp(-Eact/kT  )  /( P  -  H.j  d<t>  (2) 

Inaccurate  results  are  obtained  when  applying  this  method 
to  quiescent  combustion  chambers  where  jet  entrainment  is  of 
major  performance.  In  this  case  the  single  zone  model  cannot  es¬ 
timate  the  effect  of  number  of  orifices  in  the  nozzle.  It  cannot 
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indicate  the  lack  of  air  and  production  of  pul  lut  ing  subtances 
which  are  dependent  on  local  temperature  and  chemical  concentration. 

A  cwo-zone  model  has  been  developed  to  improve  the  previous 

model.  In  this  model  the  combustion  chamber  contents  are  divided 

into  two  parts,  the  burning  and  the  nonburning  zone.  Formulae  of 

preparation  and  combustion  rate  of  the  single  zone  model  are  used  in 

the  burning  zone  and  the  combustion  products  are  assumed  ideal  gases. 

The  nonburning  zone  consists  of  pure  air  which  is  treated  as  a 

perfect  gas.  Such  two-zone  model  has  been  developed  for  both 

quiescent"*  and  swirl^  combustion  chambers.  The  burning  zone  has  been 

assumed  as  a  cone  jet  of  known  cone  angle  whose  development  is  determined 

from  the  fuel  spray  penetration  data.  The  development  of  the  jet  of 

fuel  plus  air  prior  to  impingement  on  the  walls  has  been  investigated  by 
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numerous  research  workers  and  simple  formulae  for  the  penetrated 

spray  are  suggested.  Once  the  jet  has  struck  the  combustion  chamber 
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wall  it  is  treated  as  wall  jet  .  Such  model  does  not  deal  with 

the  motion  of  either  the  air  or  fuel  jet. 


The  major  problem  in  the  aforementioned  models  is  to  estimate 

the  rate  of  mixing  of  the  gas  outside  the  burning  zone  with  the 

gases  inside  this  zone.  Empirical  or  semi-empirical  estimates  are 

used.  High  speed  photography  of  the  spray  shows  that  when  the  gas 

flows  perpendicular  to  the  fuel  spray,  the  spray  deflects  in  the 

direction  of  the  gas  flow  and  loses  its  symmetry.  Assuming  same  drop 

size  distribution  and  mean  droplet  diameter  for  every  increment  of 

ruel  in  the  injection  process,  some  investigators attempted  to 

;se  spalding's  single-droplet  results  to  predict  burning  rates  in 

an  engine.  The  concept  of  a  finite  spray  burning  rate  based  on 
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single  droplet  was  also  proposed 


assuming  no  interaction  between 


droplets.  However,  Shipinski^  ^  evaluated  these  models  and 

concluded  that  the  overall  performance  of  the  engine  is  highly 
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dependent  on  the  burning  rate.  Recently  few  investigators  have 

tried  to  formulate  more  fundamental  methods  for  predicting  mixing 
and  turbulent  flow.  Some  of  these  models  divided  the  fuel  spray  into 
a  lot  of  small  packages  and  no  mixing  among  them  was  assumed 
i  thers  tried  to  solve  the  equations  governing  the  droplet  motion  to 
define  the  spray  development.  However,  data  for  the  droplet  size 
distribution  and  the  effect  of  various  factors  upon  it  are  required 
f  it*  the  development  of  these  models. 
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Itis  the  aim  of  the  present  work  to  review  the  fuel  spray 
models  and  to  develop  a  new  spray  model  for  the  calculation  of 
fuel  spray  trajectory  in  a  swirl  chamber  diesel  engine  based 
on  the  combination  of  the  discrete  droplet  model  with  a  multi¬ 
dimensional  gas  flow  prediction  method  .  The  model  includes 

turbulence  interaction  between  fuel  and  air  ,  variation  of  fuel 
droplet  diarn.  i er  and  consiquently  drag  force  and  hen*  transfer 
to  the  fuel  droplet  . 


5.2  REVIEW  OF  FUEL  SPRAY  MODELLING 


Overpenetration  of  fuel  spray  in  quiescent  chamber  would 
lead  to  impingment  upon  a  cold  surface  giving  consequent  fuel 
wastage.  Early  studies  of  fuel  penetration  were  mainly  based 
on  the  correlations  of  the  experimental  data.  Along  this  line 
many  investigators  have  carried  out  experimental  measurements  of 
the  penetration  of  a  single  spray  injection  into  air  using 
different  techniques  3°-31  Hay11  surveyed  these  penetration  for¬ 
mulae  and  concluded  that  Dent's  expression  of  the  following  form: 

■%  =  0.865  |  ftP/b  }*  t  d  |  *  (-^)1/4  (3) 

a  a 

can  be  recommended  for  small  densities.  However,  spray  penet¬ 
ration  in  diesel  combustion  chambers  varies  with  the  operating 
parameters  and  geometrical  design  of  injectors.  Therefore,  these 
correlations  are  limited  in  their  applications  and  a  theoretical 
calculation  approach  of  the  fuel  spray  is  needed. 
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Riehm  ,  Triebnigg  and  Sitkex  were  the  first  who  attempted  to 
calculate  the  spray  penetration.  They  treated  the  spray  motion 
as  constant  radius  spheres  leaving  the  injection  nozzle  under 
aerodynamic  resistance  forces.  They  applied  the  law  of  conser¬ 
vation  of  energy.  Mehlig38  analysed  the  experimental  results  of 


Ti  LeLx.igg  and  assumed  that  t  ho  lateral  m< •  •. i  m  o f  the  spray  re¬ 
sults  from  turbulence  which  depends  mi  the  injection  system 
uni'/.  1'his  leads  to  that  the  radial  velocity  of  the  spray  does 
not  depend  on  the  back  pressure.  lie  deduced  an  expression  for 
sprat  penetration  from  the  geometrical  relations. 

Hakki  OZ  calculated  the  spray  penetration  using  the 
differential  equation  of  free  jet  and  considering  the  influence 
Of  turbulent  friction  force  alone.  Making  use  of  the  available 
experimental  data  he  proposed  the  following  formula  for  the  spray 
penetration  : 


=  4.8! 


-0.04  p . 
e  d 


Ogasawara  and  Sarni^0  have  attempted  tv  investigate 
theoretically  the  behaviour  of  a  single  droplet  representing  the 
mean  droplet  diameter  and  moving  at  the  spray  tip.  The  equation 
of  motion  of  the  droplet  and  the  equation  of  heat  balance  have 
been  derived  during  the  delay  period  assuming  that  the  droplet 
diameter  was  constant  and  that  the  air  velocity  had  no  effect. 

A  closed  form  integration  of  the  equation  was  obtained  giving 
the  trajectory  of  the  droplet  during  the  delay  period  in  the  fol¬ 
lowing  form  : 


:  c  ,d  /',i 
o  i  o 
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1  n  |  1  +  y=-c'.-  ;  1-exp  (-37. 3  £  N,  )  :  l  5) 


"  "  \  =  **>*  d  /iJa  '  Rcf  =  Co  d/  a 

>'y  =  constant  =  0.85 

Considerable  discrepancy  appeared  between  the  theoretical  and  ex¬ 
perimental  results  for  kerosene  flame.  This  resulted  from  neg¬ 
lecting  relative  velocity.  Air  motion  contributes  to  the  magni¬ 
tude  of  Reynolds  number  and  thus  affects  the  droplet  behaviour 
and  the  heat  transfer  to  the  droplet.  Therefore, it  is  necessary 
to  take  the  air  motion  to  formulate  a  mathematical  model  based 
on  the  droplet  behaviour.  royasi  and  Knriota--  '  developed  models 
for  combustion  and  formation  of  nitric  oxide  and  soot  in  direct 


injection  diesel  engines.  In  this  model,  the  spray  injected 
into  a  combustion  chamber  from  an  injection  nozzle  is  divided 
into  a  lot  of  small  packages.  Each  package  experiences  its  own 
history  of  temperature,  pressure,  and  equivalence  ratio,  which 
govern  the  chemical  kinetics.  The  processes  occuring  in  the 
package  are  shown  in  Fig. (1) .  The  package,  immediately  after 
the  fuel  injection,  involves  a  lot  of  fine  droplets  and  a  small 
volume  of  air.  During  the  package  motion,  the  air  entrains  in¬ 
to  it  and  fuel  droplets  evaporate.  Therefore,  the  small  pac¬ 
kages  consist  of  liquid  fuel,  vaporized  fuel  and  air.  Shortly 
after  the  start  of  injection  ignition  occurs  in  the  gaseous  mix¬ 
ture  and  a  sudden  expansion  of  the  package  takes  place.  There¬ 
after,  the  fuel  droplets  evaporate  and  fresh  air  entrains  into 
the  package.  Vaporized  fuel  mixes  with  fresh  air  and  combus¬ 
tion  products  and  spray  continues  to  burn.  All  fuel-air  packages 
are  assumed  to  burn  at  the  stoichiometr ic  condition.  The  effect 
of  air  motion  and  the  interaction  of  neighbouring  droplets  have 
not  been  taken  into  account  in  this  analysis. 

The  only  studies  of  spray  injection  into  gaseous  flow 

fields  which  are  based  on  continuum  mechanics  are  those  of  Wakuri 
40  4 1 

et  al  and  Cotton  et  al  .  These  studies  enable  us  to  calculate 
the  spray  penetration,  but  neither  the  evaporation  rate  nor  the 
mixing  field  details,  such  as  velocity,  temperature  and  concen¬ 
tration  distribution,  can  be  calculated  with  these  methods.  Adler 
4  2 

and  Dyn  developed  a  theoretical  model  based  on  the  continuum 
mechanics  basis  using  the  classical  differential  equations  of 
continuity,  momentum  energy  and  diffusion.  A  differential  volume 
element  having  a  mass  equivalent  to  the  total  mass  of  air,  liquid 
droplets  and  evaporated  liquid  contained  in  the  element  was  used. 
The  element  is,  therefore,  characterized  by  its  concentration 
and  evaporation  functions. They  treated  the  transient  phenomenon 
on  a  steady  state  basis  by  taking  the  elemental  section  of  the 
spray  as  a  quasi  -  solid  body  under  equilibrium  of  the  inertia 
and  drag  forces  in  the  flow  field.  The  cross  section  of  the 
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Fig(  1  )  Schematic  diagram  of  the  spray  combustion  model 
of  Hiroyasu 


spray  was  replaced  by  a  rectangular  cross  section  with  constant 
height/width  ratio  to  formulate  a  two  dimensional  model.  An 
integral  method  was  used  to  replace  the  partial  differential 
equation  by  a  set  of  ordinary  differential  equations  which  was 
then  solved  to  determine  the  spray  centre  line.  Thereafter, 
the  velocity,  concentration,  temperature  and  density  structures 
are  built  on  this  framework. Such  a  model  cannot  determine  the 
exact  values  of  the  local  concentration  especially  when  there 
is  interference  from  neighbouring  sprays.  It  can  however  be 
generalized  for  all  types  of  combustion  chambers. 

Spray  modelling  based  on  the  combination  of  the  discrete 

droplet  model  with  a  multidimensional  gas  flow  prediction  method 

was  developed  for  the  prediction  of  turbulent  combustion  of 

4  3 

liquid  fuels  in  steady  flow  combustors  and  in  direct  injection 
4  4 

engines  .  However,  these  models  are  still  awaiting  sufficient¬ 
ly  detailed  information  about  the  droplet  size  distribution  and 
droplet  velocity  distribution. 

During  the  past  seven  years  comprehensive  works  were  done 
at  Cairo  University  about  spray  modelling.  Firstly  the  spray 
calculation  model  was  based  on  the  solution  of  the  equation  of 
motion  of  droplets  and  air.  The  interaction  between  the  droplet 
and  air  was  considered  by  averaging  the  over  flow  processes  on 
a  scale  smaller  than  the  droplet  size  and  employing  typical  cor¬ 
relations  for  the  droplet  drag,  heat  and  mass  transfer.  The  re¬ 
sults  of  this  work,  which  have  heen  reported  in  several  publica- 
15  28,29  45,46  , 

tions  were  sufficiently  encouraging  to  extend  the 

model  to  include  the  prediction  of  air  motion  and  turbulence  in¬ 
teraction  between  fuel  and  air. 
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Ga s  M o t ion  and  Mixing 

The  flow  inside  diesel  eii.j-.nt  swir i  chambers  is  a  tr.im-- 
iont>  turbulent,  two  dimensional  flow.  The  flow  is,  thus,  govern¬ 
ed  by  the  differential  conservation  equations  of  energy,  mass 
and  momentum.  The  e, nations  governing  the  flow  in  a  cylindrical 
swirl  chamber  of  a  diesel  engine  are  obtained  by  decomposing  the 
velocity,  pressure,  density  and  enthalpy  into  a  mean  and  a  fluc¬ 
tuating  value  respectively.  The  turbulent  diffusion  fluxes  which 
appear  in  the  governing  equations  are  modelled  using  Boussinesq 
approach  in  terms  of  the  mean  flow  gradients  and  eddy  diffusivi- 
ties.  These  eddy  diffusiv.it  ies  are  determined  by  solving  addi¬ 
tional  differential  conservation  equations  for  the  time  averaged 
kinetic  energy  of  turbulence  K  and  its  dissipation  rate  e  taking 
into  consideration  compressibility'  effects. 


Within  the  above  framework,  the  governing  set  of  equations 
may  be  represented  in  terms  of  a  single  equation  for  an  arbitrary 
dependent  variable  1  ,  as  mentioned  in  chapter  2  ,  as  follows  : 


ft  (^>+  r  fr  (rp  v‘*,+  r  fc  r  J 


-  ~  rf  = 

r2  3  0  1  -J 


S 


1 


where  1  stands  for  radial  velocity  V, 
total  enthalpy  h,  kinetic  energy  k  ar>d 
kinetic  energy  of  turbulence  . 


(6) 


tangential  velocity 
dissipation  rate  of 


u, 

the 


The  equation  represents  also  the  continuity  equation  by 
replacing  $  by  1.  The  definition  of  and  S ^  are  given  in 

Table  (1) . 


Droplet  Trajectory 


The  solution  of  the  fuel  droplet  trajectories,  size  and 
temperature  histories  is  necessary  for  the  evaluation  of  the 


source  terms  in  the  gas-phase  flow  equations.  The  droplet  is 
subjected  during  its  motion  to  the  effect  of  the  air  velocity 
and  therefore  the  equation  of  motion  of  a  droplet  is  given  by  : 


d 

dt 


(m,C? ) 


+  R 


s 


(7) 


The  source  vector  Rg  contains  terms  involving  the  pressure 
gradient,  Basset,  Magnus  and  gravity  forces.  The  pressure  gra¬ 
dient  and  Basset  terms  are  neglected  as  they  are  very  small 
quantities.  The  Safi  lift  and  Magnus  forces  are  also  neg¬ 
lected  because  the  droplets  are  not  in  a  high  shear  region  of 


the  gas 

flow 

• 

Substituting  for  the 

mass  and 

area  of 

droplet 

the  equation 

of 

motion  is  described 

by  the  following 

two  equa 

tions  : 

dr . 

r 

%• 

= 

-3  C;.r-I  CD 

k  a 

—  cos 

P 

c 
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> 
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r 

= 

—  sin 

0  n 

K 
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where  ; 

_  1 

6 

= 

tan  W£y  /W£x 

W£x 

= 

C£x  '  ux  =  *  -  ux 

W«y  '  Cly  '  uy  *  *  '  uy 


and 

ux  =  u  cosG  -  V  sin  9 
Uy  =  u  sin  0+  V  cos  0 

where  9  =  tan-1  — 2— 

r^-y 


The  drag  coefficient  can  be  represented  reasonably  well  by 
Ingebo's  equation47  multiplied  by  a  correction  coefficient 
3^  =  0.85  as  it  was  derived  for  a  non-evaporating  droplet*. 


CD  =  0 . 85x27/Re0' 84  (9) 

The  rate  of  decrease  of  droplet  mass  during  the  evapora¬ 
tion  period  is  given  by  : 
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Mot  ion  and  Mixing 

The  flow  inside  diesel  en,*ne  swirl  chambers  is  .1  rrun::- 
iontr  turbulent,  two  dimensional  flow.  The  flow  is,  thus,  tjovern- 
od  by  the  differential  conservation  equations  of  energy,  .mass 
and  momentum.  The  equations  governing  the  flow  in  a  cylindrical 
swirl  chamber  of  a  diesel  engine  are  obtained  by  decomposin'?  t.  :.c- 
velocity,  pressure,  density  and  enthalpy  into  a  mean  and  a  £1 'ac¬ 
tuating  value  respectively.  The  turbulent  diffusion  fluxc-s  w.u.ich 
appear  in  the  governing  equations  are  modelled  using  Roussinesq 
approach,  in  terms  of  the  mean  flow  gradients  and  eddy  diffuslvi- 
ties.  These  eddy  dif fusivities  are  determined  by  solving  addi¬ 
tional  differential  conservation  equations  for  the  time  averaged 
kinetic  energy  of  turbulence  K  and  its  dissipation  rate  *  taking 
into  consideration  compressibil ity  effects. 


Within  the  above  framework,  the  governing  set  of  equations 
may  be  represented  in  terms  of  a  single  equation  for  an  arbitrary 
dependent  variable  1  ,  as  mentioned  in  chapter  2  ,  as  £0} lows  : 
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where  i  stands  for  radial  velocity  V,  tangential  velocity  u, 
total  enthalpy  h,  kinetic  energy  and  dissipation  rate  of  the 
kinetic  energ.  of  turbulence  . 


The  equation  represents  also  the  continuity  equation  by 
replacing  t  by  1 .  The  definition  of  and  S ^  are  given  in 

Table  (1). 


i  i-  Droplet  Trajectory 


The  solution  of  the  fuel  droplet  trajectories,  size  and 
temperature  histories  is  necessary  for  the  evaluation  of  the 


(10) 


,  typg-pj  = 


-4  ■ 


dt 


where  •  g  is  the  mass  fraction  of  vapour  at  the  droplet  sur¬ 
face  and  ^  the  mass  fraction  of  vapour  in  the  free  stream.  The 
mass  transfer  coefficient  is  calculated  from  tne  Ranz  -  Marshall 
equation  :: 

.  Hh 

/  2  r,  ill) 


wnere  ; 


D  =  D  (T  /T  )  ( P  /p  ) 

c  o  s'  o  o'  a 


The  Sherwood  number  is  augemented  due  to  convection  as 
4  2 

Biro,  et  al  suggested  by  the  following  expression  : 


Sh  =  2+0.6  Re  Sc 


V3 


U2) 


Using  the  abo^e  equations  and  taking  into  consideration  that 
the  partial  pressure  of  vapour  at  the  surf _ce  is  equal  to  the 
saturation  pressure  corresponding  to  the  temperature  of  the  drop¬ 
let,  the  rate  of  change  of  droplet  radius  can  be  written  as  r 

dr.  CK  _  P 

».  -  Sh  n  s 


wnere 


dt 
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R,  T 
V  s 


Ps  =  2.39x10  Exp  (-6065/Ts) 


s,T. 


=  462  p0,043  +  0.06  (T  -673 ) ;  350  <  T  <  600  °C 


The  heat  balance  equation  for  the  droplet  is 

dT 


dr. 


dt1  =  471  r?  h(W+4lT  r?peL  dt~ 


(14) 


Thus,  the  rate  of  temperature  increase  can  he  written  as  : 


dt  r.c  i  (1g-1?.  2  r 


L  Sh  Dc 
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ip?  ! 

The  last  term  of  the  equation  is  equal  to  zero  during  the  sen¬ 
sible  heating  period.  The  Nusselt  number  is  obtained  from  the 


following  correlations'0'  multiplied  by  a  correction  factor 


p  =  z/e  -1  ”l’’ch  considers  tne  effect  of  mass  transfer 
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Nu  =  0.54  Re^ 


Nu  =  2  +  0.6  Re 


loo  <•  Re 
100  ‘  Re 


100000 


(16) 


where , 


z  =  '(cpv  ^dK*!u 

The  duration  of  the  sensible  heating  period  can  be  determined 
from  the  following  relation  ; 

T 

,  s  d  T, 

(17) 


fcs  =  3rfo^V 


T  h(Ta-Tf) 
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substituting  for  the  rate  of  change  of  droplet  radius  from  egua 
tion  (13) ,  the  equation  of  motion  (8)  can  be  written  in  the  fol 
lowing  form. 

a -  For  the  sensible  heating  period  : 

1 

'ofii  'V»*i  *2  B 
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Re 
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Re 


8.6 

0.84  ’  a 


b  -  For  the  evaporation  period  : 
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where 


x ; 

ii 

Sh 

D 

8.6  p 

2r 

l 

Rv  Ts 

Re0'84  ^ 

Y  =  3y 

Sh 

Dc  Ps 

3.6  pa 

2rZ 

P.  R  T  = 
CVS 

Be0’84  “t 

A  B 

v. 


A  = 


B  = 


B2  +  CS 


rP-y 

(x-u  —  —  + 


K  ■  V  £7' 
(y-u  £  -  v(r  -y)/K' ) 


(19) 


C 


K'  =  x2  +  (r  p-y ) 2 

Equations  (15,18  &  19)  are  used  to  determine  the  droplet 
positions  and  temperatures  along  the  trajectory. 

iii-  Source  Term 

When  the  fuel  is  injected  into  the  combustion  chamber,  a 
liquid  jet  is  formed  and  breaks  up  into  fine  droplets  after  a 
certain  period.  The  break  up  length  depends  on  various  factors 
which  will  be  discussed  below. 

In  this  model,  the  spray  which  is  injected  into  a  combus¬ 
tion  chamber  from  an  injection  nozzle  is  divided  into  five  in- 
troudction  locations.  The  injection  rate  diagram  shown  in  Fig. 

(2)  was  correlated  for  various  factors  and  an  empirical  relation 
is  given  in  the  following  fornr^ 

*  *  *2  *3 
Q  =0.9  *  +  1.3  *  -  1.2  *  (20) 

where , 

* 

Q  =  Q^/Q,  normalized  fuel  quantity=cumulative  fuel 
quantity/quantity  of  fuel  per  cycle. 

* 

$  =  $/$  normalized  injection  crank  anglesin jection  angle 

calculated  from  start  of  injection/injection  period. 


The  injection  period  is  discretised  into  a  number  of  equal  intervals 
equal  to  l°crank  angle.  The  quantity  of  fuel  injected  during  this  period 
is  calculated  from  a  detailed  atomization  model.  Dividing  the 
total  number  of  droplets  equally  on  the  five  introduction  loca¬ 
tions  and  using  the  droplet  size  distribution  it  is  possible  to 
identify  the  number  of  droplets  for  each  size.  In  the  absence 
of  a  detailed  model  for  the  droplet  velocity  it  is  possible  to 
assume  the  same  distribution  for  the  velocity  taking  into  con¬ 
sideration  the  variation  of  the  exit  velocity  due  to  the  varia¬ 
tion  of  the  fuel  and  combustion  chamber  pressures  as  shown  in 
Fig.  (3).  Thus,  every  droplet  is  assigned  different  initial  con¬ 
ditions  of  location,  size  and  velocity.  The  droplet  size  range 
is  covered  with  several  sizes  of  increment  10  jum  covering  the 


Rate  of  fuel  injection 


Fig(  2  )  The  fuel  delivery  law  and  injection 
law  at  various  engine  speeds 


Injected  fuel  quantity 


144 


range  till  the  maximum  droplet  diameter  according  to  the  prede¬ 
fined  conditions. 


The  flow  field  of  combustion  chamber  is  subdivided  into 
grids  as  shown  in  Fig.  (4)  and  each  grid  is  treated  as  a  control 
volume.  The  net  efflux  rate  of  droplet  mass  from  a  computation¬ 
al  grid  which  is  the  source  term  in  the  continuity  equation  is 
obtained  by  assumming  overall  trajectories  traversing  that  grid 
and  is  given  by  : 


where. 


m.  '  £  W  ?  °«Di 


(21) 


At  any  instant  of  time,  t,  the  $  field  of  the  variables 
is  assumed  and  the  predictions  for  a  time  increment  At  are  then 
obtained  by  solving  the  gas  flow  equations,  by  a  marching  integ¬ 
ration  algorithm.  In  this  way  the  solution  is  marched  forward 
in  time  until  the  desired  period  has  been  covered.  More  details 
on  the  solution  procedure  can  be  found  in  Ref53.  Once  the  solu¬ 
tion  of  the  gas  flow  equations  are  solved  at  any  time  then  the 
fuel  concentration  can  be  calculated  to  determine  the  spray  con¬ 
tours  as  follows  : 


The  mass  of  air  in  the  grid  at  this  moment  is  given  by  , 
Fig.  (  5  )  : 


and 


rn  =  r  d0  dr.t 
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(23) 


Some  results  of  the  calculations  of  the  fuel  droplet  trajecto¬ 
ries  in  swirl  and  open  combustion  chambers  are  given  in  Figs. 
(6)  and  (7)  respectively. 


It  is  seen  from  the  suggested  model  that  the  characteristics  of 
the  mixing  process  and  consequently  the  combustion  process  are 
mainly  a  function  of  the  droplet  size  distribution  .  Therefore, 
comprehensive  woil  is  required  to  obtain  fundamental  knowledge  for 
transient  characteristics  of  the  droplet  size  distribution  in  clear 
and  multifuel  . 


Fig(  4  ) 
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Radial  displacement  mm 


Rig(  7  )  Trajectory  of  a  droplet  inside  the  piston  bowl  at 
various  injection  pressures  . 
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CHAPTER  (, 


ATOMIZATION  OF  MULT I FUEL 

h-t  MECHANISM  OF  ATOMIZATIuN 

The  efficiency  of  combustion  systems  which  use  liquid 
fuels  depends  on  the  process  of  liquid  atomization  as  an  essen¬ 
tial  intermediate  stage  between  injection  and  combustion.  A 
spray  of  fuel  is  obtained  inside  the  combustion  chamber  either 
by  an  atomizer  or  an  injector.  The  injectors  introduce  the  fuel 
spray  into  the  combustion  chamber  with  high  velocity  creating 
suffient  turbulence  inside  the  fuel  spray  jet  to  satisfy  better 
dispersion  of  fuel  droplets.  The  atomizers  disintegrate  the 
fuel  either  internally  by  spreading  the  fuel  on  the  surface  of 
swirl  chamber  or  by  internal  blast  fluid  or  externally  by  ex¬ 
ternal  blast  fluid  and  then  disperse  the  resulting  droplets  in 
the  combustion  chamber  space.  Experimental  investigation  showed 
that  when  a  fuel  is  injected  through  a  smooth  circular  orifice 
into  a  chamber  of  quiescent  gas,  four  main  regimes  of  jet  disin- 


ii-  Atomization  of  fuel  with  higher  exit  velocity  and  aero¬ 
dynamic  resistance  force,  Fig.(8,b), 

iii-  Atomization  of  fuel  with  still  higher  exit  velocity  and 
initial-  turbulence  resulting  from  resisting  air  force,  Fig. 

(8 , c)  . 

iv-  Atomization  regime. 

These  four  regimes  of  breakup  are  encountered  as  the  injection 
velocity  is  progressively  increased.  The  mechanisms  controll¬ 
ing  the  breakup  of  liquid  jets  in  the  first  three  of  these  re¬ 
gimes  are  reasonably  well  understood.  However,  the  hreakup 
mechanism  for  the  jets  in  the  fourth  regime,  which  is  called 
atomization  regime,  is  currently  unknown.  Several  researchers 
have  proposed  mechanisms  ascribing  the  breakup  process  to  aero¬ 
dynamic  effects,  liquid  turbulence,  cavitation  phenomena,  ve¬ 
locity  profile  and  fuel  pressure  oscillations . 

On  the  bases  of  the  large  number  of  experiments  and  theore¬ 
tical  studies  carried  out  on  the  plain  jet  injectors  which  are 
dealt  with  in  this  work,  one  may  consider  that  the  atomization 
of  liquid  fuels  goes  through  the  following  stages  : 

1-  Spreading  of  fuel  jet  into  a  cone  with  a  cone  angle  depend¬ 
ing  on  the  fuel  velocity. 

2-  Entrainment  of  air  and  division  of  fuel  spray  into  concen¬ 
tric  conical  sprays,  as  shown  in  Fig.  (  9)  ,  with  various  velocities. 
The  central  spray  has  a  yelocity  greater  than  the  outer  one. 
Therefore,  radial  velocities  are  created  from  the  central  liga¬ 
ment  causing  the  deceleration  of  the  outer  ligaments. 

3-  Occurance  of  small  waves  in  the  liquid  surface  as  a  result 
of  internal  turbulence. 

4-  Tearing  of  the  surface  into  fine  ligaments  as  a  result  of 
drag,  surface  tension  and  viscous  forces. 

5-  Breaking-up  of  fine  ligaments  into  droplets  due  to  the 
above  mentioned  forces. 


6-  Further  breaking  up  of  these  droplets  by  secondary  ato¬ 
mization  . 

7-  Agglomoration  of  droplets  by  collision. 

The  relative  importance  of  these  six  stages  changes  according 
to  the  type  of  injector,  relative  velocity  between  the  air  and 
fuel,  fuel  properties,  air  pressure  and  temperature  and  air- 
fuel  ratio.  The  radial  velocity  creates  radial  turbulence  for¬ 
ming  ligaments  of  a  thickness  depending  upon  the  magnitude  of 
its  wave  length. 

The  theoretical  analysis  of  the  jet  disintegration  is 
based  on  Lord  Rayleigh's  studies"" who  considered  the  growth  of 
initially  small  disturbances.  He  showed  that  ideally  the  wave¬ 
length  which  governed  the  disintegration  was  that  for  which  the 
rate  of  growth  was  maximum.  He  found  that  the  required  wave 
length  was  given  by  : 


A  =  4.5  d  (24) 

J 

Thus,  If  the  volume  of  the  spherical  droplets  is  equal  to  the 
velocity  of  the  cylindrical  ligament  of  length  equal  to  the 
wave  length  and  radius  equal  to  the  jet  radius,  the  droplet 
diameter  is  given  by  : 

?'  D3  -  |d2  X 

or 

D  =  1.89  d  (25) 

In  practice,  Rayleigh's  breakup  of  a  plain  jet  occurs  only 
if  the  relative  yelpcity  is  low  and  the  surface  tension  is  ade¬ 
quate.  In  the  more  general  case,  drag  forces  tend  to  tear  the 
surface  of  the  turbulent  liquid  jet  into  fine  ligaments  of  wave 
length  A  •  These  ligaments  become  spheres  of  diameter  D  under 
the  action  of  surface  tension.  Hence,  the  drop  size  predicted 
by  Rayleigh  represents  the  upper  limit  of  possible  sizes. 


A  similar  analysis  of  breakage  of  a  liquid  column,  which 
takes  both  the  viscosity  and  surface  tension  into  account  has 
been  made  by  Weber .  The  time  required  for  a  column  of  diame¬ 
ter  dj  to  break  up  is  given  by  ; 

|  8  p  *5  3/2 

fc  *  C  |<— >  +6^c£  dJ  <26) 

!  t 

and  the  optimum  length  to  diameter  ratio  0  for  breakage  is 
also  given  by  : 
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N/2  n 


*5 


(27) 


However,  the  fuel  atomization  is  affected  by  the  spray  turbu¬ 
lence  in  the  fuel  nozzle,  which  depends  mainly  on  the  construc¬ 
tion  of  the  fuel  nozzle  and  Reynolds  number.  The  effect  of 

kinematic  viscosity  on  the  solid  spray  length  for  various  injec- 

5  7 

tion  pressures  was  determined  experimentally  and  the  results 
are  given  in  Fig.  (101.  The  solid  length  can  also  be  calculated 
using  the  following  formula  : 


i 

s 


372  d  M 


0.308 


-0.71 


-1.2. 


The  relation  of  length-diameter  ratio  of  the  diesel  engine 
nozzles  must  be  less  than  3/4  and  therefore  the  Reynolds  number 
may  not  characterize  the  spray  motion  completely. 

V 

York  et  al^®  had  developed  an  approach  to  study  the  pro¬ 
cess  of  atomization  in  which  they  indicated  that  the  ligaments 
are  subjected  to  internal  disturbances  which  finally  produce 
Rayleigh's  break-up  of  the  unstable  film  when  ; 


-*•  2 


Droplets  may  thereafter  breakup  during  their  motion  in¬ 
side  the  combustion  chamber,  which  is  defined  by  secondary  atomi 
zation,  as  a  result  of  drag,  inertia,  viscous  and  surface  ten¬ 
sion  forces.  The  aerodynamic  forces  imposed  on  the  droplets 


I 


Fig  (  10)  Effect  of  kinematic  viscosity  of  fuel  on  the  break 
length  for  various  injection  pressures  . 


during  its  motion  way  cause  .Jo:.' urination  of  thei  r  shape  and 
eventually  overcoming  the  droplet  cohesive  forces  and  may  in 
some  cases  break  t hem  i nt  o  smaller  fragments .  Rabine  et  al ^ 
showed  that  droplet,  hr  ear.  up  is  associated  with  a  critical  value 
of  the  Weber  number  -a f •  o i  which,  secondary  atomization  may  be 
expected.  Th i s  occurs  when  : 
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Critical  Weber  number  after  which  secondary  atomization  may 
be  expected  is  given  by  B'  rstein  et  al(’°  as  equal  to  100  and  by 
Lambrisis  et  a 1 ( ’ L  as  equal  to  65. 

Finally,  the  atomization  of  fuel  injected  in  diesel  engine 
combustors  depends  on  the  external  aerodynamic  resisting  force, 
surface  tension  force,  internal  binding  force  which  depends  on 
its  viscosity,  and  the  internal  turbulent  force  resulting  ftom 
the  initial  turbulence  created  during  the  fuel  motion.  The  ex¬ 
ternal  aerodynamic  resisting  force  increases  with  the  increase 
of  the  injected  fuel  velocity,  relative  velocity,  medium  den¬ 
sity  and  the  surface  area  of  the  spray.  The  complete  theore¬ 
tical  calculation  of  the  droplet  size  distribution  is  not  deve¬ 
loped  till  now  and  recourse  is  taken  to  experiments  for  its, deter¬ 
mination. 


Several  experimental  methods  have  been  used  for  measurement 

of  the  droplet  size.  These  methods  are  photographic  thermal, 

electric,  optical  methods  based  on  scattering  or  absorption  of 

liaht  as  well  as  physico-chemical  technique  and  slide  sampling 

technique.  A  review  of  the  published  techniques  for  drop  size 
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measurement  has  been  carried  out  by  Azzopardi 

McCreath  and  Beer6 3  pointed  out  that  for  sprays  from  liquid 
fuel  atomizers  in  furnaices,  drop  size  frequency  and  spatial  dis¬ 
tribution  control  fundamental  flame  characteristics  such  as 
radiant  heat  transfer,  flame  length,  flame  stabilisation,  smoke 


formation,  carbon  carry  over  and  formation  of  oxides  of  nitrogen. 
Droplet  size  distribution  is  required  also  for  the  calculation 
of  the  spray  concentration.  The  photographic  method  is  most 
suitable  for  the  study  of  a  single  droplet,  because  of  its  small 
depth  of  focus,  or  sprays  of  small  cone  angles.  Both  the  molten- 
wax  method  and  the  frozen  spray  method  are  not  suitable  for 
measuring  the  droplet  size  distribution  at  high  air  temperatures 
and  pressures,  since  the  progressive  development  of  the  fuel 
spray  inside  the  combustion  is  not  reflected  in  the  results. 

Amongst  the  other  techniques  it  seems  that  the  optical  ones 

show  most  promise.  If  an  overall  distribution  is  required  then 

the  diffraction  technique,  particularly  as  described  by 

6  4 

Swithenbank  et  al  appears  to  possess  admirable  sensitivity  to 
changes  in  both  the  peak  value  and  the  width  of  the  distribution. 

The  Laser-Doppler  approach  using  visibility^ ' 66  qs  still  limit¬ 
ed  to  very  small  particles  and  work  is  needed  to  extend  its 
range.  Since  the  drop  size  frequency  is  required  to  be 
determined  at  high  pressures  and  temperatures  simulating 
the  combustion  chamber  environment,  the  slide  sampling  method 
is  the  most  suitable  until  the  development  of  other  techniques 
reaches  a  stage  where  they  could  be  applied  freely.  In  this 
method  a  slide  with  thinly  coated  soot  deposited  from  a  candle 
flame  is  introduced  into  the  fuel  spray  for  a  short  time.  Arrange¬ 
ments  are  provided  to  prevent  interfering  of  the  spray,  to  satis¬ 
fy  suitable  moving  velocity  of  the  slide  preventing  the  soot  lay- 
er  distortion,  and  to  adopt  suitable  layer  thickness  minimizing 
the  flattening  coefficient  of  Crater  size.  The  impression  left 
by  the  impinging  droplets  in  the  soot  layer  are  photographed  un¬ 
der  microscope. 

The  experimental  set  up  is  shown  in  Fig.  (ii) .  It  consists 
of  a  pressure  chamber  of  diameter  150  mm  and  height  180  mm.  The 
pintle  type  injector,  Bosch  DNO  SD  211  with  a  nozzle  diameter  1  mm 
and  pin  diameter  0.55  mm,  is  fitted  on  the  top  of  the  chamber , Fig . (12 )  . 
The  fuel  is  pumped  to  the  injector  by  a  fuel  pump  driven  by  an 
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electrical  motor  using  a  gear  chain  system  with  a  reduction 
ratio  1:2.  The  rack  of  the  fuel  Mimo  ’  -- 
moving  system  to  determine  accor  d  >-  S')  ->  - 1  1 

plunger  and  consequently  the  quantity  of  deliveied  fm.  J  .  1  he 

fuel  is  fed  to  the  pump  from  a  reservoir  through  a  filter  and 
a  flow  meter  to  determine  the  rate  of  flow  of  fuel.  A  dummy 
injector  is  used  for  intermittent  operation  of  the  working  injec¬ 
tor  by  means  of  an  electromagentic  three  way  valve.  The  fuel 
pressure  at  the  inlet  of  the  injector  is  picked  by  a  piezoelec¬ 
tric  transducer  and  the  needle  lift  by  an  inductive  pick  up. 


Both  signals  together  with  a  signal  indicating  the  moving  moment 
of  the  slide  sampling  apparatus  are  displayed  on  an  osci 1 lograph . 
The  pressure  chamber,  which  is  fitted  with  a  pressure  manometer, 
is  fed  with  compressed  air  from  a  high  pressure  reservoir  through 
a  heater  and  a  pressure  regulating  valve.  The  flowing  air  is 
heated  by  a  methane  torch  capable  of  heating  the  incoming  air  up 
to  about  500  °C .  The  air  temperature  is  measured  by  a  Cu-Co  ther¬ 
mocouple  connected  to  a  potentiometer. 


A  special  sampling  apparatus  mounted  in  the  pressure  cham¬ 
ber  at  a  distance  60  mm  from  the  injector  is  used  for  measuring 
the  droplet  size,  Fig.  ( 1 3 )  .  It  consists  of  a  cylindrical  tube 
of  a  30  mm  outer  diameter  prepared  with  a  rim  slot  of  2.5  mm  width 
A  cylindrical  core  of  23  mm  outer  diameter  having  a  rim  slot  is 
mounted  inside  the  cylindrical  tube.  The  glass  slide  is  insert¬ 
ed  inside  it.  A  lever  ratchet  tension  spring  system  is ‘connect¬ 
ed  to  the  core  axis  to  rotate  it. 


Fuel  sprays  injected  from  the  injector  and  passing  from 
the  slots  of  the  fixed  cylindrical  tube  are  collected  on  the 
glass  slide,  when  the  slot  of  the  cylindrical  core  faces  that 
of  the  cylindrical  tube  during  its  rotation  as  shown  in  Fig.  ( 1 3 . t> ) 
The  rotation  of  the  cylindrical  core  is  satisfied  at  a  specified 
moment  after  the  injection  by  an  electrical  system.  Connecting 
this  circuit  by  the  main  switch,  the  electrical  current  passes 
to  both  electromagnetic  three  way  valve  and  a  time  relay.  The 
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magnetic  valve  changes  the  fuel  flow  from  the  dummy  to  the  active 
injector,  while  the  time  relay  lasts  the  current  of  the  solenoid 
valve  by  0.5  sec.  Such  delay  period  permits  _he  injector  to  in¬ 
ject  six  sprays  before  the  introduction  of  the  glass  slide  in 
the  spray  to  achieve  steady  state  conditions.  A  vibrating  switch 
which  operates  at  maximum  needle  lift  position  is  used  to  close 
the  solenoid  circuit  and  consequently  rotate  the  cylindrical  core 
only  at  this  position.  A  uniform  coot  layer  is  deposited  on  the 
glass  slide  from  a  kerosene  flame.  The  soot  layer  thickness  is 
determined  from  the  relation  between  the  droplet  Weber's  number 
and  the  impression  diametex^  in  the  order  of  1.5: 3. 5  of  maximum 
droplet  diameter.  The  soot  layer  thickness  is  taken  equal  to 
350  yum  for  diesel  fuel  and  120  ^m  for  gasolene.  The  average 
value  of  the  soot  layer  thickness  is  measured  by  a  travelling 
microscope. 

The  impinging  fuel  droplets  leave  impressions  on  the  soot 
layer.  These  impressions  are  magnified  60  times  under  an  elec¬ 
tronic  microscope,  where  representative  photographs  are  taken  at 
different  radial  distances  from  the  centre  of  the  spray.  These 
photographs  are  then  remagnified  ten  times  on  a  screen  by  a  pro¬ 
jector  to  satisfy  total  magnification  of  600.  The  number  and 
diameter  of  these  impressions  after  correction  from  Fig.  (14) 
are  determined  with  a  scale  graduated  in  0.5  mm  and  classified 
into  groups  with  size  intervals  10  yrfm.  Each  group  of  impressions 
represents  a  part  of  the  fuel  spray  at  a  specified  radial  dis¬ 
tance.  The  total  summation  of  these  groups  gives  the  total  drop¬ 
let  size  distribution.  The  accuracy  of  the  measurement  of  im¬ 
pression  is  about  +  3  pm.  Samples  of  these  droplet  photographs 
are  shown  in  Fig.  (15). 

The  same  technique  is  used  to  measure  the  droplet  size  dis¬ 
tribution  inside  an  operating  motored  diesel  engine.  Fuel  spray 
injected  from  the  nozzle  is  collected  on  a  slide  sampling  appara¬ 
tus,  which  is  fitted  in  the  space  between  the  two  discs  that  are 
fixed  over  the  piston  crown  Fig.(i6).  The  glass  slide,  which  has 


a  dimension  of  70x20  mm  is  allowed  to  rotate  around  a  pivot  fix¬ 
ed  in  the  piston  crown  under  the  effect  of  the  tension  spring 
whenever  the  stopper  moves.  The  stopper  is  allowed  to  move  un¬ 
der  the  effect  of  an  attraction  force  created  by  a  solenoid  valve. 
This  valve  can  be  actuated  at  a  difinite  controlled  angle  by  a 
push  button  which  is  closed  by  a  cam  fixed  on  the  pump  shaft. 

More  details  of  the  apparatus  and  results  are  given  in  Ref68. 

6 . 3  Physical  Properties  of  Multifuel 

The  main  physical  properties  which  affect  the  atomization 
process  are  viscosity,  surface  tension  and  density.  The  vis¬ 
cosity  of  fuel  affects  the  spray  form  together  with  the  atomi¬ 
zation  quality.  High  viscosity  causes  a  large  droplet'  size,i.e. 
low  and  poor  atomization,  and  a  high  spray  jet  penetration. 

Such  droplets  evaporate  with  a  slow  rate  forming  a  hetrogeneous 
mixture  with  imperfect  combustion.  Droplet  size  depends  also 
on  the  surface  tension  of  the  injected  fuel  as  the  surface  ten¬ 
sion  force  opposes  the  distortion  of  the  fuel  surface.  The  fuel 
density  affects  spray  penetration  as  the  kinetic  energy  of  fuel 
issuing  from  a  nozzle  is  one  of  the  main  factors  determining 
spray  behaviour. 

Studying  the  behaviour  of  the  previously  mentioned  proper¬ 
ties,  fuel  mixtures  of  diesel  fuel  and  another  fuel  as  light 
die‘sel  fuel,  kerosene  or  gasolene  are  tested.  Diesel  fuel  is 
used  as  a  hasic  fuel  and  the  other  fuel  is  blended  With  it  in 
intervals  of  10%  by  volume.  In  addition,  fuel  preheating  up  to 
70  °C  is  applied  to  all  these  tested  fuels  as  their  properties 
are  known  to  vary  considerably  with  temperature. 

Kinematic  viscosity,  surface  tension  and  density  of  the 

tested  fuels  vary  from  8.553x10  6  m3/S,  27.5x10  3  N/m  and  847  Kg/m 
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respectively  for  diesel  fuel  to  8.08x10  m  /S,  20.4x10  N/m 
and  732  Kg/ro3  for  gasolene  at  20  °C  as  shown  in  Figs. (17)  (18) 
and  Cl 9 ) -  It  is  clear  from  these  figures  that  the  viscosity  of 
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Fig(  )  Nature  of  kinematic  viscosity  variation  for  various 
blond  mixtures  . 


p,g(  19  )  Nature  of  density  variation  for  various  blend  mixtures 
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diesel  fuel  is  about  11  times  greater  than  that  of  gasolene 
while  the  surface  tension  and  density  are  about  1.4  and  1.2 
times  greater  than  that  of  gasolene.  The  kinematic  viscosity 
of  diesel  fuel  is  sharply  reduced  when  blended  with  gasolene  as 
compared  with  blending  with  kerosene  or  light  diesel  fuel. 
Similar  behaviour  can  be  deduced  for  both  surface  tension  and 
density.  It  is  found  also  that  the  physical  properties  of  the 
fuel  are  quite  sensitive  to  temperature.  The  variation  of  kine¬ 
matic  viscosity  and  density  with  temperature  of  the  tested  fuel 
is  shown  in  Figs .  (20&-21)  •  It  is  clear  from  the  results  that  the 
kinematic  viscosity  of  heavy  grade  fuels  such  as  diesel  fuel 
and  light  diesel  fuel  decreases  sharply  with  the  temperature 
increase.  A  lower  decreasing  rate  is  noticed  in  light  fuels 
such  as  gasolene  and  kerosene.  Correlations  of  the  experimental 
data  give  the  following  empirical  relations  for  the  kinematic 
viscosity,  surface  tension  and  density  : 

=  A  (T  +  30)  “b  ecF 

al  =  (1-9x10_4T)  (axF  +  KjJlO-3  (30) 

P£  =  (1-9x10_4T) (K2+a2F) 


where  the  constants  are  given  in  table  2. 

Table  2:  Constants  of  Equation  30 


diesel  fuel-  diesel  fuel-  diesel  fuel-  . 
light  diesel  kerosene  gasolere 

fuel  _ 


A 

5 . 55xlO~4 

1.66x10" 

al 

0.51 

3.0 

a2 

16.0 

56.0 

b 

1.14 

1.09 

c 

0 . 274 

1.266 

K1 

27.4 

24.9 

K2 

846.0 

806.0 

6.15x10 

7.1 
17.0 

1.1 

2.33 

20.8 

745.0 
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Variation  of  fuel  viscosity  with  temperature  rise 
for  clear  and  multifuel  . 


The  atomization  quality  of  fuel  sprays  depends  to  a  great 
extent  upon  fuel  properties.  The  atomization  of  fuel  sprays  is¬ 
sued  from  injectors  has  been  studied  by  many  authors.  Most  of 
the  previous  experiments  were  carried  out  to  determine  the  Sauter 
mean  diameter  in  atmospheric  conditions.  In  reciprocating  en¬ 
gines,  however,  fuels  are  injected  in  high  compressed  air.  More¬ 
over,  the  droplet  size  distribution  of  fuel  injected  in  such  com¬ 
pressed  air  is  required  for  the  computational  model  of  fuel-air 
flow  field.  Some  authors  ^  'reported  the  droplet  size  distribu¬ 
tion  under  the  same  conditions  of  diesel  engines.  The  liquid 
immersion  technique  consisting  of  collecting  a  sample  of  spray 
in  a  cell  filled  with  immersion  liquid,  was  employed.  The  cell 
was  located  at  a  high  distance  650  mm  from  the  injector  to  prevent 
break  up  of  large  droplets  on  the  cell  and  flow  of  small  droplets 
with  the  air  stream  resulting  from  the  motion  of  the  slit  plate. 
However,  the  droplet  size  distribution  in  this  case,  can  not  re¬ 
present  that  in  the  vicinity  of  the  injector  which  is  required 
for  the  computational  model  due  to  secondary  atomization  and  re¬ 
combination  of  droplets.  Moreover,  the  droplet  size  might  not 
be  always  constant  during  the  injection  period  but  fluctuates 
rapidly  during  the  injection  period  due  to  the  variation  of  the 
injection  pressure  and  nozzle  needle  lift.  The  transient  charac¬ 
teristics  of  the  droplet  size  distribution  was  reported  by 
Hiroyasi  and  Kodata  using  the  liquid  impersion  technique.  The 
cell  was  located  at  high  distance  160  mm  from  the  injector.  More¬ 
over,  the  used  immersion  liquid  can  be  used  only  for  kerosene. 

Therefore,  the  objective  of  the  present  work  is  to  obtain 
the  droplet  size  distribution  of  multifuel  for  various  injection 
pressures,  back  pressures,  and  sampling  positions.  Some  examples 
of  the  droplet  size  distribution  are  exhibited  in  Figs. (22)  to 
(25).  These  are  the  frequency  distribution  diagrams  of  droplet 
number.  More  details  of  the  results  can  be  found  in  '  .  The_ 

variation  of  the  frequency  distribution  diagram  of  droplet  number 
for  various  fuels  is  shown  in  Fig. (22)  for  a  pintle  type  injector,Figi2) 
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Fig  (  22  )  Droplet  size  distribution  of  jlear  fuel  sprays  injected 
by  a  pintle  type  injector  with  injection  pressure  118 
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The  samples  are  taken  during  the  period  from  0  to  0.5  ms  from 
the  start  of  injection.  It  can  be  noticed  that  the  curve  has  a 
sharp  and  high  peak  at  small  diameter  droplets  for  gasolene 
fuel.  The  curve  has  rather  gentle  slope  for  diesel  fuel  and 
the  curves  of  kerosene  and  light  diesel  fuel  are  in  between. 

This  means  that  large  droplet  diameters  are  obtained  with  diesel 
fuel.  The  droplet  size  distribution  of  sprays  obtained  b^  a 
pintle  type  injector  for  various  diesel-gasolene  mixtures  is 
given  in  Fig.  (2  3.)  .  it  can  be  noticed  that  the  curves  acquire  a 
sharp  and  high  peak  from  the  basic  distribution  of  pure  diesel 
fuel  at  the  small  size  range  of  droplet  diameters  with  increas¬ 
ing  the  gasolene  percentage  in  the  mixture.  This  behaviour  may 
be  attributed  to  the  decrease  of  viscosity,  surface  tension  and 
density  with  the  increase  of  gasolene  percentage  in  fuel. 

The  effect  of  medium  pressure  and  injection  pressure  on 
the  frequency  diagram  of  the  number  of  droplet  of  multifuel 
sprays  is  studied.  The  distribution  curves  are  found  to  be  shif¬ 
ted  towards  the  large  size  range  with  the  increase  of  medium 
pressure  and  decrease  of  injection  pressure  as  shown  in  Figs. 

(24  &  25  )  • 

The  influence  of  various  factors  such  as  injection  pressure, 
nozzle  diameter,  fuel  quantity,  injection  angle  and  distance  from 
the.  injector  is  studied  also  inside  the  cold  engine  model.  Some 
examples  of  the  frequency  distribution  diagrams  are  shown  in 
Figs.  (26  )  and  (2?)*  The  fuel  droplets  are  collected  on  the  slide 
at  20  degree  b  TDC.  The  results  show  that  the  size  distribution 
c 1  ve  is  shifted  towards  the  small  size  of  droplet  diameters  as 
the  nozzle  diameter  decreases  and  also  as  the  fuel  quantity  in¬ 
creases.  This  is  believed  to  be  due  to  the  formation  of  small 
ligament  with  small  nozzle  diameter. 

>3  ^  Mean  Diameter  of  Spray  Droplets 

The  droplet  size  distribution  is  evaluated  by  the  relevant 
mean  diameter  D  defined  by 
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Fig (  23  )  Droplet  size  distribution  of  diesel-gasolene  mixtures 


7oSn/n  °/oSn/n 


Percentage  No  of  droplets  %  dn/n 


Droplet  dvn'.eU 


Fig (  2  r’  ) 


Drop  l °t  s  i  ,-e  d  j  st.r  i  but  i '  n 
inside  a  eo  1  d  enqinc  model  i 


8  16  24  32  40  48  56  64  72  ft) 

Droplet  diameter  D  /am 


Fig(  27)  Droplet  size  distribution  for  diesel  fuel  sprays 
inside  *  cold  engine  model  for  various  nozzle 
diameters . 


max 


D<1i  <■> 


l/lq-p) 


,D 


D„ 


min 


qp 


D 


max 


°Pi§  ao 


D  . 
nun 


(31) 


The  mean  diameter  of  spray  droplets  is  used  to  describe  the 
quality  of  atomization,  trajectory  and  penetration  of  the  spray 
and  to  facilitate  correlation  with  physical  parameters.  There 
are  at  least  six  mean  diameters  which  may  be  used  to  repesent. 

the  spray  distribution  and  therefore,  the  values  of  p  and  q  de- 
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pend  on  the  field  of  application  .  Consider  the  Langmuir's 
equation,  which  is  commonly  used  for  describing  the  evaporation 
rate  of  spherical  particles  in  still  air,  a  necessary  prerequisite 
is  to  known  the  relevant  average  diameter  of  the  fuel  spray  drop¬ 
lets.  For  number  of  droplets  having  different  diameters  the  equa¬ 
tion  is 
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is  the  linear  mean  diameter. 
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Also  considering  the  efficiency  of  an  atomizing  operation, 
in  which  a  fraction  of  energy  "E"  lost  in  pressure  drop  a  is 
used  in  forming  new  ligament,  the  equation  may  be  expressed  for 
the  entire  spray  as  follows  : 
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which  can  be  arranged  in  the  following  form  ; 
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where,  is  the  Sauter  mean  diameter. 

It  is  seen  thus  that  the  evaporation  law  holds  for  and  the 

atomizing  efficiency  for  D^- 


The  Sauter  mean  diameter  is  thus  used  for  the  mass 

transfer  and  combustion  problems  since  the  surface  area  and  the 
•olume  of  droplet  are  of  primary  importance  for  evaporation  and 
ombustion.  The  Sauter  mean  diameter  can  be  defined  as  the  dia¬ 
meter  of  a  droplet  that  has  the  same  surface  to  volume  ratio  as 
that  of  the  total  spray.  It  is  calculated  from  the  formula 

D32.  =  ^  /  Z  n^D?  .  (35) 


The  influence  of  various  factors  on  the  Sauter  mean  dia¬ 
meter  for  various  multifuel  spray  is  studied,  sprays  of  the  four 
main  types,  diesel,  light  diesel,  kerosene  and  gasolene  are  test¬ 
ed  with  the  same  injection  pressure.  The  rack  of  the  fuel  pump 
is  adjusted  to  its  maximum  position.  Some  of  these  results  are 
shown  in  Figs. (28  “30  ).  It  is  clear  from  Fig. (28)  that  the 
mean  droplet  diameter  decreases  slightly  with  the  increase  of 
light  diesel  fuel  quantity.  The  rate  of  decrease  of  the  mean 
droplet  diameter  increases  with  the  increase  of  kerosene  and 
gasolene  quantity.  This  may  be  attributed  to  the  decrease  of 
viscosity  and  surface  tension  of  the  mixture.  It  is  found  also 
that  the  mean  droplet  diameter  decreases  with  the  increase  of 
fuel  pressure  due  to  the  increase  of  the  available  atomization 
energy  affecting  the  initial  disturbance  and  entrainment  as  shown 
in  Fig. (29) .  The  effect  of  medium  density  is  illustrated  in 
Fig. (30)  which  indicates  that  the  Sauter  mean  diameter  increases 
with  increasing  the  medium  density.  This  is  believed  to  be  due 
to  the  decrease  of  the  entrained  air,  and  consequently  increase 
of  ligament  diameter,  resulting  from  the  increase  of  air  vis¬ 
cosity.  It  can  also  be  attributed  to  the  agglomeration  of  drop¬ 
lets  due  to  small  penetration  of  spray  in  a  denser  medium.  It 
is  clear  also  from  the  figure  that  the  increase  of  Sauter  mean 
diameter  with  the  increase  of  medium  density  in  case  of  gasolene 
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Fin(  28)  Variation  of  Sauter  mean  diameter  with  the  quantity  of 
y  diesel  fuel  in  multifuel  at  various  injection  pressures 
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is  higher  than  in  the  case  of  diesel  fuel.  This  may  be  due  to 
further  droplet  agglomeration  because  of  the  sharp  loss  of  the 
kinetic  energy  >f  droplets  leading  to  faster  retardation  of 
droplets. 


Studying  the  variation  of  Sauter  mean  diameter  across  the 
fuel  spray  it  was  found  that  the  local  Sauter  mean  diameter  does 
not  change  in  case  of  gasolene  while  in  case  of  diesel  fuel  it 
decreases  at  the  periphery  of  spray.  The  relation  between  the 
Sauter  mean  diameter  and  other  mean  diameters  is  shown  in  Fig. 
(31) .  As  regard  to  the  maximum  droplet  diameter  for  multifuel 
sprays,  a  general  expression  is  deduced  in  the  following  form  : 
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Generalization  of  Droplet  Size  Distribution  Data 

Since  there  is  no  theoretical  analysis  until  now  for  the 
calculation  of  the  droplet  size  distribution  of  fuel  sprays,  se¬ 
veral  probability  expressions  have  heen  developed  by  several  re¬ 
searchers  for  this  purpose.  At  the  same  time  the  droplet  size 
distribution  curye  is  greatly  affected  by  fuel  properties,  in¬ 
jection  pressure  and  operating  conditions.  To  generalize  such 
different  data  it  is  required  to  normalize  the  droplet  diameter 
by  dividing  it  by  a  suitable  mean  droplet  diameter  to  transfer 
these  discrete  values  to  a  unique  curve.  In  the  present  work, 
the  mean  droplet  diameter  D^2  is  found  to  be  suitable  for  this 
purpose.  Figure  (33)  shows  such  representation  of  cumulative 
frequency  distribution  of  dimensionless  droplet  diameter. 

In  order  to  choose  the  most  suitable  distribution  function 
which,  fits  adequetly  the  experimental  data,  curve  fitting  of  the 
experimental  data  with  these  previously  suggested  probability 
functions  are  carried  out.  Curye  fitting  of  the  experimental 
data  using  the  simple  Rosin-Ramraler  distribution  function  J  gives 
the  following  expression  : 
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Fig(  32)  Comulative  frquency  distribution  of  droplet  size  for 
multifuel  sprays  with  dimensionless  droplet  diameter 
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-0.427/(D/D32)5,32 
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and  using  log  probability  distribution  function  gives  the 
following  expression  : 


4-n-  -  0.94  e~  (1  •  657y  +  0. 5 )  2  8) 

dy 

where , 

y  =  In  ( D/D  ^ 2 ) 

7  2 

and  using  upper  limit  log  probability  distribution  function 
gives  the  following  expression  : 


^  ,  0.572  e-<1.015y-0.913)2 

dy 

where , 

!  0.76  (D/D32) 

y  “  ln  j  1.99-  (D/D32)' 


(39) 


It  has  been  found69  that  the  droplet  size  in  diesel  fuel 
sprays  was  well  correlated  with  the  chi-square  distribution  func¬ 
tion  when  using  the  cumulative  distribution  of  volumes  instead 
of  numbers.  An  attempt  is  made  to  correlate  the  size  of  the 
droplet  sampled  with  the  chi-square  distribution  and  to  find 
the  degree  of  freedom  as  shown  in  Fig.  (33)  •  It  is  clear  from 
the  figure  that  the  degree  of  freedom  f=9  which  was  given  by 
Hiroyasi  cannot  fit  the  data  while  f  =  17  fits  adequatly  the 
experimental  data.  The  smaller  diameters  obtained  by  Hiroyasi 
may  be  attributed  to  the  location  of  immersion  cell  at  a  distance 
of  160  mm  from  the  nozzle.  The  chi-distribution  function  which 
fits  the  experimental  data  can  be.-  written  in  the  following  form  : 


7.b  -7 . 5 (D/D, ^ ) 

^=1953  (D/D32)  e  d(D/D32)  (40) 

An  attend  t  to  correlate  the  droplet  sampled  with  the 
Tanasawa-Tesima  distribution  function"^  is  also  carried  out. 
This  droplet  size  distribution  function  can  be  written  in  the 
following  form  for  the  normalized  curve) 
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which  can  be  written' in  the  following  straight  line  form  : 


dn/n  d(D/D32) 


(D/D32)a 


-  In  a  -  b(D/D32)e 


plotting  (1/ (D/D32 ) a )  dn/n  d{D/D32)  against  (D/D32)^  for  dif¬ 
ferent  assumed  values  of  a  and  B  one  can  get  the  most  suitable 
values  of  a  and  6  which  satisfy  the  minimum  square  deviation. 
Determining  the  values  of  a  and  Q  from  the  experimental  values 
it  is  possible  to  write  the  equation  in  the  following-  form  : 


ar,  9  4  -2.4(D/D  } 

(nd (P/d32)  *  10 (D/D^^)  a  <43, 

It  is  required  now  to  find  out  the  most  suitable  function 
from  these  five  expressions  which  fit  the  experimental  data- 
adquatly  and  allow  for  easy  use  and  extrapolation.  A  comparison 
of  these  expressions  with  experimental  data  is  given  in  Fig.  (34  )■ 
This  comparison  she  vs  that  both  the  T-T  and  chi-square  equations 
give  a  sufficiently  accurate  fit  of  the  data  while  the  R-R  equa¬ 
tion  shows  poor  fitting  of  the  data.  The  upper  limit  log  pro¬ 
bability  equation  shows  an • acceptable  fitting  of  data  but  it  has 
a  complicated  form  which  does  not  permit  direct  evaluation  of 
the  droplet  size  distribution.  Regarding  the  log  probability 
equation  it  is  found  that  it  fits  well  the  largest  droplets 
range  while  greater  deviations  are  obtained  in  the  small  diame¬ 
ter  range.  Therefore,  it  is  possible  to  conclude  that  both  the 
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T-T  and  chi-square  distribution  functions  satisfy  adequatly  the 
experimental  data  of  either  clear  or  multifuel  sprays  injected 
through  pintle  nozzle  injector  within  the  accuracy  of  the  used 
experimental  technique. 

6.7  Generalization  of  Mean  Droplet  Diameter 

The  mean  droplet  diameter  of  certain  fuel  spray  depends 
on  a  large  number  of  yariables  such  as  fuel  properties,  opera¬ 
ting  conditions  and  geometrical  dimensions  of  the  fuel  nozzle. 
Dimensionless  analysis  facilitates  experimental  determination 
of  mean  droplet  diameter  at  various  conditions.  The  dimension¬ 
less  groups  decrease  the  number  of  experiments  required  to  obtain 
an  empirical  expression.  Dimensional  analysis  gives  the  follow¬ 
ing  relation  for  Sauter  mean  diameter  : 

D32/d  =  A  Rea  Web  C®  (p  f/  ^)f(*n/d)g 


where,  A,a,b,d,e,  f  and  g  are  constants  to  be  determined  from 

the  experimental  results.  The  indices  that  fit  the  experimental 

data  by  a  smooth  curve  were  determined  and  the  general  equation 

of  Sauter  mean  diameter  is  derived  by  plotting  Re~°‘434We”1'°^(j  _1"° 
-0  12 

(p  f/  f^)  ’  versus  D-^/d  as  shown  in  Fig.  (35).  The  following 

relation  is  determined  for  the  Sauter  mean  diameter 


D32/d  =  107  Re 
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An  explicit  form  of  this  relation  with  yarious  parameters  is 
found  in  the  following  form  : 

-O  54 

o  o  r  a  n  a  n  n  a  v  » 


0.385  0.737  0.737  0.06 


D32  =  6156  v 


A  satisfactory  coincidence  is  obtained  when  the  results 
determined  by  this  relation  is  compared  with  the  experimental  re¬ 
sults  as  shown  in  Fig.  (36)* 
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Fig ( 36  )  Comparison  of  calculated  and  experimental  values  of 
Sauter  mean  diameter  . 


CHAPTER  7 
CONCLUSIONS 


Theoretical  calculation  of  the  fuel  spray  behaviour  inside  nwltifuel 
engine  combustion  chambers  which  leads  to  better  evaluation  of  heat  re¬ 
lease  requires  complete  knowledge  of  the  air  flow  field.  In  the  same 
time  theoretical  calculation  of  the  flow  field  is  very  complicated  due 
to  the  presence  of  transient,  compressible  and  turbulent  flow.  So  a 
mathematical  model  supported  by  experimental  results  capable  to  inves¬ 
tigate  clearly  the  flow  field  inside  diesel  engine  swirl  chamber,  and 
describing  the  effect  of  the  working  conditions  and  constructional 
parameters  on  it  is  required. 


Development  of  a  complete  mathematical  model  for  the  prediction  of 
turbulent  flow  field, in  swirl  chambers  has  been  carried  out  in  the 
present  work  by  considering  geometrical  and  operating  conditions.  The 
mathematical  model  of  the  turbulence  flow  has  been  obtained  by  solv¬ 
ing  the  finite-difference  form  of  the  equations  of  energy,  mass,  momen¬ 
tum,  the  turbulence  kinetic  energy  K  and  its  dissipation  rate  e 
taken  the  compressiblity  effects  in  consideration. 


The  source  term  for  the  turbulence  dissipation  rate  e  has  been  modi¬ 
fied  to  suit  the  geometrical  and  operating  conditions  as  follows  : 

S£  =  |  (C1(1-C3R£)  G  -  C2  pc} 

Where  the  coefficient  F  has  been  determined  experimentally  from  the 
local  shear  stress  and  the  turbulence  generation  as  follows  : 


.C 


1 


(t/p  K)  =  CD  F 
o 


,2  n2 

A  n 

a  r* 

°e  D 


An  experimental  set-up  has  been  built  up  to  facilitate  the  measurements 
of  air  velocity  components,  and  turbulence  fluctuating  velocities  and  to 
investigate  the  effect  of  various  constructional  and  working  conditions 
on  the  flow  field,  a  careful  comparison  of  the  experimental  results 
over  a  wide  range  with  the  predicted  values  has  been  carried  out  to 
obtain  reliable  information  about  the  coefficients  turbulence  model. 
Velocity  components  have  been  measured  by  a  temperature  compensated 
hot-wire  anemometer.  In  addition  a  microphone  condenser  pick  up  has 
been  developed  for  the  measurement  of  velocity  components  inside  the 
swirl  chamber. 


Based  on  the  experimental  and  predicted  results  carried  out  in  this 
work  the  major  contributions  can  be  summarized  as  follows  : 

1-  The  hot-wire  anemometer  is  sensible  to  the  medium  pressure  and 
therefore  calibration  of  anemometer  at  various  pressures  has  been  done. 
A  formula  for  the  relation  between  velocity,  pressure  and  anemometer 
responce  has  been  obtained  in  the  following  form  : 
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PU  =  32  (V  -  VO)1'8 

2-  The  velocity  in  swirl  chamber  of  diesel  engine  can  be  adequately 
determined  by  the  microphone-condenser  after  it  has  been  calibrated 

at  various  pressuresand  temperatures. 

3-  A  mathematical  model  has  been  developed  to  determine  the  flow 
field  in  swirl  chambers  taking  into  consideration  the  effect  of  cons¬ 
tructional  and  operating  conditions  upon  the  local  turbulence  value. 

4-  The  coefficient  CD  determining  the  local  turbulence  value  has 
been  determined  experimentally  and  it  is  affected  by  the  radial  dis¬ 
tance,  engine  speed,  port  area  and  swirl  volume.  Relations  with  these 
parameters  are  obtained  as  follows  : 

CD  =  1.16  (r/Rs)-0'08 

CD  *  1.2255  (N/Nc)-0'88 

C  =  33.359  (A  ^/A) 0,704 

U  J  o  229 

CD  =  1.303  (Vs/Vc)U'^y 

and  a  general  dimensionless  correction  factor  is  obtained  in  the 
following  form  : 

F=33.47(r/RS)~°'08{N/Nc)"0,88 (Aj /A) 0,704 (Vs/Vc) 0,229 

5-  The  effect  of  streamline  wall  curvature  and  centripetal  accele¬ 
rations  has  been  introduced  to  correct  the  length  scale  of  turbulence 
determining  e  -equation,  by  multipli.ng  the  coef f icient^C^  with  the 
Richardson  number  function  (1-CjRf ) ,  the  con? 1  ant  C3  has  been  obtained 
equal  to  0.04  from  the  comparison  of  experimental  and  predicted  re¬ 
sults  . 

6-  The  suggested  theoretical  model  gives  satisfactory  results  than 
the  previous  model  and  gives  also  logic  values  for  the  length  scale  of 
turbulence. 

7-  The  effect  of  various  constructional  and  operating  conditions 
has  been  investigated  and  the  following  results  are  obtained  j 

a)  The  port  area  has  the  most  significant  effect'  on  the  tangential 
velocity  and  turbulence  intensity  than  other  parameters. 

b)  Two  eddies  are  noticed  at  both  sides  of  the  inlet  port.  The  in¬ 
tensity  of  these  eddies  increases  during  the  compression  stroke,  and 
intensif ing  with  the  decrease  of  the  port  area  ratio  (Aj/A) ,  increase 
of  engine  speed  and  increase  of  swirl  chamber  volume  ratio  (Vg/Vc). 

c)  The  radial  velocity  has  smaller  values  than  the  tangential  ve¬ 
locity. 

d)  The  turbulence  mixing  length  l  for  a  wide  range  of  results  at 
various  constructional  and  operating  conditions  can  be*  expressed  as 
follows  t 

A/R  s  ■  const 


This  constant  takes  values  varies  from  0.1  to  0.2. 
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8-  a  spray  model  inside  the  combustion  chambers  of  diesel  en¬ 
gines  is  suggested.  The  model  is  based  on  the  combination  of 
the  discrete  droplet  model  with  a  multidimensional  gas  flow  pre¬ 
diction  method.  The  variation  of  droplet  diameter  and  velocity 
and  consequently  the  variation  of  heat  transfer  to  droplet  and 
drag  force  are  taken  into  consideration.  However,  this  model 
requires  the  determination  of  the  droplet  size  distribution  and 
mean  droplet  diameter'. 

The  droplet  size  distribution  and  mean  droplet  diameter 
have  been  determined  experimentally  for  clear  fuel  as  well  as 
for  multifuel  using  the  slide  sampling  technique.  Droplet  size 
distribution  functions  and  Sauter  mean  diameter  data  are  obtain¬ 
ed.  The  use  of  these  relations  demonstrated  substantially  better 
results  of  the  spray  model.  The  model  is  capable  of  determining 
the  spray  behaviour  for  any  fuel. 

In  the  past,  research  efforts  have  not  emphazised  the  ef¬ 
fects  caused  hy  changes  in  the  fuel.  Many  new  efforts  have  been 
started  to  determine  the  evaporation  of  clear  and  multifuel  under 
various  operating  conditions,  and  the  multifuel  droplet  behaviour. 
The  variation  of  droplet  velocity  at  different  introduction  loca¬ 
tions  is  also  required.  The  application  of  these  new  efforts  in 
spray  model  will  emphasize  the  effect  of  the  fuel-  character  on 
the  heat  release  pattern  and  facilitates  successful  prediction 
of  fuel  concentration  required  for  the  calculation  of  the  concen¬ 
trations  of  nitric  oxide  and  soot  in  diesel  engines. 
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